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Synopsis 


The Magnetohydrodynamic (MHD) model describes plasmas as a collisional, electrically 
neutral, and a highly conducting fluid. In turbulent flows, fluctuations are present over a wide 
range of length scales (wavenumbers). The nonlinear terms in the MHD equations transfer 
energ}' between the wavenumbers of different sizes. Energy is transferred between velocity 
(u) modes, between magnetic (b) modes, and between velocity and magnetic modes. In 
this thesis we investigate these energy transfers and test the predictions of various MHD 
turbulence phenomenologies. Our study is restricted to two-dimensional MHD turbulence. 

The energy spectrum in fluid turbulence is given by E{k) = C < e where 

< e > is the average rate of energy dissipation and O' is a universal constant. In MHD 
turbulence two classes of phenomenologies have been proposed: (1) The phenomenology 
by Kraichnan (1965), Iroshnikov (1963), Dobrowolony et ai.(1980) [1-3] (together referred 
as KID) suggested that the energy spectra of (= u ± b) is E'^ik) = . 

where Bq is the mean magnetic field, and are the energy flu.xes of z*. .\dditionally. 
it was predicted that energy fluxes H^ are equal. An extension of KID phenomenology 
predicts that H'^/H" = m'^ lm~ [4], where are the scaling e.xponent of the energy spectra 
E^(k) (2) .A.ccording to a different class of phenomenology [5-7], the energy spectra 

in -MHD turbulence is Kolmogorov-like with E^{k) = and the 

constants C* are believed to be non-universal. These phenomenologies have been tested in 
simulations [8-13] and Solar Wind observations [14-16] and appear to favour Kolmogorov- 
like phenomenology, but there is still a lack of consensus among researchers. One of th^ 
objectives of this thesis is to investigate this issue further. 

The interaction in MHD invloves three modes k, p, q satisfying the condition k d- p -f q 1 
The combined energy transfer from any two modes to the third mode is known. In this thesii 
we develop a formalism to describe energy traxisfer between a pair of modes by the mediatioj 



ii 

of the third mode in the triad. We have used this formalism to study some important features 
of energ}' transfer. 

The properties of energy transfer are an important component of turbulence studies. The 
total energ\' is known to cascade from low”- wavenumbers to the high wavenumbers. Some other 
aspects of energy transfer in MHD turbulence have also been studied, such as the nonlocal 
transfer between the high wavenumbers and the low wavenumbers [17-19] and the net energy 
transfer to a mode from all other modes [20, 21]. The predictions of some of these studies have 
not been rigorously tested. Additionally, several other interesting energy transfers have not 
been investigated as yet. In this thesis we have numerically studied some of these transfers, 
for example, the energy flux out of a velocity wavenumber-sphere, a magnetic wavenumber- 
sphere and various energy fluxes from velocity to the magnetic modes. We have also studied 
the fine-grained energy transfer between various velocity and magnetic wavenumber shells. 

In the discussion given below we briefly describe the contents of each chapter, which 
includes its aim and its central results. 

Chapter 1 contains a survey of the field. We also state the thesis problem here. 

In Chapter 2 we test the predictions of KID phenomenologies and the Kolmogorov-like 
phenomenology in direct numerical simulations of forced 2-D MHD turbulence. The phe- 
nomenologies are tested over a wide range of normalised cross helicity (normalised velocity- 
magnetic field correlations). To test the phenomenologies we plot the time-averaged energy 
spectra and the energy cascade rates. It is diflncult to distinguish between the spectral expo- 
nents —3/2 and —5/3 predicted by KID phenomenology and Kolmogorov-like phenomenology 
respectively. Instead, following the approach of Verma d al. [11], w'e test the predictions 
of the energy- cascade rates H^ by the KID and the Kolmogorov-like phenomenology. Our 
results of the energy Cascade rates are inconsistent with the KID phenomenology and its 
extension b}' Grappin, but are in close agreement with the predictions of Kolmbgorov-like 
phenomenology with having a small variation with the normalised cross-helicity Using 
our simulation data we study the dependence of Kolmogorov constants on cr^. 


In Chapter 3 we find a new approach for studying energy transfer in turbulence. We 
obtain a description of the energy transfer between a pair of modes in a triad by the mediation 
of the third mode. We find that the mode-to-mode energy transfer can be expressed as a 
sum of an "effective’ transfer and a "circulating’ transfer. These two have the property that 
the effective transfer causes a change in the modal energy, while the circulating transfer 
does not affect modal energy since the same amount of energy is transferred from k to p. 
from p to q, and from q back to k. W'e have derived the expression for effective transfer 
between a pair of velocity modes, a pair of magnetic modes, and between a velocity and a 
magnetic mode. In each of these, the effective transfer between the two modes is mediated 
by the third mode of the triad. By applying the constraints of rotational invariance, galilean 
invariance, and finiteness, we show that the circulating transfer is equal to zero upto the 
linear order in a series. Using the expression for the effective transfer between modes, we 
define several energy cascade rates and shell-to-shell transfer rates (energy transfer between 
two wavenumber shells). 

In Chapter 4, the energy cascade rates and the shell-to-shell transfer rates defined in 
Chapter .3 are computed by a direct numerical simulation of 2-D MHD turbulence. The 
cascade rates and shell-to-shell transfer rates are averaged over a quasi-steady state achieved 
by large-scale kinetic energy forcing. The picture of energy fluxes that emerges is quite; 
complex — there is a forward cascade of magnetic energy, an inverse cascade of kinetic 
energy, a flux of energy from the kinetic to the magnetic field, and a reverse flux which 
transfers the energy back to the kinetic from the magnetic. The energy transfer between 
different wave number shells is also complex — local and nonlocal transfers possess opposing 
features — for example, energy transfer between some wave number shells occurs from kinetic 
to magnetic, while this transfer is reversed for some other shells. The net transfer of energy 
is from kinetic to magnetic. The results obtained from the flux studies and the shell-to-sheiJ 
energy transfer studies are consistent with each other. Comparison of the simulation results 
with earlier analytical studies [18, 19] and numerical simulations [20] are also made. ! 



IV 


The evolution of global quantities are presumed, by most researchers, to depend only 
on the initial value of the global quantities and their spectra. In Chapter 5 we numerically 
stud}" the effects of subtle changes in initial conditions, keeping the global quantities and their 
spectra unaltered, on the evolution of global quantities in two-dimensional MHD turbulence. 
We find that a change in the initial phases of complex Fourier modes of while keeping 
the initial values of total energy, cross helicity, and Alfven ratio (ratio of kinetic energy to 
magnetic energy) unchanged, has a significant effect on the evolution of cross helicity. On 
the other hand, the total energy and Alfven ratio are insensitive to the initial phases. 

The thesis is concluded in Chapter 6 by a brief smnmary and the future research 
directions. 
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Chapter 1 
Introduction 


Turbulent flows are ubiquitous. They occur over various scales, from the scale of a labor- 
atory to that of a cluster of galaxies, spanning a range of the order of 10^®. Turbulence in 
non-conducting fluids can be generated and studied in a variety of e.xperimental setups in 
the laboratory, and also observational studies of turbulence can be carried out in naturally 
occurring flows such as the atmosphere. The Reynolds numbers required for generating MHD 
turbulence however are not usually achie\'able in a terrestrial laboratory. MHD turbulence 
exists in a variety of astrophysical flows such as in the Solar Wind, interstellar medium, 
stars, etc. The study of turbulence is important for these systems. For example, the molecu- 
lar clouds form 20% — 50% of the interstellar gas in the milky way, and they are the sites of 
star formation. Turbulence in the molecular clouds is believed to play an important role in 
star formation [22]. In the Solar Wind, turbulence plays an important role in the evolution 
of energy, temperature, pressure, etc. [23]. Turbulence is believed to play an important role 
in the generation of large-scale magnetic field in the galaxy and other astrophysical objects : 
[24]. 

.A complete description of any of the flows discussed above (for example, the Solar! 
Wind) requires inclusion of features like compressibility, inhomogeneity, shear, anisotropy, 
etc. However, in this thesis we will concern ourselves with the simplest model — the homo- 
genous. incompressible MHD equation without any complications like shear and anisotropy.' 



1.1 Problem definition 
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and study some of the generic properties at the small scales. This mode! .^how.s a rich hcha- 
viour which will be useful for understanding more realistic flows. 

1.1 Problem definition 

One of the most celebrated results in fluid turbulence is due to Kohm)goruv( 1!I41 ! [‘jr);. He 
predicted that the energ}' spectrum in fluid turbulence is proportional to k * . hor Mill) 
turbulence. Kraichnan(1965) [1] and Iroshnikov(1963) [2] gave the first plienoiuenuktcy and 
predicted both the kinetic and magnetic energy spectrum to be proportional to Doliro- 

wolny et ai [3] generalized Kraichnan-Iroshnikov arguments and still predicted a spec- 

trum Recently, several researchers have argued for a Kolmogorov-like energy sped rum 
for MHD turbulence as well [5-7, 11, 13, 26, 27]. The Solar wind observations {14- 16]. and 
numerical simulations [11-13] appear to support Kolmogorov-like phenomenology rather than 
the phenomenology of Kraichnan-Iroshnikov-Dobrowolny (KID). In this thesis we have rlone 
an extensive numerical simulation to resolve the above controversy. To this end, we calculate 
both the energy spectrum and the energy cascade rates. Our simulation i.s based on the 
pseudo-spectral method and done in two-dimensions. 

The interaction between wavenumbers in turbulence involves triads. In this thesis we 
have deri^•ed a very useful formula for calculating energ}- transfer rate from one mode of the 
triad to another, with the mediation of the third mode. This has been done for both fluids 
and MHD Using these formulae we have defined various energy fluxes and energy transfer 
rates between 'shells’ These quantities give information about the energy transfer between 
different scales of the velocit}- field, the magnetic field, and between the scales of the velocity 
and the magnetic field. M'e have numerically computed these quantities from our simulation 
data. Our ntunerical results shed important light into the amplification of magnetic energy 
and locaJit}' of interactions. 

Many researchers have studied the evolution of total energy and Cross-Helicity in terms of 
the initial Alfven ratio {f u^dx/ fh^dx), and normalised cross-helicity {2 J {u.h)dx/ J{u^ -t- 
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b^)c^x. where u and b are the velocity and the magnetic fields respectively. There exists 
another quantity that can be varied keeping the initial Alfven ratio and normalised cross- 
helicity fixed. This quantity is the Fourier ‘phase’. We have studied the evolution of global 
quantities in terms of these quantities and found that the evolution of normalised cross-helicity 
sometimes depend on the initial phases, though the evolution of energy does not. This result 
is discussed towards the end of the thesis. 

-A.fter the above brief discussion of the problem definition of our thesis, we turn to a 
description of some of the definitions and basic results of MHD turbulence which are related 
to our problem. 

1.2 MHD equations 

MHD is a fluid model for plasmas. In this model the plasma is collisional. electrically 
neutral, and highly conducting [28]. The simplest MHD flows are incompressible. There are 
two independent variables in the incompressible MHD equations — the velocity field U. and 
the magnetic field B. The MHD equations governing the evolution of U and B are 


^ + (U.V)U = --Vp + ^(B.V)B + :.V2U. (1.1) 

at p Airp 

— + (U.V)B = (B.V)U + ^V2B. (1.2) 

at 

V.U = 0 (1.3) 

V.B = 0 (1.4) 


where u is the kinematic viscosity, p. is the magnetic diffusivity (= c^/(47rcr) with a being the 
conductivity of the plasma and c being the velocity of light), p is the mass density, and p is the 
pressure which is a sum of the thermal pressure pth and the magnetic pressure, pt = (Srr).: 

Eq. (1.3) is the equation of continuity for an incompressible fluid. For incompressible flows. 
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the pressure adjusts to maintain a constant mass density, and hcnrc is not at. n iit .»< jjdf'nt 
quantity. The pressure p can be determined from the poisson s ecjuation 

= V.(B.V)B - V.(U.V)U -1.51 

In the equations (1.1)-(1.2), the viscous and resistive terms i/V^U and p V=®B dissipast- enerp\. 

The MHD equations (1.1)-(1.4) are often written in a non-dimeusioiial fonn. If liu- 
characteristic length scale is £, the characteristic velocity scale is W. the charact<*risiic .scak* 
of the magnetic field is B, and the characteristic density scale is po. then the variables in 
equations (1.1)-(1.4) can be written in a non-dimensional form using U' = I lU. B' = B B^^. 
f = UtjC, p' = p/po, and p' = pJ^poW). The Eqs. (1.1)-(1.4) written in terms of these 
dimensionless variables are 


dV' 

dt' 

+ (U'.V')U' = -^Vp -f ^(B'.V')B' -f ^ 
p' UL 

n.(ii 


^ -h (U'.V')B' = (B'.VOU' + ^ V'^B', 

(1.7| 


V'.U' = 0 

(l.S) 


V'.B' = 0 

(1-9) 


where Ca = Bq/ y/Airpo and /? = The parameters W£/z/ and UCj p are the mechanical 

and the magnetic Reynolds numbers, respectivel}^ For geometrically similar situations with 
identical Reynolds number and initial conditions the dimensionless variables U'(r'. t') and 
B'(r',t') are identical. At large Reynolds numbers (> 10^ usually) the flow is turbulent. 

In a turbulent flow U and B have random fluctuations over a wide range of spatial 
and temporal scales. It is useful to separate the fluctuations into a ‘mean’ (Uo,Bo) and a 
•fluctuating’ part (u, b), i.e., U = Uq -h u and B = Bo + b. The mean quantities can vary 
over long time scales and length scales. In the simplest situation the mean can be assumed to 
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be constant in space and time. A constant Uq can be eliminated from the MHD equations by 
shifting to a frame of reference moving with a velocity Uq. However, a constant Bo cannot be 
eliminated from the MHD equations by a change of reference frame, and it, therefore, plays 
a d\'namically significant role. In Section 1.4 we will describe some of the phenomenologies 
which attempt to describe the effect of Bq on the statistics of u, b. In a reference frame of 
the mean flow (Uq = 0) the fluctuations u, b satisfy the equations 

^ - ^(Bo.V)b + (u.V)u = -lvp+ -^(b.V)b + r/V"u, (1.10) 

Ot iirp p Aitp 

^ - (Bo.V)u + (u.V)b = (b.V)u + ^V^b, (1.11) 

V.u=0 (1.12) 

V.b = 0 (1.13) 

The magnetic field can be expressed in Alfvenic units where we make the replacement 
b/ \/4~p b. The MHD equations in terms of the new u and b are 

(1.14) 

(1.15) 

(1.16) : 

V.b = 0 (1-lf) : 


^ - (CA.V)b + (u.V)u = --Vp + (b.V)b + lyV^u, 
ut p 


dh 

dt 


- (Ca.V)u + (u.V)b = (b. V)u + 


V.u = 0 


where b is in Alfvenic units in the above equation, and Ca = Bo/\/47rp 
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In terms of the Els^ser variables = u ± b, the MHD equations can he px[>re.-.ni as 

— q: (C^.V)z=^ 4- fz^.V)z=^ = --Vp + i/+V^Z=^ 

dt ^ ^ P 

V.z=^ = 0 

where = {u ± fj,)/2. The z^*" are the normal modes of the incompressihkx in\-i«ci(i f - U. 
pt = 0) MHD equations; z"*" propagates anti-parallel to Bq. and z" propacate'- parallel tu 
Bo with the propagation speed of Ca- The nonlinear terms in the MHD equations vaiiish if 
either z'^ or z~ is zero, i.e., if the system has pure modes. A turbulent flow is a mixture of 
z"*" and z~. The nonlinearity modifies the Alfvenic fluctuations z^ and creates tur!mlenct> at 
high Reynolds numbers. 

1.2.1 Important global quantities in MHD 

The inviscid MHD equations possess some conserved global quantities. These roiis<‘rvet! 
quantities are 

1 . the total energy per unit mass 

2 the velocity-magnetic field correlations (cross-helicity) 

H, = J {u.h)dx = ^ J {{z+f -{z-fjdx-, ( 1 . 21 ) 

3. the magnetic helicity in absence of Bq 


! l.lhl 



( 1 . 22 ) 


where a is the magnetic vector potential, i.e.,b = V x a. hf; = + 2Bo.A„ is conserved 
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in presence of Bq [15]. In 2-D, Hm is equal to zero, but the mean-square magnetic vector 
potential is conserved (in absence of Bq) 

A = Ja^dx. (1.23) 

The z'^ energy. = (1/2) f(z'^)^dx and the z~ energy, E~ = (1/2) /(z~)^dx are also 
conserved quantities but they can be expressed in terms of B and He- Note that the kinetic 
energy B^ = (1/2) f u'^dx and the magnetic energy Bb = (1/2) j b'^dx are not conserved. In 
the MHD equations (1.1)- (1.2) the term (U.V)U conserves the kinetic energy and (U.V)B 
conserves the magnetic energy. The terms (B.V)B in Eq. (1.1) and (B.V)U in Eq. (1.2) 
exchange energy between the velocity and the magnetic field by stretching the magnetic field 
lines (see reference [24]). The following dimensionless quantities are often studied in literature 
: the normalised cross helicity defined as (Tc = 2HclB = (B'^ — B~)I(B'^ + B~), and the 
Alfven ratio defined as = Bu/Bb. 

The evolution of the global quantities defined above has been studied in simulations [S. 
9, 29-34]. The fields u and b usually show a tendency to align with each other in numerical 
simulations, i.e., ctc increases [8, 9, 29, 30] — this is called dynamic alignment. However, the 
tendency for dynamic alignment is weaker for low values of <Jc and large Reynolds number 
[8]. The Alfven ratio is found to decrease indicating the preferential decay of the magnetic 
energy compared to the kinetic energy — this is known as selective decay. Ting et al. [31] 
performed an extensive study over a range of initial values of and CTc and found distinct 
regimes of evolution. However, it has been believed so far that the evolution depends only 
on the initial value of the global quantities and their spectra. However, in Chapter 5 we will 
show that the evolution can be affected by some subtle features of the initial state. The values 
of ac and can have an influence on the energy spectra (defined below) — this is a subject: 
of investigation of Chapter 2. 
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1.2.2 Spectral analysis in MHD turbulence 

Spectral analysis is very important in turbulence research. The velocity iirld uix t ran be 
expressed in terms of the Fourier components defined as 

where d is the real space dimension. The relationship between the modal c*iit*ig\ in 

real space and energy density £*(1^) in the Fourier space is given by 


- f u^{x)dx= f E^{k)dK 

2 Juntivol 


1 1 . 25 ) 


where £’*‘(k) = (l/2)|u(k)p. In case of isotropic turbulence, the directionally inteprated 


E'^ik) 


E'^ik)^ [ £l“(k)dk, 

Jlk\=.k 


(1.26) 


is a useful quantity. Similarly, the magnetic energy spectrum can be defined a.s 


E\k)= [ £l^(k)dk, (1.27) 

2|k|=A: 

where £’^(k) = ( 1/2) jb(k)P, and the energy spectrum of the inward and theoutw'ard propagat- 
ing Alfven waves, i.e., fluctuations are defined as 


= f (1.28) 

where E^{k.) = (l/2)|z^(k)p. In Chapter 2, we will study the energy spectra £'^(k) using 
data from our numerical simulations. 
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The evolution equation for the kinetic energy spectrum E'^{k) is obtained from the MHD 
equations (1.14) 

+ Y. :jS~(k|p,q)+ x: is-^Cklp.q) (1.-29) 

k+p-|-q=0 k-f p'fq=0 ^ 

where the nonlinear term 5““(kjp, q) is 

5““(klp, q) = -Re (i[k.u(q)][u{k).u(p)] + z[k.u(p)][u(k).u(q)]) { 1.30) 

‘S’“^k|p,q) = -i 2 e( 2 [k.b(q)][u(k).b(p)] + f[k.b(p)][u(k).b(q)]) (1.31) 

The interaction involves a triad k, p, q satisfying the condition k + p + q = 0. The quantity 
S^^lklp. q) is the combined energy transfer from the modes p and q to mode k [35. 36]. The 

energv’ transfer individually from the modes p to k and from mode q to k are not known. 

In this thesis we will formulate an approach for calculating the energy transfer between two 
modes of the triad with the third mode mediating this transfer (see Chapter 3). 

The nonlinear terms transfer energy between the various modes but conserve the overall 
energ}'. The loss of energy from a sphere in Fourier space to the modes outside is termed as 
the energy cascade rate (or energy flux)C The study of energy cascades due to each of the 
nonlinear terms of Eqs. (1.14)-(1.15) is a major objective of this thesis. We shall discuss the 
known features of energy cciscades later in this Introduction. 

1.3 Fluid turbulence 

Man}’ current theories in turbulence have their foundation in the phenomenology proposed by 
Kolmogorov [25] and Obukhov [37] in the context of fluid turbulence. Other recent ad^•ances 
in fluid turbulence have also made an impact on the study of MHD turbulence. We will 
briefly discuss some of the aspects of fluid turbulence that are relevant to our work. : 

Hhe terms ‘energy flux’ and ‘energy cascade rate’ are synonymous and will be used interchangeably. 
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1.3.1 Scaling laws in fluid turbulence 

Kolmogorov proposed a set of hypothesis [25] for fluid turbulence using uhiih tiic eneig\ 
spectra is obtained. The h 3 ''pothesis can be stated as follows . 

1. The energy flows from the energy-containing low wavenumbers to the energ> ing 

high wave numbers. There exists a range of wavenumbers, intermediate to the (>nerg\- 
containing wavenumbers (L), and the energy-dissipating wavenumber.- iliat sta^s 

in statistical equilibrium — it is called the ‘inertial range’. 

2. In the inertial range fluctuations of approximately same length I ~ {/.■“’• participate 
in nonlinear energy transfer, i.e., energy transfer is local in wavenumber spac(‘. It 
also follows that the inertial range statistics are independent of statistics of the energ\'- 
containing low wavenumbers and the energy-dissipating high wavenumbers. 

3. The statistical properties in the inertial range are isotropic. 

4. In the inertial range, the energy flux !!(/:) is independent of the waveimmbei and erpial 
to the dissipation rate, < e >. The energy spectrum E(k) in the inertial range* dt'petids 
only on the mean energy dissipation rate < c> and the wavenumber k. 

Lnder the above set of h\*pothesis, the energy spectrum in the inertial range is obtained 
using dimensional analysis. Since [£’(/c)] = [< e >] = and [A-] = cm~K it 

follows that 

£(/c) = c < e >2/3 ifc-5/3 

where C is the universal Kolmogorov constant. According to the Kolmogorov's phenomeno- 
logy, the energy spectrum in the above equation is universal, i.e., E{k) has the same scaling 
exponent in all types of turbulent flows with arbitrary geometries. In experiments (see the 
references in [38]) no significant departure from the E{k) ~ fc-s/a has been found. 
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A tacit assumption in Kolmogorov’s phenomenology is that the dissipation rate has a 
uniform distribution in space. This is reflected in the use of mean dissipation rate < e > 
in Eq. (1.32). However, the dissipation rate has large spatial and temporal fluctuations — 
a phenomenon called ‘intermittency’. The scaling exponent of E{k) which is equal to —5/3 
in Eq. (1.32) is expected to get modified due to intermittency. Intermittency is frequently 
studied by means of the longitudinal structure functions of order p defined as 

‘5'p(r) = (|[u(x + r) - u(x)].fr) . (1.33) 

where the term ‘longitudinal’ implies that the velocity component along the line joining the 
two points is being considered. According to Kolmogorov’s phenomenology, at inertial range 
separations, l/k^ « r « L, the structure functions Sp have a power law behaviour 5^(r) ~ 
with Cp = p/3. The spectral exponent m in E{k) ~ k~^ is related to C 2 as m = C 2 + 1 (see 
[38]). Experiments show that Cp deviates from p/3, and the deviation increases with p (see 
Erisch [38]), an effect that is believed to be due to intermittency. Several models have been 
proposed to model the intermittency in fluid turbulence (refer to Frisch [38] for a discussion 
of these models). 

The elementary interactions in Navier-Stokes and MHD equations involve a wavenumber 
triad k, p, q with k + p + q = 0. Researchers have investigated the following question: given 
a wavenumber k, which triads contribute most to the interactions ? (see Domaradzki [39]). 
One implicit assumption in most models is that the interactions are local, i.e.. the triads 
with /l- ~ p ~ g contribute most to the energy transfer to wavenumber k. The nonlocal 
triad interactions are those where k p » q. Locality has been investigated in several 
direct numerical simulations [39-43]. One of the conclusions from simulations is that the; 
contribution of the nonlocal triads (A: ~ p >> q) to the interactions is not small, but in such' 
triads the dominant exchange of energy occurs between wavenumbers of similar magnitudes., 
i.e., between k and p [39, 40]. We will point out in Chapter 3 that the arguments used in 
these papers are ambiguous. We derive an expression for the 'effective energy transfer' that; 
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eliminates the ambiguity. In Chapter^ we have investigated the localit} 


• assumption u.suig 


our new formalism. 

Note that in Kolmogorov's phenomenology, the energy was assumed to flow from the 
low to the high wavenumbers. The kinetic energy is a conserved quantit} in imisciti three- 
dimensional Navier-Stokes equation. The two-dimensional Navier-Stokes equation has an 
additional conserved quantity — the enstrophy where a; = V x u is the* vortiritw 

The consequence of vorticity conservation is important. Assume that encrg}' i.s be*ing fed at 
wavenumber kj. Using Absolute equilibrium theories (see [44] for a discussion of these theor- 
ies). Kraichnan [45] argued that there is an inverse cascade of kinetic energ\- to wa\-(*nun!!>ers 
less th an kj, and there is a forward cascade of enstrophy to wavenumbers greater than Ay. 
These predictions were thereafter confirmed in numerical simulations and experiment- isee 
Frisch [38] and Nelkin [46] for the references). Kolmogorov-type analysis applied to en.-trojihy 
cascade leads to a k~^ in the wavenumber region kj « k << kd^ while the usual analysis 
applied to the energy cascade yields a k~^^^ spectra in the wavenumber region k << kj 
[45. 47]. The k~^^^ spectrum has been observed for small wavenumbers in direct numerical 
simulations and also in recent experiments (see Frisch [38] and Nelkin [46] for references). In 
summary, the cascade of conserved quantities is a crucial component in the phenoinenolog}’ 
of turbulence. For further discussion refer to Frisch [38] Lesieur [44], Kraichnan [48]: and 
Nelkin [46] provides a useful guide to references. 


1.4 MHD turbulence : Phenomenology, Simulations, Solar 
Wind observations, and Analytical results 

In this section we will discuss the MHD turbulence phenomenologies proposed so far and 
current status of simulation, analytical, and observational studies. For a review on this 
subject refer to Verma and Dar [27]. 
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1.4.1 Energy spectrum: phenomenological results 

The effect of the mean magnetic field on the energy spectra has been a matter of considerable 
debate. Broadly, there exist two view-points described below. 

Kraichnan [1] and Iroshnikov [49] suggested that the mean magnetic field, or the magnetic 
field of the large-scale eddies, affect the energy transfer to the small scales, and that there 
is an equipartition of kinetic and magnetic energy at small-scales, and they have a k~^'^ 
energy spectra. These arguments were restricted to the <7^ = 0 condition. Dobrowoln\' 
et al. [3] generalized the arguments to arbitrary <7^ and again found E^{k) ~ k~'^^~. The 
phenomenology proposed by Kraichnan, Iroshnikov. and Dobrowolny et al. has been referred 
to as the KID phenomenology in our thesis. A different view-point was proposed by Marsch 
[5]. Matthaeus and Zhou [6], Zhou and Matthaeus [7] according to which the energy spectra 
k — this is the Kolmogorov-like phenomenology. 

For any mode with length scale I ~ k~^, we can define three time scales in the MHD 
equations: the Alfven time scale ta = {CAk)~^, and the two nonlinear time scales = 
{kz^)~^ . The two nonlinear time scales are not equal due to the structure of the nonlinear 
interactions. The KID and the Kolmogorov-like phenomenology differ in the choice of the 
interaction time scale. 

Recall from Section 1.2. that and z~ are two independent modes traveling anti- 
parallel and parallel to the mean magnetic field. In a turbulent flow z"*" and z~ interact 
via the nonlinear term (z^.V)z^ leading to a exchange of energy between the modes. In 
the phenomenologies the interactions are assumed to be local in wavenumber space. Let 
denote the time scale over which the interacting eddies z'^{k) and z~{k) remain correlated. 
The variation Sz'^lk) in the amplitudes of the eddies over the time interval is 

(1.34); 

Assuming random interactions, the amplitude variation in N such interactions will be 

Az^ ^ y/N{5z^). (1.35) ! 
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Therefore, the number of interactions N* taken to obtain a variation cinal ,o ,ts 


amplitude zf is 


A’ 


'± 






and the corresponding time T* — is 








The time scale for energy transfer is assumed to be r±. The energy fluxes Ii=^ of 


fluctuations, out of a sphere in k-space, can be estimated as 


T± 


and being oppositely traveling fluctuations, KID phenomenology assumes that 
they will get de-correlated in an Alfven time scale. Hence n ~ ~ {('a^tK Since 

{zff ~ kE^{k), we get 

n+ ~ n- w ^E+ik)E-{k)k^, ( l.:i9l 

Bo 

for arbitrary (Tc value. For cTc = 0 (£''■ ~ E~), we get 

E^{k)K{BoUY^^k-^^\ (I.-IO) 


It is assumed that in the absence of mean magnetic field, the magnetic field of the largest 
eddy would play the role of Bq. 

In the Kolmogorov-like phenomenology the nonlinear time scale, r^i are taken to be the 
interaction time-scale for the c* eddies, i.e., r* ~ ~ (fcz^)“T Then, from Eq. f 1.38) the 

fluxes are 

n* =» (rrf f n? )i. (1.41) 


The above equation can be inverted to give 

E^{k) = - 5 / 3 _ 


(1.42) 
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Hence the spectral exponent is the same as the Kolmogorov’s —5/3 scaling exponent in fluid 
turbulence. This phenomenology was given by Marsch [5]. Matthaeus and Zhou [6], Zhou 
and Matthaeus [7]. In contrast to the Kolmogorov constant in fluid turbulence, the constants 
C’*’ are expected to depend on the cross helicity and the Alfven ratio [50]. 

The KID phenomenology predicts that the fluxes H"^ and H" are equal, whereas in rhe 
Kolmogorov-like phenomenology they are unequal. In order to construct a model within the 
KID phenomenology where n"*" ^ 11“ Grappin tt al. [4] proposed an extension to the KID 
phenomenology. They assume that the in Eq. (1.37) is to be determined by including 
non-locality in the analysis. The nonlinear time scale if interactions are 

local. Assuming non-locality, Grappin et al. replaced rv£ in Eq. (1.37) by an average 
(tvl) — ! /q/A: obtained rri^ -k m~ = 3 and H'^/H" = m'^/Tn~. where 

E^{k) ~ A:-'"*. 

The Kolmogorov-like phenomenology is also a limiting case of a more generalized phe- 
nomenology constructed by Matthaeus ajad Zhou [6], and Zhou and Matthaeus [7], which 
yields 

Bo + ^/kE^k) 

where .\ is a constant. Here the wavenumbers satisfying ykE^{k) » Bo follow —5 3 
spectrum and the wavenumbers iJkE^{k) « Bq follow a —3/2 spectrum. 

In both 2-D and 3-D MHD turbulence, there is a forward cascade of energy. The con- 
sequence of this result is that the phenomenologies discussed above for MHD turbulence 
based on energy cascade are valid for both 2-D and 3-D flows. Recall that in the earlier 
discussion of Section 1.3 we pointed out that 2-D fluid turbulence is fundamentally different 
from 3-D fluid turbulence. 
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1.4.2 Simulation results 


Earlier theoretical [4] and numerical [8-10] work in this area clainu'd sup|»ort fm KID 
phenomenology and its extension to Grappin. The relationship rir. v.-ith 

YJI+ fYi~ = 3 was found to be present in numerical calculations of the Edd_\ Dtiinju'ct Quasi- 
Normal Markovian approximation (EDQNM) closure scheme [4] (an introdurtivui !<< LDQN'M 
can be found in the book by Leslie [35] and Lesieur [44]). It was reported in nuniericai sim- 
ulations of Biskamp and Welter [ 8 ] that Eik) ~ for low cr, values [ 8 ]. The e.xpt.ufut'^ 
T 77 + and 777 “ were found to be unequal and the relationship m+ ■+ = 3 was fuund to hold 

appro.ximately in numerical simulations [9, 10]. However, the errors involved in these simula- 
tions were too large to reach a definite conclusion. A clearer picture has been mneruiim from 
the work done during the last decade. 

The KID phenomenology and the Kolmogorov-like phenomenology differ in their pre*dio- 
tions regarding the scaling exponents of energy spectra and also in the relationship ht'tween 
the fluxes and H". While H''' = H" in the KID phenomenology [see Eq. (1.39)]. Il’*' ^ 11“ 
for Kolmogorov-like phenomenology’ [see Eq. (1-41)]. Verma ei al. [ 11 ] exploited this dis- 
tinction between the two phenomenologies to make a comparison between them. In their 
simulation they found that 7 ^ 11 “. thus raising a doubt on the validity of KID }>h<‘n(unen- 
ology. In Kolmogorov-like phenomenology, at low cTc values the constants CW and C~ in 


Eqs. (1.42) are expected to be equal and E'^IE~ = (n+/n“)^. These equalities were found 
to hold in Verma et al. ’s simulations [ 11 ] of decaying 2 -D MHD turbulence. Grappin et al.'s 
[4] extension of KID phenomenology’, where H+ZH” = was not tested by Verma ef 

al. in their simulations. Verma et al s simulations were done for decaying turbulence and the 
energy’^ spectra and energy cascade rates were given without averaging. In our simulations we 
will average the fluxes and spectra over a steady-state which will reduce the errors, and we 
will also test Grappin et al.'s predictions. 

These earlier simulations were not done at high enough Reynolds numbers. Con- 
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sequently, the inertial range obtained was not large enough for an accurate estimation of 
the spectral exponents. Much higher Reynolds numbers have been achieved in some of the 
recent numerical simulations done by Politano et ai [12], and Muller and Biskamp [13]. In 
addition, the phenomenon of Extended Self Similarity (ESS), which we describe below, is 
now used by researchers to extract the exponents with greater accuracy. 

In the inertial range, the structure functions Sf[r) = follow a power law 

behaviour Sf[r) ~ r'^p, where 5z^{r) = (z*(x + r) — z^(x))). By plotting 5*(r) against 
5*(r). a power law ~ was found to e.xist over a much wider range [51]. This is 

called E.xtended Self Similarity (ESS) and it was earlier reported in fluid turbulence [52-55]. 
ESS has been found in simulations of MHD turbulence [12. 51] and Solar Wind observations 
[56, 57]. In fluid turbulence, (3 = 1 is an exact result [38. 58], so the ESS plot gives the 
e.xponents Cp directly. In MHD turbulence (3 is not equal to one [12. 13]. So, only the ratio 
the exponents Cp/Cs can be obtained by plotting Sf{r) against S*(r). 

Fortunately can be calculated using mixed correlators L^{r) =< {Sz^[r))^Sz'^{r) > 
for which an exact result is known. Politano et al. [59, 60] derived the following exact 
relationship, 

< [Sz'^YSz^ > —2 < z^(x)r*(z)r^(x + r) >= —Cdt^r (1-44) 

which can be simplified for large magnetic field to 

< {6z^Y5z^ >= —Cdt^r (1-45) 

The above relationship is analogous to ]du(r)p ~ r obtained from the Navier-Stokes equation 
for homogenous, isotropic fluid turbulence. The relationship (1.45) has been verified in a 
numerical simulation [12], A power law 5^(r) ~ [.^'*(f)]^p was found to extend well beyond 
the inertial range [12], and it was used to obtain the exponents C*- 

The exponent in E^{k) ~ A;”'"* and the exponent (2 cf second order structure 
functions 5* are related to each other by = (* + 1. Using the relationship between 
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and rr?*, it can be seen that (2 = 1/2 according to the KID pheiunuf'snijf > 1 ') bnt cnual to 

2/3 in the KoImogorov-Iike phenomenology. In addition, Cr = P 

phenomenology and (p = p/4 for the KID phenomenology. Hence, it nia\ 'uv ijj^j 

C 3 = 1 for Kolmogorov-like and (4 = 1 for KID phenomenology, {‘.sing KSS the . \jH(iie;;t>. 
were obtained for 2-D MHD turbulence by Politano ci al. [12] with gre:it<'r accuracy tiian 
was possible earlier. They found that = 0.69 (rn'*' = 1.(591 and <,T {j.7 =: ],7j 

[ 12 ]. Therefore, these results support the Kolmogorov-like {thehonK'nelo*.'’. . which predirt? 
that 772^ = 5/3 [Eq. (1-42)]. when we consider 

Politano ef al. also observed that (4 ~ 1 , leading them to conclude that KID phetioinen- 
ology is a valid description of MHD turbulence. However, we differ in onr vitnv r rear dim; the 
implication of their results for two reasons: ( 1 ) KID phenomenology i.s not consistent with 
the Verma et a/, s observation that n+ ^ H* and £+/£“ ~ (H^/n-/* which are ron.si.stent 
with the Kolmogorov-like phenomenology; ( 2 ) (^2 found in Politano fi als siimilation i.*' clo.se 
to the value expected by Kolmogorov-like phenomenology and the ob.served dc\ iation is small 
{Sm+ ~ 0.04). A deviation of this magnitude can be expected to arise due to interniiifenry 
coriections. In this light, the result (4 ~ 1 may be a fortuitous occurance arising In' a 
deviation from the Kolmogorov value of (4 = 4/3 due to intermittency. 

In another recent simulation of 3-D MHD turbulence performed by Miiller and Biskamp 

[13], It was found that m± ~ 5/3. They concluded that their results .supjiort the validity of 
Kolmogorov-like phenomenology to MHD turbulence. 

MHD turbulence is found to be intermittent. Numerical simulations [12, 1 3. 5 1 ] and Solar 
Wmd observations [61-63) show that the exponents C" are nonlinear funct.ons of p. and the 
deviation front both p/3 and p/4 increases with p. In simulations, the probability distributions 
of [(fb(r)l = |b(x-f r) -b(x)| is found to have a tail [51, 64], which is longer for smaller 
scale separations r, [64] indicating small-scale intermittency. The probability distribution 
I u(r)l ju(x -hr) - u(x)| has a smaller tail and appears to be less intermittent than 
[<fb(r)|. Some models have been proposed to account for intermittency on the lines of models 
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earlier proposed for fluid turbulence, but all these models are under the framework of KID 
phenomenology. The predictions of She-Leveque (SL) model [65] for fluid turbulence compare 
well with the exponents obtained in 3-D simulations of MHD turbulence [55]. However, 
the same predictions do not agree with the exponents The (p’s predicted by both the S-L 
model for fluid turbulence and its modification to MHD turbulence [66, 67] do not have a 
close agreement with the Cp’s found in 2-D MHD turbulence simulations. It is clear that none 
of the present phenomenologies can explain the intermittent behaviour of MHD turbulence. 
In this thesis, we will not investigate intermittency. Rather, we will focus on resolving 
the controversy between the applicability of Kolmogorov-like and KID phenomenologies by 
testing their predictions for the energy spectra and the cascade rates. 

1.4.3 Solar wind observations 

The Solar wind provides a unique opportunity for in situ measurements of MHD turbulence 
properties. Researchers have used observational data to verify various turbulence theories. 
From the arguments of Matthaeus and Zhou [6], Zhou and Matthaeus [7], and Marsch [5], we 
see that the arguments of KID phenomenology are applicable when Bq» yJkE^{k) (when 
k is in the inertial range); on the other hand Kolmogorov-like phenomenology is applicable 
when Bo « yJkE^{k). In the Solar Wind, Bq >> ykE^iJ^ for k in the inertial range, 
and hence KID phenomenology should be applicable according to these authors. However, 
the energy spectra in the solar wind are found to be nearly proportional to k~^^^ [15. 16, 14], 
The Fig. 1.1, taken from reference [15], shows the total energy spectra to be proportional to 
^,-i.69±o.08 AU^. Hence, the Solar Wind observations show that there is some inconsistency , 
in the phenomenological arguments given above. Verma [26] has attempted to resolve this 
contradiction by applying renormalization group treatment to MHD equations. 

The Solar Wind is found to be intermittent [61-63], The usual signatures of inter-, 
mittency, namely, spectral exponents having a nonlinear dependence on p and a tail in the 

-An astronomical unit (AU) is the mean distance between the Sun and the Earth which is approximately; 
equal to 1.5 x 10* Km 
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Figure 1.1: Total energy spectra E{k) = E'^{k) + E {k) in the Solar wind at 1 AT (.shown 
by the solid line) taken from Matthcieus and Goldstein [15]. 

probabilit}' distribution of |^v(r)| are found in Solar Wind data. The SL model for MUD is 
found to agree quite well with the Solar Wind data [61]. However. Grauer [66] has cautioned 
that the errors in the spectral exponents from the Solar Wind are too large to be certain 
about the agreement between SL model and the Solar Wind data. 

We will discuss below the results of the analytical theories applied to MHD turbulence. 

1.4.4 Analytical results 

The Renormalization group (RG) approach (see McComb [68] for an introduction) was ap- 
plied to MHD turbulence by Fournier et al [69], Camargo and Tasso [70], Verma [26]. 
Fournier et aL found that RG of MHD is quite complex. There are various types of regimes 
where the transport coefficients, viscosity i/ and resistivity /x, get renormalized either due 
to kinetic small scales or due to magnetic small scales. They also find a strong dimensional 
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dependence. Later, Camargo and Tasso [70] applied similar techniques, and showed that the 
effective (renormalized) viscosity and resistivity are, in general, non-negative. They mapped 
the effective transport parameters as a function of initial parameter. 

Verma [26] constructed a self-consistent renormalization group scheme. Within this 
scheme the mean magnetic field was found to get renormalized, scaling as while the 

magnetic energy scaled as The Renormalization group (RG) calculation demonstrated 

that the magnetic field acting on the eddies I w k~^ is not the magnetic field of the largest 
eddies but a renormalized field. This gives rise to a k~^l^ spectrum instead of the 
spectrum. These calculations were done by replacing the term Ic.Ca in Eq. (1.14) with an 
isotropic term kCA and also for cTc = 0 and = 1. These restrictions should be lifted in 
future for a more realistic calculation. However, this calculation provides support for the 
Kolmogorov-like phenomenology. 

Kraichnan’s Direct interaction approximation (refer to Leslie [35] for details) was applied 
to MHD by Verma and Bhattacharjee [50]. Assuming E{k) ~ k~^l^ they computed the 
Kolmogorov constants C* [see Eq. (1.42)] and found that C* are not universal but depend 
on the Alfven ratio. In our simulations we compute (7* for different ac and show that they 
have a slight dependence on ac- Assuming a k~^/^ spectra, these constants had also been 
calculated earlier [11, 71] in a decaying simulation of 2-D MHD turbulence; however, the ; 
errors involved were large. 

The EDQNM scheme was applied to MHD turbulence by Grappin et al. [4], and the 
spectral exponents m* were obtained by numerically solving the EDQNM equations. They ! 
found that m'^ ~ m~ when (Tc ~ 0, but as cr^ — > 1, the exponent m'^ — >■ 3, and m~ O, ! 
and for any intermediate cTc the m'*’ -|- m~ = 3 and H’^/H" = EDQNM thus gives 

results consistent with the phenomenology proposed by Grappin. ^ 

i 

To conclude this section, numerical simulations [11-13], Solar Wind data [14-16], and* 
also Renormalization group calculations [26] appear to support Kolmogorov-like phenomen-' 
ology over the KID phenomenology, atleast for low <Tc. At high cTc too, indications are that 
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KID phenomenology is not applicable [11]- However, further support is required !»>■ r-ur,..fe 
measurements of spectral exponents or alternatively by testing the flux i' m - h.p^ :n >‘i lit Ipq 

by the phenomenologies, and also by stronger theoretical arguments. 


1.5 Energy cascades in MHD turbulence 

In fluid turbulence the dynamics of energy transfer has been well studied. Wv had discussed 
in Section 1.3 that in 3-D fluid turbulence the kinetic energy is transferred from large scales 
to small scales, whereas in two dimensions there is an inverse cascade of kinetic t*nergy from 
small scales to large scales [44]. Contrary to fluid turbulence, the direct numerical simulations 
(DNS) of MHD turbulence show that the total energy is transferred from large scales to small 
scales both in 2-D as well as 3-D turbulence [11, 71], and that there is a constant cascafle rate 
of energy in the inertial range [11]. There have been theoretical predictions of the magnitude 
and directions of only a few of the various fluxes [17, 18, 72]. 

In Section 1.2.1 we had stated the conserved quantities for inviscid 2-D and 3-1) MHD 
equations. If turbulence is being forced at the wavenumber k/, then it is observed that in 
2-D MHD turbulence, the mean-square vector potential A [Eq. (1.23)] has an inverse cascade 
to wavenumbers less than ky, and in 3-D MHD turbulence the magnetic helicity [Eq. (1.22] 
has an inverse cascade. These cascades are predicted by Absolute Equilibrium theorie.s [73], 
and closure calculations [17, 18], and have also been observed in numerical simulations [74, 
75]. One consequence of the inverse cascade of A is that the magnetic energy spectrum 
E (k) ^ k / for A; < ky [75, 64]. This spectrum has been observed in numerical simulations 
[75]. In our simulations we do not force the magnetic energy (but force the kinetic energy), 
and hence the inverse cascade regime is not found. 

The velocity field evolves due to two nonlinear terms (u.V)u and (b.V)b in Eq. (1.14), 
and the magnetic field evolves due to the terms (u.V)b and (b.V)u in Eq. (1.15) [see Sec- 
tion 1.2]. The term (u.V)u transfers kinetic energy between different scales of velocity fields; 
(b.V)b and (b.V)u transfer energy from the velocity field to the magnetic field and vice 
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versa; (u. V)b transfers magnetic energy between different scales of magnetic fields. Recently 
Ishizawa and Hattori [20] studied some of the kinetic and magnetic energy transfers arising 
from the non-linear terms. In this thesis (Chapter 4) we have carried out a detailed in- 
vestigation of the various energy cascades and shell-to-shell energy transfers in 2-D MHD 
turbulence. 

Pouquet et al. [17] studied the energy transfer between large and small scales using 3-D 
EDQNM closure calculations. They found that the large scales of the magnetic field gain 
energy in the presence of small-scale residual helicity (which is the difference between kinetic 
helicity and magnetic helicity)^. The large scale magnetic field, in turn, was found to enhance 
the energy exchange between the small-scale magnetic and the velocity fields. This energy 
exchange resulted in equipartition of small-scale kinetic and magnetic energies. By drawing 
an analogy between the MHD equation for the magnetic field [see Eq. (1.2)] and the vorticity 
equation"*, Batchelor [72] argued that the transfer between kinetic and magnetic energies takes 
place primarily at small scales. A similar conjecture was made by Pouquet and Patterson 
[21]. They also conjectured that an inverse cascade of energy from small-scale to large- 
scale magnetic field is the mechanism responsible for the ehancement of large-scale magnetic 
energy. Contrary to the prediction from 3-D EDQNM of Pouquet et al. [17], simulations 
of decaying turbulence showed that it is the magnetic helicity and not the residual helicity , 
which is important for the growth of large-scale magnetic field [21]. 

In a 2-D EDQNM study, Pouquet [18] obtained eddy viscosities for MHD turbulence. 
She found that the small-scale magnetic energy acts like a negative eddy viscosity on the, 
large-scale magnetic energy. The inverse cascade of .4. and hence the ehnancement of large- 

I 

scale magnetic energy, was conjectured to arise due to the destablization of the of large-scale 

magnetic field by the small-scale magnetic field. She also found that the small-scale kinetic 

energy has the effect of a positive eddy viscosity on the large-scale magnetic energy. Ishizawa! 

^The kinetic helicity is defined as / v.wdx and magnetic helicity is given by Eq. (1-22) 

^The equation for the vorticity evolution can be obtmned by taking the curl of the Navier-Stokes equation 
(see Lumley [76]) 
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and Hattori [19] resolved the calculations further and obtained eddy \-iscosity diu' 1c> eacli of 
the nonlinear terms in Eqs. (1.1)-(1.2) of Chapter 1. They found that the edrly visrositN' due 
to (b.V)u is positive leading to a transfer of energy from large-scale magnetic field to small- 
scale magnetic field. In their calculation, energy was also found to be transferrt'd from the 
small-scale velocity field to the large-scale magnetic field. Both, Bouquet's [L^; and Ishizawa 
and Hattori's [19] calculations give a non-local energy transfer from small-scale \elority field 
to the large-scale \-elocity field. 

In a recent work Ishizawa and Hattori [20] employed the wavelet basis to investigate 
energy transfer m 2-D MHD turbulence. There are some similarities between the observations 
made by us in this thesis and those in the work of Ishizawa and Hattori^ Using our formalism 
we do a thorough investigation of all the energy transfers (Chapter 4). The wavelet basis 
used by Ishizawa and Hattori [20] is a useful tool for investigating spectral propcrl ie.s and 
energy transfer properties in different regions of the flow. In Ishizawa’s work, tht* flow was 
divided into a turbulent and a coherent region which are distinguished by the level of vorticity 
fluctuations in that region — turbulence is known to have high levels of vorticity fluctuations 
[76]. They found that the energy transfer between velocity/magnetic scales occurs more 
efficient!}' in the turbulent region and less efficiently in the coherent region. 

In some of the above mentioned studies [17, 21] a distinction was not made between 
the energies trctnsferred to a mode from the kinetic and the magnetic fields. Also, the energy 
transferred into a mode k from different wave niunber regions was not separately considered — 
only the net energy transfer into a wave number k was computed [20, 21]. The EDQNM 
closure calculations dealt mainly with nonlocal energy transfers, i.e., between large scales 
and small scales [17-19]. In another study, Frick and Sokoloff [77] solved a shell model of 
MHD turbulence and calculated only those fluxes transferring the kinetic energy between the 
velocity modes, and the magnetic energy between the magnetic modes. 

In our simulations we investigate the (1) various energy fluxes arising within and between 
^However, we became aware of Ishizawa’s results only after the completion of our work 
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the velocity and magnetic fields, and (2) the fine-grained (considering many wave number 
shells) energy transfer between magnetic and kinetic shells. This will give a more informed 
picture of the physics of energy transfer in MHD turbulence. 

In 3-D MHD turbulent flows, beyond a critical magnetic Reynolds number, the magnetic 
field can grow and reach a steady state [21, 78-80] — this is called the dynamo effect and 
is believed to be the mechanism responsible for the generation of magnetic fields within 
astrophysical objects like the sun, earth and galaxies. The non-linear energy transfer in 
MHD turbulence is important for understanding the growth process of magnetic energt' in a 
dynamo. The magnetic helicity is believed to play a crucial role in the generation of large-scale 
magnetic field [81, 24]. In two-dimensional flows there is no magnetic helicity; the magnetic 
energy cannot be sustained in two-dimensional MHD turbulence. Although, it is not possible 
to indefinitely maintain a steady state magnetic field in two-dimensional turbulence [82'. a 
steady state can be achieved for a finite period of time [18]. This, coupled with the similarity 
in physics of the E'^ and E~ energy cascades in 2-D and 3-D, allows us to probe the non-linear 
energy transfer in real MHD flows through 2-D simulations. In Chapter 4, we have calculated 
the energy flux into a magnetic-mode sphere. We have investigated the contribution to this 
flux from the fluid-mode sphere, from the modes outside this sphere, and from the modes 
outside the magnetic-mode sphere. Our results shed light into the amplification of magnetic 
energy. The methods used here are completely generalizable to the 3-D ca.se and can be used 
to directly study the dynamo problem at a later date. 

1.6 Topics beyond the scope of this thesis 

In this thesis we investigate the energy spectra and the energy cascades for the simplest 
example of incompressible MHD turbulence. We do not include a mean-magnetic field, 
shear, current sheets, etc. in our simulations. 

Most astrophysical plasmas (for example, the Solar Wind) evolve in presence of a mean 
magnetic field. Turbulence with a finite Bo is anisotropic due to a preferred direction in 
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space making it simultaneously harder to study and richer in its physics. The niean rna-metic 
field affects both the energy cascade and the energy spectrum. I he Solar wiim .c,-..; ...noa^ 
show that the intensity of fluctuations transverse to the mean magnetic field iarver than 
the intensity parallel to the mean magnetic field [83]. Thus, the energ.v sp.-rtra for the 
fluctuations transverse and parallel to Bo will not be the same. The energy spectra vvd! also 
have an angular dependence. The energy cascades have been inve.stigated tc> sunie deuree and 
interesting observations have been made. Numerical simulations and perturl>ati\ e c„iniliitions 
show that energy cascade occurs more efficiently between the modc.s tran.sv(*r .''0 to Bji and 


less efficiently between modes parallel to Bq [84, 85]. For more details on the iiiv<>viigatioii 
of turbulence with Bo, the reader is referred to Shebalin ei ai [84], Oughton if al. 
Montgomery and Turner [86], Oughton et al. [87], Mattaeus ei al. [88]. and Kinn<‘\’ and 
Mcwilliams [89]. 


1.7 Outline of the thesis 

In Chapter 2 we will compare the predictions of various turbulence phenoitM‘uu!ogir«- with 
the numerical results of our forced, 2-D MHD turbulence simulations, .'\fter .shinviug that 
the Kolmogorov phenomenology is preferred over KID phenomenology, we will .study tiie 
dependence of the Kolmogorov’s constants on the normalised cross-helicity values. 

In Chapter 3 we will devise an approach for describing the energy transfer between a 
pair of modes of the triad where we will introduce the idea of an ‘effective’ tran.sfer. We will 
show how the effective transfer can be used to re-define shell- to-shell energy transfer [39]. and 
also define energy cascade rates. 

In Chapter 4 we will numerically compute the shell-to-shell energy transfer rates and 
energy cascade rates defined in Chapter 3 for forced 2-D MHD turbulence. 

In Chapter 5 we will study the evolution of global quantities for different initial condi- 
tions. We will show that the evolution of global quantities can sometimes depend on subtle 
features of the initial state. 
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In Chapter 6 we will conclude with a summary of our findings and state some of the 
possible future investigations. 




Chapter 2 

A numerical validation of MHD 
turbulence phenomenologies. 

2 . 1 Intro duction 

In Chapter 1 we described the existing MHD turbulence phenomenologies. Kriiichnaii, Irosli- 
nikov, and Dobrowolny et al. (KID) [1-3] predicted that the energy spcctnun rxp. is 
-3/2 (Eq. (1.40) of Chapter 1). Their phenomenology predicts = H" independent of 
E+/E' ratio. In contrast to this phenomenology, Marsch, Matthaeus and Zhou. Zljon atid 
Matthaeus [5-7] argued that the energy spectra E^{k) ~ (Eq. (1.42) of chapt<*r 1 J. .An 
extension of the KID phenomenology due to Grappin et al. [90] predicted that the fitixes 0=^ 
and the exponents in E^{k) ~ satisfy the relationship 

As discussed in the introduction, the solar wind observations indicate that the spt'ctral 
indices are closer to —5/3 than to —3/2 [14, 15]. Various attempts using numerical simu- 
lations have been made to test which one of the phenomenologies is correct. In simulations 
of Pouquet et al. [9], Politano ef al. [10], and Verma et al. [11], reliable estimate.s of 
the spectral exponents could not be obtained. Biskamp and Welter [8] concluded that the 
-3/2 scaling is valid. However, Verma et al. [11] studied the energy fluxes and found that 
n+ differs significantly from H', in contradiction to the predictions of KID phenomenology 
[Eq. (1.39) of Chapter 1]; For small <jc they also found that E+ /£- ~ (n+ /E~f, as predicted 
by Kolmogorov like-phenomenology [from Eq. (1.42) of Chapter 1]. From these observations 



2.2 Energy spectra and energy flux expressions 


29 


Verma et al. concluded that Kolmogorov’s phenomenology is preferred over KID’s phenomen- 
ology. Recently large Reynolds number simulations were done by Politano et al. [12] for 2-D 
MHD turbulence, and by Muller and Biskamp [13] for' 3-D MHD turbulence. They found 
that the spectral indices are closer to —5/3 than to —3/2. However, these results are limited 
to (Tc o; 0 fluctuations. 

Analytically Grappin et al. [90] applied EDQNM closure scheme to MHD equations 
and found that n''‘/n“ = rn^ Jm~ . These predictions have not been observed in the solar 
wind. Recently. RG calculations of Verma [26] shows that Eq. (1.42) of Kolmogorov-like 
phenomenology is a consistent solution of RG equations, but KID's formula [(1.40)] is nor. 

In this chapter we will investigate the energy flu.xes (also termed as energy cascade rate) 
and energy spectra to test which of the phenomenologies is closer to numerical results for 
large cases. From the spectra and fluxes we will also calculate some of the constants. Our 
method is the same as that of Verma et al. [11]. However, we time-average the energy spectra 
and the fluxes over a steady state which improves the quality of data and also test Grappin 
et al.'s phenomenology. 

The plan of this chapter is as follows: In Section 2.2 we define energy spectra and energ}' 
fluxes. In Section 2.3.1 we describe the simulation method used in this thesis. In the present 
chapter and in Chapter 4 we study forced turbulence. The description of the forcing employed 
in the simulations is given in Section 2.3.2. Section 2.4 contains the results. The discussion 
follows in Section 2.5. 

2.2 Energy spectra and energy flux expressions 

In this section we will define the quantities which we compute and study in this chapter.; 
Energy flux is defined as the energy lost by a sphere Fourier space. In this section we will 
give the expressions for the fluxes corresponding to z* energies. We had described in Chapter 
1 that the energy densities of z*, denoted by E'^, are defined as 
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Z JumU'ot 

can be expressed in terms of the Fourier components z |k ) as 

B* = \J |z*(I<)Pdk, 

where we define the quantity E^{k) = |z=^(k)|V2 as the modal eno:cy s|*r.-:r.i. a::.) 


E^{k) = lf |z*(k)!‘dk 

2 J|k|=:Jt 

IS the omni-directional energy spectra. In this chapter our interest is in ’ ni’.i.i. 

We have performed our simulations in two-dimensions with periodic hir.;:;'':.'; tondition.s. 
W’e choose our box-size to be 27r. For such flows the wavenumber spare i.** v r.lr 
and ky = {2ir/L)ny, where n^, Uy are integers. For discrete wave mmiber spa<'e. the omni- 
directional energy spectra in Eq. (2.3) can be written as 

k+l/2 

E £±(k). (2.JI 

k^l/2 

The energy flux of can be obtained from the MHD equation using its tlefiniiion. 
The equation for the evolution of E±(k) can be derived from the MHD equal ions ( 1 . of 
Chapter 1. 


dE^{k) 


1 k+p-i-q=o 

+2jAk^E^{k) + 2un^E^ik)=^-t Y: /?e(i[k.z^(q)][z±(k).z^(p)]).f2.5) 


Integrating the above equation from k - 0 to k = K. the right hand side of the equation 
gives energy flux out of the sphere K to the modes outside. Hence, the fluxes 11^ are equal 
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to 

n*(A') = (i[k.z^(q)][z*(k).z*(p)]) . (2.6) 

k P 

where k + p + q = 0. We will compute the above energy fluxes in our numerical simulations. 

In this chapter we will discuss the properties of the energy spectra E^{k) and the energy 
fluxes n^(A'). Then we will calculate the Kolmogorov’s constants using these quantities. 
Our analysis will be carried out for a range of normalised cross-helicity cTc values (defined in 
Section 1.2.1 of Chapter 1). 

2.3 Simulation details 

2.3.1 Numerical method 

In our simulations we use the MHD equations given in terms z"^ and z~. i.e., 

^q:(Bo.V)z±+(z^.v)z’t = -Vp + i/±V2z± + z/:pW 

+ V^z± + V^z^ + F±, (.2.7) 

V.z^ = 0. (2.8) 

The terms are the forcing functions. The corresponding forcing functions for the u and the 
b fields are related to F^, i.e., F“ = (F‘'' + F“)/2, and F*" = (F'*‘ — F“)/2 respectively. In our 
simulations we force only the velocity and not the magnetic field (i.e., F*^ = 0). Consequently, 
we get F"^ = F“ = F“ = F. 

In the numerical simulations we time evolve Eq. (2.7) in spectral space. Taking the 
Fourier transform of the above equations, we write 

^ = Mf-D± + Ff. (2,9) 


k.zj = 0, 


( 2 . 10 ) 
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where 

IS the Fourier transform of the nonlinear terms, and 




V 







' 2 . 12 ] 


IS the Fourier transform of the dissipative terms. The variables zf . {<z'* .T” 7, ' jn, are Jhf 
Fourier transforms, pk is determined using the incompressibility const r.iin?. 

The Eq. (2.9) is time evolved using the Adam-Bashforth .scheme fur tiie tiun’ine.ii terms 
and Crank-Nicholson for the dissipative terms. Both the.se scheme.'' haw an acc'.n'Hcy of 
OiAt^). Using the two schemes we obtain at time step n + 1 from z* at imie ''tep u and 
n — 1 using 


±(n+l) ±{n) 

It 




At 




n+l) 


0±i"i 


)-F^ 


i, i I 


.13) 


The nonlinear terms are usually computed using the p.seudo-spect ra! method (iU) which 
we illustrate below. For the 2-D case, u and b do not have any component along the 
axis, and there is also no variation along z-axis, i.e., u = {u^(x,y, f). t/. n. 0). and 

b = (6^(x.i/.t).6j,(a-,y,t),0). The x and y-components of Elsasser variable.s in a 2-1) MUD 


equation are 


dt 




) + — 




+l'±lkl,) ve-j + {u^/k't,) w.-j + F* 

(2.14) 
dv 

+ l'±K) V.-J + + cy 

(2.15) 
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In Fourier space the nonlinear terms involve a convolution sum. For example, the term 
expressed in the Fourier space as the convolution sum (k — p) 

(also see the Fourier space MHD equation written in the introduction). To compute this term 
for each k by direct summation requires 0{N'^) operations, where N is the number of grid 
points in Fourier space. The number of operations is reduced to 0{NlogN) by using the 
pseudo-spectral method. The are transformed from Fourier space to real space using 
inverse Fast Fourier transforms. The products etc. are formed in real space. The 

resulting product is transformed back to the Fourier space using a forward fast Fourier 
transform and then the Fourier transformed derivative d/dx is obtained by multiplying the 
transformed product by ik^. The number of operations required is thus reduced to 0{NlogX ) 
h In the above procedure the aliasing error is known to modify the convolution sum by adding 
the contribution of higher wavenumbers to the lower wavenumbers [91]. In order to remo\ e 
the aliasing errors arising in the pseudo-spectral method we use square truncation wherein 
all the modes with > N/3 or > N/3 are set equal to zero [91]. 

In direct numerical simulations hyperviscosity vk'^/kl^ and hyper-resistivity gk’^/kl^ is 
added to the equations. It damps the higher wavenumbers while not modifying the lower 
wavenumbers. In this way, it becomes possible to simulate higher Reynolds number turbu- 
lence. In general higher orders of hyperviscosity / k'^'^ are often used. However, much 

higher orders of hyperviscosity need to be avoided as it has been pointed out by Biskamp tt 
al. [92] that the effect of hyperviscosity on the inertial range is large for high \^lues of p. We 
have used the ralue of p = 1 in all our simulations. 

The simulations are carried out in a periodic box of size 27r x 27r. .\11 the quantities are ! 
non-dimensionalised using the initial total energy and length scale of 27r. 

Hhe vorticity- vector potential form of the 2-D MHD equations which is seemingly simpler, is infact more' 
expensive computationally. A larger number of FFT’s are required to compute the nonlinear terms in the 
vorticity-vector potential formulation 
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2.3.2 Steady state 

The normalised cross helicity (T, is defined as ctc = 2HJE. where H, = /:’ - /:’■ i> tie- ru^s 
hehcity and E is the total energy. In decaying turbulence the tt)tal tuiergy due lo 

the viscous and resistive dissipation, whereas cr^ often shows the tendency to increase- with 
time. We shall discuss the issue of increase of cross-helicity in greater d<‘tai! in (.’hapter o. 
In simulations it is observed that the decay rate of the majority species E* ( if rr • . p j 
larger than the decay rate of the minority species E~ . Therefore. //.- ilerreases with time. 
To obtain a steady state we inject kinetic energy and cross-helicity (//i-j at a rate such that 
both E and Ue are maintained at a statistically steady value. 

.At each time step we construct F as an uncorrelated random function that is diwigc-nce 
free fi.e.. V.F = 0). The T-component of F is determined at every time ste}> by 

F,= 12.161 

where is equal to the average energy input rate per mode, and phase q is a uniformly 
distributed random variable between 0 and 27r. The y-component of the forcing function is 
obtained by using the condition of zero divergence, i.e., 

= - ( j) (2.17) 

We force all the modes in the annulus 4 < A; < 5. The value of is determined from the 
a\-erage rate of the total energy input. 

The simulation method described above is common to all our simulations described in 
this thesis. The specific values of r/’s, //’s, keq, and the forcing parameters for each simulation 
will be stated in the next section. 

2.4 Results 

In this section we will study the energy spectra £±(fc) and the energy fluxes n± over a 
range of a, \alues. Using these results we will carry out a comparative test of KID and the 
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Kolmogorov phenomenologies and compute the appropriate constants. 

We have obtained a steady state for the approximate cTc ^'alues listed in Table 2.1. The 
.A.lfven ratio in he steady state are also shown in the table. Note that in our simulations higher 
(7e states also have a higher r^. We have performed all the runs at a tmit Prandtl number 
{u = ix). The values of v and the hyper- viscosity parameter k^q are also shown in the table. 
The fluxes and energy spectra are averaged in the steady state over a time inter\‘al tavg- 


Runs 

o-c 

1/ 

h 

‘^eq 

ta 

dt 

^avg 

R1 

0.1 

5 X 10-® 

16 

0.2 

5 X 10-* 

6 

R4 

0.4 

5 X 10"® 

16 

0.5 

2.5 X 10-“ 

6 

R7 

0.7 

5 X 10-® 

16 

0.65 

2.5 X 10-“ 

6 

R9 

0.9 

5 X 10"® 

14 

0.8 

2.5 X 10-“ 

6 

R95 

0.95 

5 X 10-® 

14 

0.9 

2.5 X 10-“ 

3 


Table 2.1: The parameters for our simulations and the steady state values of cTc and 

The compensated energy spectra , E^{k)k°’ with a = 5/3 and 3/2 have been plotted in 
Figs. 2. 1-2.5 for the runs R1-R95. From these plots of E^k‘^ we are unable to distinguish 
between the —3/2 and —5/3 exponents, as the inertial range is not wide enough. Even though 
the exponents m* cannot be accurately determined from these plots, it is clear that m"*" and 
m~ are close to each other, i.e., m'^ ~ m~. For low cTc fluctuations, the spectral exponents 
have been measured by Politano et al. [12] to be close to —5/3, with m’*" = —1.69 and 
m~ = —1.7, as predicted by Kolmogorov-like phenomenology. In the same simulations of 
Politano et al. [12], the structure function of fourth order, i.e.. |Az*(r)|‘‘ were found to be 
proportional to r, in accordance with the predictions of KID phenomenology. However, we 
have argued in Chapter 1 (see Section 1.4.2) that the predictions of KID phenomenology ' 
regarding the fluxes are not in agreement with simulations, and hence the agreement of 




k 


Figure 2.1: E^{k)k°‘ versus k for run R1 with ac = 0.1. The solid [E"^) and c 
dotted line (£’“) corresponds to q = 5/3, and dashed (£''*') and widely sjiacer 
{E~ I correspond to a = 3/2. 
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Figure 2.8: The fluxes 11^ for run R95 with steady state <7c — 0.95. 


in KID phenomenology and simulations observed by Politano et al. [12] may be fortuitous. 
The small deviation of m* from —5/3 observed in Politano’s simulation can arise because of 
intermittency. 

For high cTc fluctuations attempts have been made to measure from the spectral 
plots [9-11], but no conclusive statements were made from these measurements. Instead of 
measuring the spectral exponents directly, Verma et al. [11] had tested the predictions of 
KID and Kolmogorov phenomenologies regarding the fluxes H''" and 11“. We will follow their 
approach. 

We have plotted the fluxes for runs R1-R95 in Figs. 2. 6-2.8. We find that 11“ 7^ 11“ for 
all the (Tc’s listed in Table 2.2. They become increasingly unequal with an increase in <7c. -As 
an example, for cTc = 0.4, the ratio n‘'‘/n“ = 1.08, and for ac = 0.95, n'^/n“ = 4. These 
observations had earlier been made by Verma et al. [11] in their simulations. Since TV' 
they had concluded from these results that KID phenomenology is not applicable to MHD 
turbulence. However, an extension of KID phenomenology due to Grappin et al. [4] predicts 
inequality of fluxes. This phenomenology was not tested by Verma et al.. According to 
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Grappin's phenomenology [4] m/U' = m+/m-. However, we u^.^ervr Huit uti: 
do not support this prediction, since m'^ ~ 7V~ but ^ 11 


Tims, tiif' f!u\ >■ 


•tiOI!,*. 




obser\ed in the simulations cannot be explained succes.sfully f>y hii\' uf ti;-' ex;*.;!;!.' m.-deis 
based on KID phenomenology. 

Inequality of fluxes 11+ and II' is consistent with the prediction.^ of Ki'd- V piie- 

nomenology [5-7]. For low values of Uc, the constants C+ and C (in I'.q. * i . tJ‘ uf t !!;i|t!er 

1) should be expected to be nearly equal. It follows from Eq. (!.42) of ('hapten I that the 
relationship £’+/£'" = (n+/n”)^ will hold. From Table 2.2. where we* hine* tabuiaU'ii these 
ratios, we find that our simulations confirm this prediction of the Kolniogoren like :;u:;a‘L- 

ology. For highest crc(= 0.95) in our simulations (see Table 2.2) we find that and 

(n+/n~)' differ approximately by a factor of two. The predictions of Koiniouico*. -like* plu*- 
nomenology will be consistent with our simulation results if C+ and ('~ are unenjual for hitrii 
cTc cases. Earlier Verma and Bhattacharjee [50] had applied DLA. to calcuiate* the KwlnMuuttA 
constant for MHD turbulence. They show in their calculations that the C'* depend on .'Mfeen 
ratio and argue that it will also depend on < 7 ^. 

On the basis of the results of reference [12, 11] and our own results, there i.s justification 
for supporting the view that E^(k) follow Kolmogorov’s -5/3 scaling at a!! value.'; of rr,.. 

In the Kolmogorov-like phenomenology for MHD turbulence the energy .sjjertruin for 
is given as 

E^(k) = (2.18) 


where C are the Kolmogorov's constants. From the above equation we can write as 

E=^(k)k^/^ 

(n+)V3(n=F)-2/3- (2-19) 

In fluid turbulence the Kolmogorov constant are found to be universal. We had pointed 
out abot-e that in MHD turbulence the constants may depend on cr^. Now we will study 
the dependence of G+ on (7^. Verma et al. [11] had previously obtained from deca}'ing 
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simulations. However, because of the large fluctuations in the energy spectrum and the fluxes 
there were large errors in the values of these constants. In our simulations and 

will be averaged over the steady state and (7* then calculated from Eq. (2.19). 


^^1 




IBB 

H- 

IBI 


C-/C+ 


(n+/n- 

R1 

0.1 

0.48 

0.41 

0.041 

0.038 

4.0 

4.1 


1.2 

1.2 

R4 

0.4 



0.06 

0.042 

3.3 

4.0 

1.2 

1.7 


R7 

0.7 



0.07 

0.03 

3.0 

4.4 

1.4 

4.1 

5.4 

R9 

0.9 

1.18 


0.06 

0.017 

3.3 

4.2 

1.3 

5.1 

12.4 

R95 

0.95 

1.25 


0.048 

0.01 

3.3 

5.7 

1.7 

13.9 



Table 2.2: The E^{k)k^^^ and H^ obtained from the simulations, and the Kolmogorov's 
constants C* calculated from these values. The errors in due to the fluctuations are 
indicated in text. 

The constants were obtained using Eq. (2.19) by first performing a fit of a constant 
function, i.e., f{x) = constant to E^k^^^ vs. k and to n*(/j) vs. k. The fit was performed 
in the inertial range over a wavenumber interval typically equal to 20. The best estimates of 
C'* thus computed are given in Table 2.2. The worst error estimate of C* were computed 
by considering the maximum error in E^k^!^ and n*(A:). The error in C* calculated in 
this manner was roughly in the range of 0.4 to 0.8 for the various «Tc cases. The following 
observations can be made from Table 2.2. 

The constant C'*' and C~ are unequal at high values of cr., i.e., ^ C~. and C~ is 

larger than C'^ — these results are consistent with those from Verma et al.’s [11] decaying 
simulations. At ~ 0 (R1 in Table 2.2), C"*" ~ C” = C as expected with C ~ 4. The 
constant C~ does not show a tendency to vary upto (Tc = 0.9. However, beyond = 0.9. 
C~ has a dependence on cr^. It increases from a value of 4.2 ± 0.5 at ctc = 0.9 to a value 
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0 . 7 ± 0.8 at cTc = 0 . 95 . On the contrary varies at low <7^ value,*. i! ;ai !u .5 - y .5 
(7c — 0.1 and decreases to 3.3 ± 0.3 at ctc = 0.4. At higher c,. values. ( “ liu*-**. lu.’ anv 
dependence on (7c- and is roughly 3.3 ± 0.8. Therefore. C~ is iuilepe: • »>: r },)» ;,av p 

and C'^ is independent of Uc at high <Tc’s. Verma d al. fllj had earliej foiiiid th.i! lu: ]uw rr. 
case C'*' ~ C“ ~ 5 — 7, and for high ac case C'^ 4 — 5 and C" ti: 8 - Id. 1 hen Iesu]J^ 

differ from our present results; this may be due to large error bars in the:: 'hu, . 

The ratio E'^ / E~ can be obtained from the definition cr- = i E'^ - F~ < ( /• " - /■, • i. 7 |jp 
\alue E*IE~ corresponding to the cTc’s of Run R1-R95 are listed in Table H.'J, The ratio 
C~IC'^ IS also shown in the table. In Fig. 2.9 we have plotted {'’- /(’•*■ / " /;■ foj. 

(7c = 0. to 0.7. From the figure there appears to be a linear relation.*^}!!}) between ' ' f and 
E^ IE~ at low- values of cr^. The data points for higher values of cr. do not appear tt; follow 
this linear dependence between C“/C+ and E^^fE' (not shown in the figure |. 1 he best fit 
between the two parameters is found to be C"/C+ = (0.09 ±0.01 )(£'■*•/£■" ; -♦ fft.'ld - y y:}). 



the computed*^uls^fo(*7 To f to 0 7 afd f diamonds are 

points [C-IC^ = (0.09 ± 0.01)(£7/n“ (^93 ±7 o3^ T “ ‘”7 

from (Tc = 0. to 0.95 ^ ^ * O.uSjj. The figure inset shows the data 
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2.5 Discussion 

In this chapter we studied the numerical results of energy spectra E^{k) and fluxes 11* 
in 2-D MHD turbulence with kinetic forcing. From the energ}' spectra, it was difficult to 
distinguish between -3/2 and -5/3 exponents as predicted by the KID phenomenology and the 
Kolmogorov-like phenomenology, respectively. The spectral exponents m'^ and m~ however 
were found to be close to each other. The predictions regarding fluxes were earlier tested by 
Verma et al. in a decaying simulation, and KID phenomenology was found to be inconsistent 
with the numerical results. We performed the same test in a forced simulation and found KID 
phenomenology to be inconsistent with the numerical results, as expected from earlier results 
of Verma et al. [11]. Note, however, that our simulations were done in the steady state. We 
have time-averaged the energy spectra and the flu.xes resulting in less errors compared to 
Verma et al. [11]. We found that Grappin et al. ’s extension of KID phenomenology was also 
not consistent with the numerical results. We concluded that KID phenomenology [1-3, 90] 
cannot explain the numerical results. However, as observed by Verma et al. [11], we also 
find that E'^/E~ ~ (H'^/n")^, upto a factor of atmost two. Therefore, we tend to believe 
that Kolmogorov-like phenomenology is more suitable for MHD turbulence. Our belief is 
reinforced by recent high Reynolds number simulations by Politano et al. [12], and Muller and 
Biskamp [13]. Unfortunately, these high Reynolds number simulations have been done only for 
small (Jc- .A. similar analysis should be carried out for large (Jc- Recently. Verma [26] has shown 
that the spectra predicted by Kolmogorov-like phenomenology is a consistent solution 
of RG equations of MHD. This approach is based on renormalisation of mean magnetic field. 
Unfortunately, this approach also assumes ctc = 0, = 1. Generalization of the above RG 

argument is needed for (Tc ^ 0 cases. 

The Kolmogorov’s —5/3 scaling had recently been observed in high Reynolds number 
simulations [12, 13] for low ctc cases, and a theoretical reasoning had also been provided in 
RG calculations in support of —5/3 scaling [26]. For high cr^ cases our simulation results 
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preferred Kolmogorov-like phenomenolog>' over KID type pheno^u•Ilok)pH•^. l-rum ^1] the 
earlier evidences and our simulation results we concluded that >ikf p.i' 

logy probably provides the correct description for MHD turbulence. 

We computed the Kolmogorov’s constants for various cases from 0.1 ti- O.iir,. We 
found that C"*" varies at low values of but stays constant beyond (t- 0.4, v. hih* ( stays 


constant upto (Tc — 0.9 and varies beyond that. The value of C of C in tledi ■ 
legimes were approximately equal to 3.3 i 0.8 and 4.2 i 0.5 respecli\el) < the v.tusf error 
estimates are quoted). These constants were earlier evaluated by Xernia ft al ! i in a 
decaving simulation. However, in their simulations the fluctuations in 11* ami /■. * ■ /r 
were too large to provide the trend of the constants with variation in cr^. In our 'Iri.uhttkuis. 
\ve have averaged the energy spectra and fluxes over the steady slate maintained h\ forcing 
and obtained a clearer picture of the variation of the constants with <7^. 

The Kolmogorov's constants are also expected to depend upon the Alfven ratio j.'H)]. It 
Ks desirable to obtain both ac and dependence of C*. However, we have iu>t studied this 
dependence in our numerical simulations. In all our numerical simulation.s =■■ 0.2 “ (i.9 and 
we could not vary it independently of ac- The reason for this is the difficulty in {>l>tainiiig the 
desired steady state Alfven ratio in simulations. It is not clear to us at the moment whether 
It i'. possible to maintain a desired level of Alfven ratio by external forcing. 

In this chapter we studied some properties of the energy fluxes (cascade rates) of z* 
eucMgies and using these properties we attempted to see which of the MHD turbulence phe- 
iiomenologies is correct. In our analysis we ignored various energy fluxes and energy transfer 


lalrt related to u and b. These fluxes are very important for magnetic field amplifif ation 
ob^el\ed in astrophysical objects. In the next chapter we will develop a formalism for studv- 
mp; energy transfer properties between the various scales of u and b fields. In the following 
chapter we will use that formalism for studying energy fluxes (cascade rates) and shell-to-shell 
energy transfer rates in numerical simulations of 2-D MHD turbulence. 



Chapter 3 


A new approach to study energy 
transfer in turbulence 

3 . 1 Introduction 

In fluid and MHD turbulence eddies of various sizes interact anaongst themselves; energy 
is exchanged among them in this process. These interactions arise due to the nonlinear- 
ity present in these systems. The interactions involve wavenumber triads (k,p,q) with 
k -f p -|- q = 0. The complete nonlinear energy transfer term consist of all possible com- 
binations of wavenumber triads that satisfy this triangle equality. 

The combined energy transfer computation to a mode from the other two modes of 
a triad has generally been considered to be fundamental. This formalism has played an 
important role in furthering our understanding of locality in fluid turbulence [39-42], and 
also in analysing subgrid scale eddy viscosity [93, 94]. In this formalism, Batchelor [72] first 
suggested a formula for the energy transfer in fluid turbulence between two regions in Fourier 
space, and Domaradzki and Rogallo [39] calculated shell-to-shell energy transfer rates. Taking 
energy energy transfer between shells as an example, we will show in Section 3.3, that this 
formalism has certain limitations. 

In this chapter we present a scheme to calculate the energy transfer between two modes 
in a triad interaction. Using this scheme we can calculate the energy transfer between two 
shells, cascade rates, and also study locality of energy transfer. 
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One of the main objectives of this thesis is to study cascade rato aiic Ma ’,; ttvshdl 
transfer rates in MHD turbulence. The formulae for the combined eiu-rp>' transfer Uii,. were 
used to study energy transfer in MHD turbulence [11. 17-20]. Ibnuner. ral (a>rade 
rates cannot be calculated using this formalism. Our formalism of "isu-de-t^. tnudr" tran.sfw 
enables us to calculate all these energy transfer rates between the spl.c:.'. and i’n»* ^he!i^. 

This chapter is organized into the following sections, in Section we will .stale the 
known results regarding the combined energy transfer between velocity !su»dt*s itf a triad. 
In Section 3.3 we will present a discussion on the formulae which are User] fur .;>ii 


energ}' transfer between shells in fluid turbulence. In Section 3.4 we will furimuatf* a wa_\ of 
studying energy transfer between a pair of modes within a triad, with the third iiiudr' of the 
triad meditating the transfer. This new method is used to re-define shell-lo-sheli tran‘:f«’r in 
Section 3.5. The known results of energy transfer in a MHD triad are discussed in S<>ctiun 3.5 
and the "mode-to-mode” transfers in MHD in Section 3.7. The re.sults of .Section 3.7 are used 
in Section 3.8 to define shell-to-shell transfer rates and the energy ca.scade rute.s in Mill) 
turbulence. The conclusion to this chapter follows in Section 3.9. 


3.2 Energy transfer in a triad: Known results 


The Xavier Stokes equation in real space are written as 


.^ + (u.V)u = -Vp + i.VV (.111 

where u is the velocity field, and v is the fluid kinematic viscosity. In Fourier space, the 
kinetic energy equations for a Fourier mode is 

+ 2 i5“(k|p,q), 

k+p+q=o ^ 


(3.2) 
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where £'“(k) = |u(k)p/2 is the kinetic energy in Fourier space for mode k. The nonlinear 
terms 5““(k|p, q) is 

^““(kip, q) = -Re (z[k.u(q)][u(k).u(p)] + i[k.u(p)][u(k).u(q)]) (3.3) 

This term represents the combined transfer of kinetic energy from mode p and q to mode k 
[44]. Note that the wavenumber triad k, p, and q should satisfy the condition k + p -f q = 0. 
This nonlinear term satisfies the following detailed conservation properties 

5““(k|p, q) + F““(p|k,q) + 5““(q|k,p) = 0. (3.4) 

The quantity 5““(k|p.q) represents the nonlinear energy transfer from the two modes 
p and q to mode k. It would be useful to know the exact energy transfer between any two 
modes, say from p to k. In Section 3.4 we will explore this idea further. 

3.3 Shell-to-Shell energy transfer in fluid turbulence 

We had mentioned in the last section that the combined energy transfer from modes p and q 
to mode k is given by 5““(k|p, q). Using this result researchers [39] have discussed energy 
transfer between two shells. The quantity 

r;:t = 5EE'y““(k|p.q)- (s-s) 

^ kern pen 

is usually interpreted as the rate of energy transfer from shell n to shell m [39, 72]. Note 
that k-sum is over shell m, p-sum over shell n, and q = — k — p. However, Domaradzki [39] 
points out that it may not be correct to interpret the formula (3.5) as the shell-to-shell energ}' 
transfer. The argument is as follows: 

In energy transfer between two shells m and n, two types of wavenumber triads are 
involved in the energy transfer (see Fig. 3.1 shown below). In triad of type I the wavenumbers 
p, q are located in one shell, and k is located in the other. In type II, wavenumber k is in one 
shell, p or q in the other shell, and the third wavenumber is located outside the two shells. 
In Eq. (3.5) the summation is carried over both kinds of triads. However, it is easy to to see 
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Figure 3.1: two types of triads involved in transfer between shell m and .shell n. Mode q 
could be either inside a shell or could be located outside the shells. 

that the energy transfer from shell n to shell m takes place through both the k-p and k-q 
legs of triad I but only through the k-p leg of triad II. Hence Batchelor’s and Dunjatadzki's 
formalism do not yield correct shell-to-shell energy transfers (as pointed out by Batchelor 
and Domaradzki themselves). 

3.4 “Mode-to-Mode” energy transfer in a triad 

The nonlinear interaction in the Navier-Stokes and the MHD equations are intrinsically three- 
mode interactions. The expressions for the energy transfer to one mode of the triad from 
the other two were presented in Section 3.2, In this section we shall explore the possibility 
of obtaining an expression for the energy* transfer between any two modes within a triad — 
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we will call it the mode-to-mode transfer. We emphasize that this approach is still within 
framework of the triad interaction — that is, the triad is still the fundamental interaction 
of which the mode-to-mode transfer is a part. The energy transfer between two modes of a 
triad by the mediation of the third mode is sought here. 

We shall first consider only the Navier-Stokes equation. In a later section, the discussion 
will be generalised to the MHD equation. 

3,4.1 Definition of Mode-to-Mode transfer in a triad 

Consider a triad k, p, q shown in the figure below. The quantity ^““(k|p|q) in Fig. 3.2 


u(q) 



Figure 3.2: The mode-to-mode energy transfers sought to be determined. 


denotes the energy transferred from mode p to mode k with mode q playing the role of a 
mediator. We wish to obtain an expression for 

"^““’s should satisfy the following relationships : 


1. The sum of ^““(klplq) and "^““(klqlp), which represent energy transfer from mode p 
to mode k and from mode q to mode k respectively, should be equal to the total energy- 
transferred to mode k from modes p and q, i.e., 5'““(k|p.q) [see Eq. (3.3)]. Thus. 

r“(k|p|q) + r“(k|q|p) = 5““(k|p,q), 








(3.6) 


3.4 “Mode-to-Mode” energy transfer in a triad 


52 


^““(pjkjq) + 1S"“(pjq|k) = ^-npik.q). 3,7, 

r”(q|k|p) + r'(qlp|k) = ^--'q^k.p). ,:),S| 


2 . The energ}' transferred from mode p to mode k, i.e., ^“'‘(kipiqi, will ht* equal and 
opposite to the energy transferred from mode k to mode p i.e., 'K“’'>p k qi. i.tKu i-ee 


Fig. 3.3. Thus, 

r“(k|p!q) + r“{p|klq) = 0. .3 9, 

r“(k|q|p) + r“(qlk|p) = 0. 

1^““(p|q|k) + :K““(q|p|k) = 0. . 3 . 11 , 


u(q) 



R““(plklq) 


Figure 3.3: Relationship between the Trajisfer u(k) ^ u(p) is physically the same as 


There are six equations and six unknowns. However, the value of the determinant formed 
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from the Eqs. (3.6)-(3.11) is zero. Therefore we cannot find unique given just these 

equations. However, it is reasonable to expect that there is a definite amount of energy 
transfer from one mode to the other mode, say from mode k to mode p given u(k), u(p), 
u(q). Since Eqs. (3.6)-(3.11) do not yield a unique we need to use constraints based 


on invariance, symmetries, etc. to get a definite using Eqs. (3.6)-(3.11). 

3.4.2 Solutions of equations of mode-to-mode transfer 

One solution of Eqs. (3.6)-(3.11) is 

5 ^““(k|p|q) = -i?e(i[k.u(q)][u(k).u(p)]), (3.12 1 

From the definition of j^““(k[p|q), it directly follows that ^““’s satisfy the following conditions. 

^““(k|p|q) + 5 ^““(k|q|p) = 5““(k|p.q), (3.13) 

r“(p|k|q) + r“(p|q|k) =5““(p|k.q), (3.14) 

r“(q|k|p) + r“(q|p|k) = 5-(q|k.p). (3.15) 

Using the triad relationship k + p + q = 0, and the incompressibility constraint [k.u(k) = Oi. 
it can be seen that ^^““'s, etc. also satisfy the following conditions 

r“(k|p|q)+r“(pik|q) = 0, (3.16) 

r“(k|q|p)+r“(q|k|p) = 0, (3.17) 

r“(p|q|k)+r“(q|p|k) = 0. (3.1SI 


Comparing Eqs. (3.13)-(3.18) with Eqs. (3.6)-(3.11), it is clear that the set of is one 
instance of the :|l““’s, i.e., '|?.““(k|p|q) = ^““(k|p(q). However, this is not a unique solution. 

If another solution "^““(klplq) differs from <?““(k|p|q) by an arbitrary function i-e.. 
)^““(k|p|q) = j?““(k|p|q) + Xa, then by inspection we can easily see that every possible 
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solution of Eqs. (3.6)-(3.11) must be of the form 


^““(k|p|q) = ;”(k|p|q) + -Vi 

(3.19) 

:!t"“(k|q|p) = ^““(kiq|p) - .Vj 

d3.20) 

?t-“(p|klq) = r“(P|k|q) - -Va 

(3.21; 

1l“*(p|q|k) = r“(Piq|kl + .Va 

''3.22' 

?t““{q|k|p) = r“(q|kiPi + .Va 

(3.23) 

r“(q|p|k) = r“(qip|k) - .Va 

i3.2-b 


Note that can depend upon the wavenumber triad k, p, q, and iht* Fourier cijrapcii;. 
ents u(k), u(p), u(q). Xa needs to be determined from other symmetry and invariance 
arguments to uniquely fix 

In Appendix A we have attempted to determine We construct A'j,. by i»b-i>ivin^ 
that it can depend on k, p, q, u(k), u(p), u(q): and it must satisfy rotational im.iri.iiiiT, 
galilean invariance, and it should be finite. In our procedure we write down a general scalar 
function of k, p, q, u(k), u(p), and u(q) with undetermined coefficients. The ctadficients 
can be written as a series in the scalars of the type k.p, k.u(p). u(k).u{p) with unknown 
constants. By imposing invariance, and finiteness constraints we have shown that, upto the 
linear order in these scalars, the coefficients must vanish, giving A';:^ = 0. I’nfortunatcdy, we 
have not been able to determine the higher order coefficients. Hence, from th<‘ <inal\ si.s of 
Appendix A we cannot show for definite that Xa will vanish. However, vve belie\’e that .Aa 
is equal to zero. When Xa = 0, we get a simple relationship, 

r“ = (3.25) 

giving the energy transfer between a pair of modes in a triaxl with the third mode mediating 
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the transfer, as simply But if Xa 0 then the energy transfer picture is more comple.x 
and the energy transfer between a pair of modes will be given by Eqs. (3.19)-(3.24) where 
.Va is unknown. 

We will now give a physical interpretation to the two parts of ^““’s, i.e., s and Aa- 
We see from Fig. 3.4 that <^““(k|p|q) + ATa gets transferred from p to k, ^““(q|k|p) + Aa 


u(q) 



Figure 3.4: Mode-to-mode transfer can be expressed as a sum of circulating transfer and 
effective mode-to-mode transfer. 

gets transferred from k to q, ^““(p|qlk) + Aa gets transferred from q to p. The quantity 
A'a flows along p -> k q — ^ p. circulating around the entire triad without changing the 
energy of any of the modes. Therefore we will call it the Circulating transfer. Of the total 
energy transfer between two modes, 5?““ d- A a, only can bring about a change in modal 
energy. Aa transferred from mode p to mode k is transferred back to mode p via mode q, i.e. 
the mode p transfers Aa directly to mode k, and mode p transfers Aa back to k indirectly 
through mode q. Thus the energy that is effectively transferred from mode p to mode k 
is just r“(k|p|q)- Therefore r"(k|p|q) can be be termed as the ^^effective mode-to-mode 

transfer" from mode p to mode k. 
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In summary, we have attempted to obtain an expre.'^slon for the enrrsiy transfer rate 
between two modes in a triad. To leading order in a series exp.;--:"’ we* ran -<ho 5 v 
that But the ambiguity in remains b<’cause of the luck of a * 

proof. However, the important conclusion is that the mode-to-nio«le <’ii(*rcy tran.-f«T ran he* 
expressed as = 5?“^^ +Xa where 5?““ is the '‘effective mo<le-rt>-nn><!e tra:!>f<*r' and A.^ k 
the ‘circulating transfer’. In the next section we will further devt'liip the nutni:) uf efh'njve. 
transfers by applying it to energy transfer between shells. 

3.5 Shell-to-Shell energy transfer and Cascade rates us- 
ing “mode-to-mode” formalism 

In Section 3.3 we had pointed out the problem in treating the exf)re'-iot; in K<j. a.' 
the shell-to-shell transfer. This expression, commonly taken to be tht* shell- ri/-.hel! transfer, 
actually gives the energy transferred to shell m from both shell n ami the ijuxies q. In thi- 
section we will make use of the idea of “circulating” and “effective mode-to-mode transfer" 
(see Fig. 3.4) to redefine shell-to-shell transfer. 

Consider energy transfer between shells m and n in the figure below. The mode k i.s in 
shell m. the mode p is in shell n, and q could be inside or outside the shells. In terms of the 

mode-to-mode transfer :^““(k|p|q) from mode p to mode k, the energy transfer from shell 
n to shell m can be defined as 

rr, = EEr“(kip|q) 

kcm pen 

where the k-sum is over sheU m, p-sum is over shell o, and k+p+q=0. The quantity r' 
can be written ae a sum of an effective transfer r“(k|p|q) and a circulating transfer A'a. 
We know from the last section that the circulating transfer does not contribute to the energ;' 
change of modes. From Fig. 3.5 we can see that Aa flows from shell m to shell n and then 
ck to m indirectly through the mode q. Therefore the effective energy transfer from 
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Figure 3.5: The circulating transfer and mode-to-mode transfer from shell n to shell m with 
the arrows pointing towards the direction of energy transfer. 

shell m to shell n is just ^““(k|p|q) summed over all k in shell m and all p in shell n. i.e., 

r~ = EEr“(kip|q). (3.2r) 

k£m pen 

where 5 ?““(k|p|q) is the effective mode-to-mode transfer. 

If = 0 then the effective transfer between two shells given by this equation will be the 
same as the same as the total energy transfer between them given by Eq. (3.26). However, if 
-Ya is non-zero then the two are not equal. The effective shell-to-shell transfer [Eq. (3.27)] 
differs from the total transfer [Eq. (3.26)] by the contribution of the circulating transfer. We 
have argued earlier that the circulating transfer does not result in the energy change of any 
of the modes. Hence, from the physical point of view the effective shell-to-shell transfer is a 
very relevant quantity to study. In the next chapter we have studied the effective shell-to-shell 
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transfers in numerical simulations. 


3.5.1 Energy cascade rates 

The kinetic energy cascade rate (H) (or flux) in fluid turbuh'nce is .Iriii.,-’ .»r ratf of joss 
of kinetic energy by a sphere in k-space to the modes out.sid<\ In Hifru!-!:**. 
in fluid turbulence has been computed using 5““(kjp.q) [-‘rx /' • diMux tha! 

by direct integration of Eq. (3.2) that the kinetic energy los! from a -pi.f-if . -f ladio^ h' by 
nonlinear convective transfers is 


n(/v) = - Y1 E 


<3.2?^ I 


rui ><■ n- jl' 


.Although the energy cascade rate in fluid turbulence can be found by the ahioi' furniula. the 
mode-to-mode approach provides a more natural way of looking at the enerKv tranyfers. In 
later sections we will obtain expressions for mode-to-mode tran.sfer in Mill) nubiih'tire and 
define various fluxes which are not accessible to the conventional approach 

Here, we will obtain an expression for the flux in terms of the efferiive uiofle-to-mode 
energy transfer. Since "^““(klplq) represents energy transfer from p t(/ k witlj q as the 
mediator, we may alternatively write the energy loss from a sphere m 

n(A') = - E E ^““(k|p|q)- < 3 . 29 ! 

M<K |p|>A' 

where mode k is inside the sphere, p is outside the sphere, and q = -k - p. 'Fhe niodf*- 
to-mode transfer '|?.““(k|p|q) consists of a circulating part and an effective {)art. i'Vom Fig. 
3.6 we see that the net amount of circulating transfer leaving the sphere is zero. Thu.s the 
circulating transfer makes no contribution to the the energy flux from the .'-phere. Therefore. 

it can be removed from Eq. (3.29). Hence, the resultant expression for the flux can be written 
as 


n(A') = - E E )^““(k|p|q). (3.30) 

|kl</C |p|>K 

with a summation over the effective mode-to-mode transfer. As expected, the expressions in 
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Figure 3.6: The circulating transfer Xa does not contribute to the energy flux out of the 
sphere K. 


Eqs. (3.28) and (3.30) are equivalent as proved in Appendix B. 


3.6 Energy transfer in a MHD triad: known results 


The MHD equations in real space are written as 


^ + (u; V)u = -Vp + (b. V)b + i/V’u. (3.31 ) 

at 

and 

^ + (u. V)b = (b. V)u + (3.32) 

where u and b are the velocity and magnetic fields respectively, and z/ and p are the fluid 
kinematic viscosity and magnetic diffusivity, respectively. In Fourier space, the kinetic energ}’ 
and magnetic energy evolution equations for a Fourier mode are 


3.6 Energy transfer in a MHD triad: known results 


60 


^E^ + 2uk^E-{k)= E i5““{k|p,q)+ E 

dt !c+p-n=:a- 




gg^(k) 

dt 


+ 2nk^E\k) 


^ is“(k|p.q)+ E :i.;t4, 


k+P+q=0 “ k+p-^q=0 

where £^“(k) = |u(k)|V2 is the kinetic energy, and ^^(k) = bik::* 2 is 'ii»‘ :iiau;.<-’ .• 
The four nonlinear terms 5““(k|p,q), ^“‘(klp.q) 5‘*(klp.q) an<i k p. c}- aja 




5“’^(k|p, q) = -Re (f[k.u(q)][u(k).u(p)] + ifk.u. p ->( k !.u q . . :l;iT 

5'’^(klp,q) = -i?e(i[k.u(q)][b(k).b(p)] + i[k.u|p}’|b|k‘;,biq!’ . 

5“*’(k|p,q) = (j[k.b(q)][u(k).b(p)] + i[k.b{pi][ui k ).bi q r ' . lil.'JT) 

5'*’“(k|p,q) = i?e(?[k.b(q)][b(k).u(p)] + i[k.b{p)l[biki.tnqf'i . 

These terms are conventionally taken to represent the nonlinear tran.sfer fruni inode- p and q 
to mode k [36, 44]. Note that the wavenumber triad k, p, and q should saH.sfv the fuiidition 
k + p + q = 0. The term ^““(klp, q) represents the net transfer of kinetic eiu*rgy from mode.s 
p and q to mode k. Likewise the term 5“*(k|p, q) is the net magnetic er«*rgy I raiisferred from 
modes p and q to the kinetic energy in mode k, whereas 5*’“{k|p. q) is the net kinetic energy 
transferred from modes p and q to the magnetic energy in mode k. The term .s'*'*’(kip. q) 
represents the transfer of magnetic energy from modes p and q to mode k. .411 these transfer 
terms are represented in the Fig. 3.7 shown below. 

Stanisic [36] showed that the nonlinear terms satisfy the following detailed conservation 
properties: 


5“(k|p, q) + S“(p|k, q) + 5"“(q|k, p) = 0. 


(3.391 
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Figure 3.7: The combined energy transfers in a triad, ^““(klp. q): u(p) and u(q) -> u(k]; 
5*“(q|k, p): u(k) and u(p) b(q); 5*’*’(k|p,q): b(p) and b(q) ->• b(k). 


5“(k|p, q) + 5“(p|k, q) + 5“(q|k. p) = 0. (3.40) 


and 

S’-'lklp, q) + 5"‘(p|k, q) + 5“‘(q|k, p) +5‘”(k|p, q) + 5‘“(p|k, q) + 5‘"(q|k, p) = 0. (3.41) 

The Eq. (3.39) implies that kinetic energy is transferred conservatively between the velocity 
modes of a wavenumber triad, and the Eq. (3.40) implies that magnetic energy is also trans- 
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ferred conservatively between the magnetic modes of a wa'.cr.i. tria.l, I h(|. ( 34 ]j 
implies that the cross transfers of kinetic and magnetic > ik p. ({ ^ ■‘ik p.qj 

within a triad are also energy conserving. 

The quantities 5““(k|p,q), 5“‘(klp,q). 5^^(kip.qk atnl ^'‘■'lk p. tj rrpr-' -i.: tlr non- 
linear energy transfer from the two modes p and q to mode k. It v.i»vh''. !'»' !■;; tu know 

the energv transfer between any two modes, say from p to k. Ihi> vv;,t Uf ' ct*.. ;j-, 

the next section. 


3.7 “Mode-to-Mode” energy transfers in MHD Etitia- 
tions 


In Section 3.4 we had presented a new approach to describe the energy tra^^fer between two 
modes of a triad with the help of the third mode, in the Navier-Stokes eqn.i*;* n; In the same 
spirit we will now attempt to find the mode-to-mode energy transfers in MHD. in MHD 
turbulence there will be three kinds of mode-to-mode transfer within a wavenumber triad 
k, p, q: kinetic energy transfer from u(p) to u(k); magnetic energy tran.sfer fnnn b|p) to 
b(k); and transfer of kinetic energy from u(p) to magnetic energy in b(k). We deiH)te these 
three transfers by "^““(klplq), 1Z.**’(k|pjq), and ■^*“(k|p|q) rosperlivelv, where the index q 
indicates that the energy transfer between modes p and k is mediated by the mcxie q. .Since 
the nonlinear interactions fundamentally involve three-modes, the energy transfer between a 
pair of modes should depend on the third mode. The mode-to-mode transfer.s and 

are schematically illustrated in the figure below. 

In this section we will obtain a description for each of these transfers, 

3.7.1 Velocity mode to velocity mode energy transfers 

In Section 3.4 we discussed the mode-to-mode transfer, between velocity mode to ve- 
locity mode in fluid flows. In this section we will find 7?*“ for MHD flows. The transfer 
of kinetic energy between the velocity modes is brought about by the term (u.V)u both in 
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u(q) 




Figure 3.8: The mode-to-mode energy transfers sought to be determined. 

the Navier Stokes equation and MHD equations. Hence, the expression for the combined 
kinetic energy transfer to a mode from the other two modes of the triad is also same for the 
two. i.e., the combined transfer to u(k) from u(p) and u(q) is given by 5’““(k|p, q) [see 
Eq. (3.3) and Eq. (3.35)]. Consequently, ’s for MHD will satisfy the constraints given 
in Eqs. (3.6)-(3.11) for the corresponding /?““ ’s for fluids. As a result. /?““(k|p,q) in MHD 
can be e.xpressed as a sum of a circulating transfer and the effective transfer ^““(k|p|q) 
given by Eq. (3.12), i.e.. 


r“(klp|q) = r“(k|p|q)+X^ 


(3.42) 
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A detailed discussion can be found in the earlier Section 3.4.2. 

3.7.2 Magnetic mode to Magnetic mode energy transfers 


This section contains a discussion on the magnetic energy tran>ff'r ffin i 
in the triad (k,p,q) [see Fig. 3.8]. This transfer which is dfiuitfi in 
should satisfy the following relationships (refer to Fig. 3.9 for a piCni ..i! ^ 

relationships): 

Ko.lo b(p) to b(k) 

p q , which 

■' ' ■ ‘y*’* of ’hpse 

1. The sum of the mode-to-mode energy transfers from b(p) to b|'k) 

b(k) should be equal to 5^*’(k|p, q) [Eq. (3.36) in .Section 3.6). wiiii 

energy transfer to b(k) from b(p) and b(q). Thus. 

?t“(k|p|q) + ?^“(k|q|p)=5“(k.p.qi, 

.in*; fru:;j bfq) to 

*3.4:}) 

1t“{p|k|q) + ?t“(p|q|fc) = .S'“(pik,q), 

t3.44) 

1f“(q|k|p) + ?t“(q|p|k) = 5“(q|k. p). 

13.45) 

and 


2. the energy transfer from b(k) to b(p), ^*’^(k|p|q). and the transfer from b(p > tu b(k). 
^^*’(p|k|q). should be equal but opposite in sign. 

^ V|p|q) + ie%|k|q) = 0, (:U6) 

^''(klq|p) + 1i‘'(qlk|p) = 0, 

(3.47) 

^%|q(k)-F^%|p|k) = 0. 

(3.48) 


The above equations cannot uniquely determine 1t“(k(p|q) since the value of the do 
terminant formed from these equations is zero. While discussing the mode-to-mode energ}- 
transfers in Navier Stokes equation in Section 3.4 we got the same result for ^“(klp|q). 
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u(q) 





b(k) b(p) 



Figure 3.9: Relationship between the Transfer b(k) b(p) is equal and opposite 

to the transfer b{p) — >• b(k) since both are physically the same. Sum of transfers from 
b(q) -4 b(k) and b(p) -4 b(k) is equal to the combined transfer. 


Following the reasoning in that section, we can extract physically meaningful information 
from Eqs. (3.43)-(3.48). 

The combined energy transfer to b(k) from b(p) and b(q) is given by Eq. (3.36) of 
Section 3.6. We denote the first term on the right hand side of that equation by ^^*’(k|p|q). 


i.e.. 


5 ?'’^(k|p|q) = -Re(f[k.u(q)][b(k).b(p)]) , 


(3.49) 
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By replying ^(klplq) 1t“(k|p|q) in Eqs. (3.43)-(3.48) we fin<l tha. JT' ,k , 

solution of these equations, i.e., 

5 ?^^(k|p|q) + ^**(k|qlp) = ■S'**(i<|p. q). 




^''(p|k|q) + ^^^p|q|k) = 5''(pjk. q) 




,y“(q|k|p)+Aq|p|k) = 5“(qik,p!. 


( *152! 


Ak|p|q) + Ap|klq) = 0, '.{.1:11 

Ak|q|p)+Aq|k|p) = 0. i.5.14i 

Ap|q|k) + ,?''(q|p|k) = 0. i:U>i 


Thus, 5 ?^*’(k[p|q) is a solution of the equations (3.43)- (3.48). By insp*‘rtii)a ii i-an be st'on 


that all solutions of the equations can be expressed as 

1t''(klp|q) = ^“(k|p|q) + K^, 

^>|q|p) = ^^‘(k!q|p) - Vk, CUT) 

^%|k|q) = ^ Vikiq) - 

1^^‘(p|q|k) = i?^(p|q|k) -b Ka, !3..)9) 

% (q|k|p) = 5? (q|k|p) -f i",!, (3.()0) 

1^*^q|p|k) = ^?^^(q|p|k) - (.3.61) 


where is an arbitrary scalar function dependent on the wavenumber triad k, p, q and 
the Fourier components u(k), u(p), u(q), b(k), b(p), b(q) at those wavenumbers. To 
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determine the unique form of this function we need to impose additional constraints, similar to 
those used in Appendix A for determining a similar function of mode-to-mode energy transfer 
between velocity modes in fluid turbulence (see Section 3.4.2). We will impose symmetry, 
invariance, and some other constraints to determine I'a (see Appendix B for derivation). The 
exact form of Va is however not crucial to our formalism. 


u(q) 



Figure 3.10: Energy transfer between a pair of magnetic modes in a triad can be expressed 
as a sum of circulating transfer and effective mode-to-mode transfer. 
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The solutions in Eqs. (3.56)-(3.61) have been schematira!!} 
see from the figure that Ia is transferred from b(p) — t b'ki b'tj 
i.e., it circulates around the triad without causing any change ni 
a circulating transfer. The magnetic energy effectively 
^^'’(k|p|q)Si.e., 


in Fie. 3.10. We 

■* 4j|*; K tij 

:* ;> a thuj 
; b{pl *.1 b(ki !s just 






The effective mode-to-mode transfer ^**(kjp|q) from b{p) to b{k) njeefj, tj.,. 
mode u(q). 

In the next section we will use to calculate effertivt* enemy trioi'-fej- i.etwfvn 
the magnetic modes inside two shells and the effective cascade rate uf :'.,i r * ' * cni'tep 

3.7.3 Velocity mode to Magnetic mode energy transfers 

In Section 3.7.1 and Section 3.7.2 we discussed mode-to-mode energy tran-fer between a pair 
of velocity modes and between a pair of magnetic modes in a wa\eniin'.leT triad, v nc. ,.!• 
In the same spirit, this section is devoted to finding the energy *rar:-.f-T lietwren the u(p) 
and b(k), '^“*(k|p|q), within the triad (k, p, q) illustrateil in Fig. 3.S. VW* ;vii: fuliiav tlie 
same sequence of steps as in the two sections 3.7.1 and 3.7.2. 

’s will satisfy the following relationships: (see Fig. 3.11 for an illustnition of these 
relationships): 


1. Since % (k|p|q) and %. (k|q|p) are the mode-to-mode energy transfers from b(p) to 
u(k) and from b(q) to u(k), the sum of the two should be equal to .S-'‘''(kjp. qj. the 
combined energy transfer to u(k) from b(p) and b{q). Therefore, we get 

1J“‘(k|p|q) -f :^“‘(k|q|p) = 5“*(k|p, q). (3.63) 


^“'(P|k|q) +K'(p|qlk) = 5’"*(plk,q), 


(3.64) 
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■!^“‘(q|k|p) + ?E“‘(q|p|k) = S“‘(q|k, p). (3.65) 

?:‘“(k|p|q) + ?t‘“(k|q|p) = 5‘“(k|p. q), (3.66) 

?j‘"(p|k|q) + ?e‘"(p|q|k) = S‘“(p|k. q), (3.67) 

?‘"(q|k|p) + ie‘"(q|p|k) = 5‘”(q|k, p). (3.6S) 

2. 'i^“*’(k|p|q) denotes the energy transfer from b(p) to u(k). '!^^“(p|klq) is physically 
the same transfer but seen as a transfer from u(k) to b(p). Therefore, 

r‘(k|p|q) + 1t‘”(p|k|q) = 0. (3.69) 

r‘(k|q|p) + 1J‘”(q|k|p) = 0, (3.70) 

r‘(p|q|k) + 1t*"(q|pik) = 0. (3.71) 

?t‘*(k|p|q) + r‘(p|k|q) = 0, (3.72) 

y‘“(k|q|p) + r‘(q|k|p) = 0, (3.73) 

•!t‘“(p|q|k) + r‘(q|p|k) = 0. (3.74) 


The solutions of these equations are not unique and the expression for the energy transfer 
between modes cannot be obtained from the above equations alone. We have encountered 
this situation earlier while discussing energy transfer between a pair of velocity modes in 
Sections 3.4 and 3.7.1, and a pair of magnetic modes in Section 3.7.2. Like in those sections, 
we will now explore the solutions of the above equations. 

W'e define the following quantities : 

r‘(k|p|q) = -i(e(i[k.b(q)][u(k).b(p)l) , (3.75) 


)?‘“(k|p|q) = -fle(ilk.b(q)l[b(k).u(p)]), 


(3.76) 
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Figure 3.11: Relationship between the Energy transfer fronj biq) - utk) i> 

physically the same as the transfer u(k) -f b(q). Sum of transfers from ui kl -* biqj and 
u(p) b(q) is equal to the combined transfer. 

Replacing and 5 ?*’“ for and respectively in Eqs. (3.63)-(3.74) we find that 
are a solution to those equations, i.e., 

r'(k|p|q) 4-r'(kiqlp) = 5“'(k|p,q), (3.77) 


5^“^P|k|q) + 5?“^(piq|k) = ■S“^(p(k,q), 


(3.78) 
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r‘(qik|p) +r‘(q|pM = S"‘(q|k,p). 

(3.79) 

^‘•(k|p|q)+^‘-(k|q|p) = 5‘*(k|p,q), 

(3.80) 

^‘•(p|k|q)+^“(p|q|k) = 5‘*(p|k.q), 

(3.81) 

^‘“(q|k|p) + ^““(qlplk) = 5*“(q|k, p). 

(3.82) 

and 


^“‘(k|p|q) + ^‘“(p|k|q) = 0, 

(3.83) 

r‘(k|q|p)+^(q|k|p) = 0, 

(3.84) 

^•‘(plqlk) + ^‘“(q|p|k) = 0, 

(3.85) 

^‘"(k|p|q)+r‘(p|kfq) = 0. 

(3.86) 

^‘“(k|q|p)+^“‘(q|k|p) = 0, 

(3.87) 

^‘“(p|q|k) + r‘(q|p|k) = 0. 

(3.88) 

The 's are just a single instance of the the "^“^’s. It can be seen 

solutions can be expressed in the form : 

by inspection that all 

:!J‘“(k|p|q) = ^‘"(k|p|q) + Za. 

(3.89) 

?t‘“(k|q|p) = Ak|q|p)-ZA. 

(3.90) 

^““(plklq) = /“(p|k|q) - Za, 

(3.91) 

?^‘“(p|q|k) = ^(p|q|k) + Za, 

(3.92) 


?t‘“(q|k|p) = /“(q|k|p) + Zi, 


(3.93) 




._J2 
^ 3 , 94 ! 
' 3 . 95 } 

.' 3 .^ 19 * 

.i.ioin 


where Za is an arbitrary function dependent on k, p, q, u(k), u(p}. u(q), b{k), b(p). 
b(q). These solutions are pictorially represented in Fig. 3.12 below. 

We are already familiar with such solutions from the earlier disfUhvjon on nu>d»*‘to-nu)(ip 
transfers between the velocity modes and between the magnetic b'rorii b’ig. .3.12 we 

see that Za transfers energy from u(p) b(k) ->■ u(q) b(pl -4' u( k ! -» bitj i -* and bark 
to u(p) without resulting in a change in modal energy. Hence, follovviu;- the diMusMoib in 
Section 3.7.1 and Section 3.7.2, we can interpret the quantit} Za a.^ a circjU.i'ing transfer. 
The can be interpreted as the effective mode-to-mode transfers. For «'\air!j)!<% ^^“{kipiq) 
is the effective transfer from u(p) to b(k),i.e. 

?S/(k|p|q) = ^(k|p|q), (:U01) 

which is mediated by the magnetic mode b(q). 

In the next section, we will use the effective mode-to-mode transfer to define, (a) shell-to- 
shell transfers between the velocity and the magnetic modes, and (b) cascade rates between 
the velocity and the magnetic modes. 
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u(q) 



Figure 3.12; Energy transfer between a velocity and a magnetic mode in a triad can be 
expressed as a sum of circulating transfer and effective mode-to-mode transfer. 


3.8 Shell-to-Shell energy transfer and cascade rates in 
MHD turbulence 


In Section 3.5 we had defined the effective shell-to-shell transfer and the cascade rates in fluid 
turbulence in terms of the effective mode-to-mode transfer rate between a pair of modes in 
a wavenumber triad. In this section, we will extend the formulation of Section 3.5 to define 
shell-to-shell transfer rates and cascade rates in MHD turbulence. 

In this thesis we will use the term w-shell (6-shell) to represent a shell in ^-space con- 
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taining velocity (magnetic) modes. The term tr-sphere (^-sphere^ will tu 

sphere in )fc-space containing velocity (magnetic) modes. Thu.s. th(' en<*rg> 
u-sphere (6-sphere) is kinetic (magnetic) energy. 


with a 


3.8.1 Shell-to-Shell energy transfer rates 

In Section 3.5 we had defined the effective shell-to-shell tran.sfcr rates !>efwcen t*.vu •; '.in'Il.*. 
Fourier space. The expression for the effective mode-to-niodc transfer he ueej; j-au ’.(•;(» ity 
modes is the same for both Fluid and MHD flows [see Section 3. 1.J and 3. . . ; lieui •*, ;|j<. 
effective shell-to-shell transfer rates between two u-shelis in MH!) wi;I he de!;:;fr| |)\ 
Eq. (3.27) of Section 3.5, i.e., 

^rn = EEr“{kip|q). -3.102. 

kcm 


The quantity '^*^(k|p|q) gives the energy transfer rate from b(p) to b(k) ineriiaterl hy 
q. Hence the energy transfer rate from 6-sheIl to the fr-shell can he obtained In 
summing "^^^klpjq) over every p in the the 6-sheII and over every k in the m*'* 6-shell, 
From Section 3.7.2 we know that = ^**’(k|p|q) •+ 1^, where ^^(kjplqi is the effective 
mode-to-mode transfer from b(p) to b(k), and is a circulating transfer. Thus. > j is 
transferred from b(p) in shell n to b(k) in shell m and then back to b(p) via the niude 
b(q). Therefore, the effective shell-to-shell transfer transfer rate from the /i*'*' to the 
6-shell is 

= E Et^^^’Wplq)- (3.1031 

kcm pcTi 

We can also define shell-to-shell transfers rates between a w-shell and a 6-shell. k jp!q) 
is the energy transfer rate from u(p) to b(k) mediated by mode q. Hence the energy transfer 
rate from the the n** a-shell to the 6-shell can be calculated by summing ^*“(klplq) over 
every p in the n w-shell and over every k in m** 6-shell. We have shown in Section 3.7.3 
that (kjpjq) — ^ (kjpjq) -|- Z^, where '^*’“(k|p|q) is the effective mode-to-mode transfer 



3.8 5hell-to-Shell energy transfer and cascade rates in MHD turbulence 


75 


from the mode u(p) to the mode b(k) in the triad k, p, q, and Za is the circulating transfer 
which is transferred along u(p) -> b(k) u(q) — )• b(p) — >• u(k) -4 b(q) -4 u(p). Hence, 
Za is transferred from the mode u(p) in shell n to b(k) in shell m and then flows back to 
u(p) via the modes u(q) -4 b(p) -4 u(k) -> b(q) -4 u(p) -4 u(p). Therefore, we interpret 
the quantity obtained by summing j?*“(k|p|q) over every p in the u-shell and every k in 
the 6-shell as the effective shell-to-shell transfer, i.e., 

(3.1041 

kcm pen 

is the effective transfer from the the u-shell to the 6-shelL The energy transfer rate 
from the 6-shell to the u-shell, = —T^^. 

In Appendix A-C, we have argued that the circulating transfer = 0, I'a = 0. and 
Za = 0 to the first order in the expansion coeflBcients. If ATa = 0, Ya = 0, Za = 0. then the 
effective shell-to-shell transfers will be the same as the total shell-to-shell transfers. We had 
pointed out is Section 3.5 that the circulating transfer, if nonzero, does not contribute to a 
change in energy of any shell, and hence from the physical point of view the effective shell- 
to-shell transfers is of greater significance than the total shell-to-shell transfers. We have 
numerically computed the effective shell-to-shell transfers given by Eqs. (3.102)-(3.104) in 
numerical simulations of 2-D MHD turbulence and gained important insights into the energ}- 
transfer process. The results of the simulations will be discussed in the next chapter. 

3.8.2 Energy cascade rates 

There are various types of cascade rates (energy fluxes) in MHD turbulence. We ha\-e 
schematically shown these transfers in Fig. 3.13. In this section we will derive the formulae 
to calculate these cascade rates within the framework of effective transfers. 

In Section 3.5.1, we derived an expression for kinetic energy flux in terms of the mode-to- 
mode transfer [Eq. (3.26)]. We showed that the circulating transfer Xa does not contribute 
to the flux. Hence, the kinetic energy flux can be expressed by Eq. (3.30) in terms of the 
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effective mode-to-mode transfer ^““(k|p|q). The kinetic energy flux 11“^ (the energy transfer 
rate from the w-sphere to outside the same sphere) in MHD will also be given by the same 
equation, i.e., 

n;;(A-) = E E r“(klp|q). (3-105) 

IpKS- |k|>K 

The magnetic energy flux nj>(A') is defined as the rate of energy lost by the Asphere to 
the modes outside the 6-sphere. Since ■^*’^(k|p|q) is the mode-to-mode transfer from b(p) to 
b(k).ng(A') can be obtained by summing "^^^(klplq) over every mode p inside the 6-sphere 
and every mode k outside the 6-sphere, i.e., 

= E E 1^“(k|p|q)- (3-106) 

|pi<A'|k|>if 

Like for the kinetic energy flux, the circulating transfer will not contribute to the magnetic 
energy flux in Eq. (3.106). We briefly explain this as follows: let us take the mode b(q) to be 

outside the 6-sphere. Then if the mode b(k) loses Ya to mode b(p) then it will also gain La 

from mode b(q), and hence Ya will not contribute to the magnetic energy flux. Therefore, 
we can write the magnetic energy flux in Eq. (3.106) as 

ng(A')= E E i^''(k|p|q)- (3-i0') 

IpKA |k|>A' 

where ^^*(k|p|q) is given by Eq. (3.49) of Section 3.7.2. 

The energy transfer between the kinetic energy and magnetic energy can be described 
by four types of fluxes. We shall define these fluxes below. 

There is a transfer of energy from a w-sphere of radius A' to the 6-sphere of the same 
radius. The rate of loss of energy from the u-sphere to the corresponding 6-sphere can be 
calculated by by summing "^‘“(klplq) over every mode inside the 6-sphere and the u-sphere. 
Following the definition of the effective shell-to-shell transfer we can define the effective 
flux from the u-sphere to the 6-sphere as 

= E E ^^“(k|p|q)- (3.108) 

|p|<A |k|<A 


3 8 Shell-to-Shell energy transfer and cascade rates in 
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by summing the effective mode-to-mode transfer. Ju'Hi mp) b!ki, 

modes inside the u-sphere and the 6-sphere. 

There is a transfer of energy from the a-sphere to th<> magne;!.- -t.r . ,<rr<v 

ponding i^sphere. The rate of loss of the the energ}’ from the f.- :ie n,. 

the 6-sphere can be obtained by summing ie'’“(kipiq) over every iiunie 

and every mode outside the 6-sphere. The corresponding effective- flux ran . .ii. uhitrd 

using 

n,y(A')= E E ^‘"(klpiqi- ■i nw. 

lp|<A- |k|>A' 

Similarly, there is a transfer of magnetic energy from a u> whtruy ji.ude^ 

the u-sphere. The rate of loss of energy from the 6-.sphere to the nuHie> unlade the «-spi»w 
can be calculated by 

Tiib/ 


nh(A’) = E E r‘(k|p|q). 

|p|<A |k|>A 


(•■f.llO) 


There is a transfer of energy from modes outside the H-.s])h<*r«‘ to the ino(i«*s o:its!<ii* the 
6-sphere. The rate of loss of energy by the modes outside the w-sphere tti tin* iiuide*' (ditside 
the 6-sphere can be obtained by 

n,t(A')= E E ^'(k|p|q)- 

|p 1>A’ |kl>A^ 

The total effective flux is defined as the total energy (kinetir-f mugiu-ti* ) lost by the 
A'-sphere to the modes outside, i.e., 

ntot(A') = n:<(A0 -b ng(A') + ut<(A') + ni<(K). cunb 


A schematic illustration of the effective fluxes defined in Eqs. (3.105)-{3.111) was given 
in Fig. 3.13. If the circulating transfer is zero, then each of the effective fluxes will be the 
same as the corresponding total fluxes. In the next chapter we will present the results of the 
detailed study of these fluxes in numerical simulations of 2-D MHD turbulence. 
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3.9 Conclusion 

In literature we find a description of energy transfer rates from two modes in a triad to the 
third mode, i.e., from u(p) and u(q) to u(k). Here we have constructed new formulae to 
describe energy transfer rates between a pair of modes in the triad (mode-to-mode transfer), 
say from u(p) to u(k). The third mode in the triad acts as a mediator in the transfer 
process. 

The mode-to-mode energy transfer in our formalism can be expressed as a combination 
of an “effective transfer” and a “circulating transfer”. In Appendices A-C we have shown 
by imposing galilean invariance, symmetry and some other constraints that the circulating 
transfer should be zero. A completely general proof is not available at present. Howe\'er 
even if the circulating transfer happens to be nonzero, it will not result in a change of modal 
energy, since the amount of circulating transfer gained by mode k from mode p is also lost 
by mode k to mode q (see Figs. 3.4, 3.10, 3.12). Only the effective transfer is responsible 
for modal energy change. As the circulating transfer does not have any observable effect on 
the energy of the modes, it may be correct to ignore it from the study of energy transfer. 

Using the notion of effective mode-to-mode transfer and the circulating transfer we 
defined effective shell-to-shell energy transfer rates and energy fluxes in fluid turbulence [see 
Eq. (3.27) and Eq. (3.30)] and in MHD turbulence [see Eqs. (3.102)-(3.104) and Eqs. (3.105)- 
(3.111)]. Some of these energy transfers can be obtained only using the “mode-to-mode" en- 
ergy transfer formulae. In the next chapter we will present the results of our numerical study 
of the effective shell-to-shell energy transfer rates and energy fluxes in 2-D MHD turbulence. 




Chapter 4 

Energy transfers in two-dimensional 
magnetohydro dynamic turbulence 

4.1 Introduction 

Turbulence is characterised by the existence of a range of spatial scales. Energy is transferred 
between these scales. Some of the features of energy transfer arc well c.stabli.slH'd. while a 
few other have been a matter of controversy, and some others have not been investigated at 
all. A detailed account of these issues was given in Chapter 1. 

The following features of energy transfer in MHD turbulence ha\'e been conclusively 
established. The total energy (kinetic+magnetic) cascades from the large scales to the small 
scales in both 2-D and 3-D MHD turbulence [11]. Fluxes associated with r"*" and energies 
are also found to cascade to the small scales [11]. 

We know from our earlier discussion that the magnetic energy of a mode e\'olves due 
TO two nonlinear terms in the MHD equations — [b.(u.V)b] and [b.(b.V)u] — the first 
exchanges magnetic energy between different scales of the velocity field, and the second 
exchanges magnetic and kinetic energy between different scales. The kinetic energy similarlj 
evolves due to two nonlinear terms — [u.(u.V)u] and [u.(b.V)b] — the first one exchanges 
kinetic energy between different scales and the second exchanges energy between the magnetic 
and the velocity fields. The kinetic and magnetic energy transfers arising from each of these 
non-linear terms above have not been investigated in detail either analytically or numerically 
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till very recently [20]. It is the purpose of this chapter to do a detailed investigation of energy 
transfer in 2-D MHD turbulence. To this end we will numerically compute the effectix'e 
energy fluxes and the effective shell-to-shell energy transfer rates using Eqs. (3.105)-(3.111) 
and Eqs. (3.102)-(3.104) of Chapter 3. Recently, Ishizawa and Hattori [20] studied some of 
the energy transfer properties in 2-D MHD turbulence using wavelets. They computed the 
net energy transfer to a wavelet from all other wavelets. In contrast, we have calculated 
the energy cascade rates and detailed shell-to-shell energy transfers. While discussing our 
simulation results we will point out results common to both our simulations. 

Non-local aspects of energy transfer in 2-D and 3-D MHD turbulence have been studied 
using EDQNM closure calculations [17-19]. For 2-D MHD turbulence, EDQNM calculations 
of both Pouquet [18], and Ishizawa and Hattori [19] give a transfer of kinetic energy from the 
small-scale to large-scale velocity field. Pouquet [18] also predicted magnetic energj' transfer 
from small-scale magnetic field to large-scale magnetic field, and another transfer from large- 
scale magnetic field to small-scale velocity field. However, the conclusions of Ishizawa and 
Hattori are exactly the opposite [19], and are also confirmed by means of a direct numerical 
simulation [20]. The reeison for the difference in the conclusion of these EDQNM results was 
discussed in Chapter 1 (see Section 1.5). Pouquet et al. [17] have also studied 3-D MHD 
turbulence using EDQNM closure calculation, as described in Chapter 1. 

In some of the studies [17, 21], distinction was not made between the energies transferred 
to a mode from the kinetic and the magnetic fields. Also, the energy transferred into a mode 
k from different wave number regions were not separately considered [20, 21] — only the 
net energy transfer into a wave number k is generally computed. The EDQNM closure 
calculations [17-19] dealt mainly with coarse-grained energy transfer (between large scales 
and small scales). In another study, Frick and Sokoloff [77] solved a shell model of MHD 
turbulence and calculated the kinetic energy flux resulting from the interaction of the velocity 
modes, and magnetic energy flux arising from the interaction of the magnetic modes. In 
our simulations we investigate the (1) various energy fluxes arising within and between the 
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velocity and magnetic fields, and (2) the fine-grained (considering nianv vsavr shells) 

energy transfer between magnetic and kinetic energies. This vviL give a ihfurnwd picture 
of the physics of energy transfer in MHD turbulence. 

This chapter is organized as follows. In Section 4.2 we will tie.'-i riiie the lamieriral 
technique to compute the fluxes from the velocity to the magnetic field, anti ili«> ■ hcH lu-shell 
energy transfers. The numerical method to solve the MHD ecjUHtion.s ha.*- been 'aabed 
in Chapter 1. In Section 4.3 we will list the parameters for our inimerira! .‘simulation. We 
perform the numerical simulation by forcing the large-scale ve}ocit\' ficdd. I he hteady state 
thus obtained in our simulations is discussed in Section 4.4. The variou.s aiui .shell-to- 

shell energy transfer rates are averaged over this steady state. Re.sults of the- fluxo obtained 
in the simulations is given in Section 4.5. In Section 4.6 we present the r(‘sults for the 
shell-to-shell energy transfer rates. A discussion of the results follows in Section 4.7 


4.2 Numerical computation of fluxes 

In this section we describe the numerical technique used by us to compute the* fluxes in 
Eqs. (3.108)-(3.111) of Chapter 3. To compute the fluxes we employ a metliod similar to that 
used by Doraaradzki and Rogallo [39]. We outline this method below using as an 

example. In Eq. (3.110) of Chapter 3 we substitute the expression for ^*“(kjp|q) [Eq. (3.75) 
of Chapter 3] : 

nt>(A')= -^e(i[k.b(q)][u(k).b(p)]) . (4.1) 

|k|>A- |p|<A- 

A straightforward summation over k and p involves 0{N^) operations, where A' is the size 
of the grid, and would thus involve a prohibitive computational cost for high Re\'nolds num- 
ber simulations. Instead, the pseudo-spectral method can be used to compute Eq. (4.1) in 
0{MogN) operations. It involves the following procedure. 
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We define two ‘truncated’ variables and as follows 

'0 if|kr<A''i 


;(k) ;iiki> A' 


(4.'2) 


and 


b<(p) 


_ / b(p) if IpI < K 
[0 if jpj > K 


(4.3) 


The Eq. (4.1) written in terms of and b"' reads as follows 

nK(A-) = &(.[k.b(k - p)][u>(k).b<(p)i) . 


14.4 i 


k P 


The above equation may be written as 


ng(A-) = Re 


5:.V?(k)X;6,(k-p)6<(p) 


L k 


! 4.-3 I 


The p summation in the equation above can be recognized as a convolution sum. The right 
hand side of Eq. (4.5) can be conveniently and eflficiently evaluated by the pseudo-spectral 
method, using the truncated variables and b"^. This procedure has to be repeated for 
every value of wavenumber K for which the flux needs to be computed. The rest of the 
fluxes defined in Eqs. (3.108)-(3.111) and also the shell-to-shell transfer rates defined in 
Eqs. (3.102)-(3.104) of Chapter 3 are similarly computed. 


4.3 Simulation parameters 

In our numerical simulations, the large-scale velocity field is forced with white-noise forcing. 
The simulations are performed with hyperviscosity and hyper-resistivity characterised by the 
wavenumber keq- The parameters for the simulations are given below: 

• Grid size = 512 x 512. 

• iy = fi = 5x 10 -®. 

• kgq — 14. 

• dt = 5 X lO”"*. 



4.4 Generation of steady state 


• Forced wavenumbers kj : 4<fc/<5. 

• Average total energy input rate ~ 0.1. The f ii: thv function 

[Eqs. (2.16)-(2.17) of Chapter 2] is fi.xed by this rate of <“iierg\ input 

• Magnetic energv' input rate = 0. 


• Steady state values of 


rA 


/u^dx 

Jb^ 


Cii 0.5, 


ac = 


2 f(u.h)dx ^ 
f(u^ + b^)dx' ~ 


iUj 


(4.T) 


Usually the cross-helicity has a tendency to grow. In order to keep tile cros.s-helicity low 

in oui simulations, we had to adjust the injection rate of cross-hclirit\ during the cour.se of 
our simulations. 


4.4 Generation of steady state 


The computational time required to obtain a statistically steady state on a grid of size 512^ 
IS enormous. So to obtain a steady state in simulations on this grid w(* proceed in .stages. 
We run on a grid of size 64^ till a seemingly statistical steady state is achieved. This steady 
state field is then used as the initial condition to achieve a steady state in a .simulation on a 
grid of size 128^ and so on to grids 256^ and then 512^. 


It is theoretically expected that in 2-D, in the long term, the magnetic energy- u-ill decay 
2J even if tie kinetic energy is stationary - this is called the anti-dynamo theorem. In our 


smmlations we found that the fields remain steady for significantly long period after which 
it will decay (Fig. 4.1). Since the N = 512^ simulations are computationally expensive, we 
have instead shown the decay of magnetic energy for a IV = 128* simulation. 
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Figure 4.1: Evolution of the total kinetic energy and the magnetic energy for simulations 
on grid sizes 512^ and 128^. For 512^ a quasi-steady magnetic energy is obtained over the 
period of the simulation. It is demonstrated in the 128^ simulation that a quasi-stead}- 
magnetic energy eventually decays - a quasi-steady magnetic energy is obtained from t = 60 
to 100. 

The magnetic energy is found to have a quasi-steady state over the time interval extending 
from approximately t = 20 — 30 and it decays over the remaining period, as is theoretical!}- 
expected. For the N = 512^ simulation too, after the period of stationarity. the simulations 
should show decay of magnetic energy. So. strictly speaking, our steady state is only a quasi- 
steady state. In Fig. 4.1 we show the magnetic and kinetic energy in a quasi-steady state. 
The Alfven ratio (the ratio of kinetic to magnetic energy) is found to fluctuate between the 
values 0.4 and 0.56. Hence over the steady state, magnetic energy dominates over the kinetic 
energy. The value of normalised cross-helicity is low (close to zero) in the steady state. We 
compute the fluxes and the shell-to-shell transfer rates over this quasi-steady state once in 
every unit of non-dimensional time. The fluxes and shell-to-shell transfer rates shown in this 
chapter are averages over 15 time units. 
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4.5 EnGrgy fluxes in 2-D MUD simulations 


In our numerical simulations we have computed all the energ} Iluxt d* l**i< u iii I.cjs. !,? 
(3.111) of Chapters. These fluxes are schematically shown in Fin. .3.13 of (’haptn .1 bthis 
we section we discuss these energy fluxes (cascade rales). 

In Fig. 4.2 we show all the fluxes. The total flux Iltoi is po.'^iti'e nwit »i!.3ig tUat there is 
a net loss of energy from the A^-sphere to modes outside for all A . in the wavenuiuhcr region 
2-5 < K < 50, the total flux is seen to be approximately constant. This veaw nninber region 


is the inertial range. 

The net transfer from kinetic to magnetic energy is a sum of the fiux<*s II"(^ i A ). 
n 4 >[= (A')], and nj^(A'). In our 512^ simulations, the splien* of radius A„ar = 241 

encloses all the modes. We observe from Fig. 4.2 that the fluxe.s { A'mar ) • i i A n.ar 1 and 

( A'^ax) are zero, as there are no modes outside this sphere of maximum radius A‘„,aj-. The 
flux nj^(A',„aa,) is found to be positive (see Fig. 4.2), indicating that there is a net transfer 
from kinetic energy to magnetic energy; since the magnetic modes are not being forced, it is 
this transfer that keeps the magnetic energ\' constant in quasi-steady .state. 

We now describe features of the fluxes 11“^, 11“^, n“^, and 11^^ observed in 

oui simulations and plotted in Fig. 4.2. We remind the reader that the ii-niode.s within a 
wave niunber sphere are called the u-sphere and the Amodes within the wave number sphere 
are called the Asphere. First, we discuss the fluxes that transfer energy between the u-modes 
and the A-modes. We find that the flux n^^{A') is positive — hence, kinetic energy is lost by 
a n-sphere to the corresponding 6-sphere. The flux 11“^ (A') is also positive. It means that 
a u-sphere loses energy to the modes outside the 6-sphere. We find that the flux ng(A') is 
negative, which implies that the 6-sphere gains energy from modes outside the M-sphere. Thus, 
all these fluxes result in a transfer of kinetic energy to magnetic energy. However, n^>(A') 
is negative, implying that there is some feedback of energy from modes outside the 6-sphere 
to modes outside the u-sphere. The net transfer, however, is from kinetic to magnetic. 
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Figure 4.2: The simulation results of fluxes. The following fluxes have been plotted : the 
total flux (Iltof), the kinetic energy flux from u-sphere to outside w-sphere 11“^, the magnetic 
energy flux from 6-sphere to outside 6-sphere Ilj^, the energy flux from u-sphere to 6-sphere 
the energy flux from u-sphere to modes outside the 6-sphere 11“^, the energy flux from 
6-sphere to modes outside the ^-sphere 11^^, the energy flux from modes outside u-sphere to 
modes outside the 6-sphere , and the net flux out of a 6-sphere -f have been 
plotted in this figure. 

We shall later show that the flux (A') plays a crucial role in driving the kinetic- to-kinetic 
flux n“<(A'). 

The energy gained by the 6-spheres is found to be constant for the approximate range 
20 < A" < Kmax- This can be seen from the plot of n*^(A') + (which is the negative of 
the total transfer of kinetic energy to the modes within the ^sphere). The constancy of the 
flux implies that a 6-sphere of radius K and that of radius of K + AK get the same amount 
of energy from the u-modes. Thus there is no net energy transfer from the a-modes into a 
6-shell (of thickness AK) beyond approximately AT = 20. We therefore conclude that the net 
energy transfer from a-modes to the 6-sphere occurs within the A” < 20 sphere. 
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We find that there is a kinetic energy^ gain by the n-sphrrr from a 
the fact that flux n“<(A') is negative — this is consistent with tiie ntnnrrir.i'i hLJ.-ior.sof 
Ishizawaand Hattori [20]. This behaviour of kinetic energy is kni>wn a,- an ' !n\ rascade' 
in literature and is reminiscent of the inverse cascade of kinetic energy in 2 I ) {bnd tnritulence 
[44] and of mean square vector potential in 2-D MHD turbulence :1S. Tri]. lluwe\'er. note 
that the kinetic energy in 2-D MHD turbulence is not an inviscid invariant, f':p. 4.2 shows 
that the inverse cascade of kinetic energy exists in the wave number range K < 00. Even the 
modes that are being forced (4 < A- < 5) gain energy from higher wave rtuinlwrs u-niodes. 
The source of this energy is the flux n‘^(A')[= -nj>(A')] which transfers etiergy from higher 


Amodes to the higher u-modes and thus effectively forces them. 

We observe from Fig. 4.2 that there is a loss of magnetic energy from lh<‘ h N[)here to 
the Amodes outside [see n*>(A’)] — this was also observed in the numt*rical .simulations 
of Ishizawa and Hattori [20]. Our observ'ation of a forward cascade of magnetic energ,v 
contradicts the result of the EDQNM closure calculations which predicts an inverse cascade 
of magnetic energy [18]. 

The flux nj>(A') is found to be constant in the inertial range. Using similar leasoning 
to that given above, we can conclude that the net magnetic energy transfer to a wave number 
shell m the inertial range is zero. However, in this case it would imply that the aitin energv 
gained by the shell (of thickness AA) from the Asphere of radius K is lost completely to 
the modes outside the region K -f- AA . Thus the magnetic energy cascades down to higher 
u ave numbers independent of the transfers between kinetic and magnetic energy. This flux 
of magnetic energy is sustained by the fluxes nj:<(A') and n^>(A'), both of which transfer 
kinetic energy into the small Aspheres. 



4.5 Energy fluxes in 2-D MHD simulations 


89 



Figure 4.3: The schematic illustration of the directions and the magnitudes of the flu.xes 
plotted in Fig. 4.2 (also see Fig. 3.13 of Chapter 3). Also shown are the magnitudes of the 
kinetic energy input rate due to forcing, and the total dissipation of kinetic and magnetic 
energy. The fluxes are shown for K = 20 but are representative of the entire inertial range. 
All quantities have been time-averaged. The fluctuations of the fluxes (except ) and the 
dissipation rate are approximately equal to 0.005. The fluctuations in are higher and is 
approximately 0.01 . 
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Figure 4.3 schematically illustrates the energ}' fluxes of F ig. 4.2 for a h ''jiherr of 
K = 20, which is within the inertial range. The energ}’ input <iu<* to for' inc a;ifl jjjjjjj 
inverse cascade [n“^(A')] into the u-sphere from higher a-mode.s. jtro-. ;i!( • tho t'!irri<y input 
into the u-sphere. This energy is transferred into and outside the sphere h) II,’*''' /{'i 
n 4 ;^(A'), the latter transfer being the most significant of ail transfers (see 4/j aj;,i ^ _^j 
The energy transferred into the 6-sphere from the u-sphere [n“^( A')I. and a nechpiitje input 
from the modes outside the u-sphere [— n^^(A')], cascades down to the higher wave mnnber 
6-modes. In the higher 6-modes, this cascaded energy {IIj^(A')]. tugethtT with the transfer 
from the w-sphere, is partly dissipated and partly fed back to the high wa\'e imiuber iMoodes, 
This feedback to the kinetic energy is mostly dissipated, though a small inv('r.‘ie ra.st-ade takes 
some energy back into the w-sphere. This is the qualitative picture of the eiierg\’ transfer in 
2-D MHD turbulence. 

The net transfer to each of the four corners of the Fig. 4.3 sum tc) zero within the 
statistical error (which is computed from the standard deviation of the saiujflrd data). This 
is consistent with a quasi-steady-state picture. 

The results presented here for the quasi-steady-state in a forced turbulence remain qual- 
itatively valid even for a decaying case - the direction of the various fluxes for the decaying 

case are identical to that for the forced simulation; but as the energy decays, the magnitudes 
of all the fluxes reduce. 

The fluxes give us mformatiou about the overall energy transfer from a wave number 

sphere or outside-sphere to another sphere or outside-sphere. To obtain a n.ore detailed 

account of the energy transfer, energy exchange between the wave number shells are now 

studied. In the following section we present a discussion on the sheU-tr^shell energj. transfer 
rates in MHD turbulence. 
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4.6 Shell-to- Shell energy transfer-rate studies 

Significant details of energy transfers are revealed by calculating the shell-to-shell energ}' 
transfer rates T““, , and between different shells defined in Eqs. (3.102)-(3.104) 

of Chapter 3. We partition the k-space into shells at wave numbers kn{n = 1,2.3....) = 
1, 16, 19.02, 22.62, 2^"^'^^^^/^. The first shell extends from Arj = 1 to = 16 — a division 
of the wave number space into smaller shells at the large scales will contain too few modes; 
the second shell extends from k 2 = 16 to k^ = 19.02, ..., the shell extends from km to 
Thus, the effective shell-to-shell energy transfer rate from the u-shell to the nr' 
u-shell [Eq. (3.102) of Chapter 3] can be written as, 

TZ= E E Er“(k|plq). (4.S1 

and the effective shell-to-shell energy transfer rate from the 6-shell to the 6-shell 
[Eq. (3.103) of Chapter 3] can be written as 

r“„= E E El?“(k|p|q). (4.9) 

km<k<km+l kn<P<kn^l Q 

and the effective shell-to-shell energy transfer rate from the u-shell to the rn}^ 6-shell, as 
defined in Eq. (3.104) of Chapter 3, can be written as, 

E E E^‘”(k|p|q). (4.10 1 

km<k<km^l fcn<p</:n+l ^ 

We categorise the energy transfer between u-shells and the 6-shells as being homologous 
if the transfers are between the corresponding shells (of the same wave number range A:„ < 
k < kn+i. say). Transfers between different shells are therefore non-homologous. Further, 
transfers (kinetic-to-magnetic, kinetic-to-kinetic, and magnetic-to-magnetic) involving shells 
which axe close in wave number space are caJled local transfers, as is the convention. Two 
shells are considered close if the wave number ratio of the larger shell to the smaller shell is 
less than 2 (in this study this would include the shells between n 4- 4 to n — 4 from the 
shell). Transfers involving shells more distant are called non-local. 



4 6 Shell-to-Shell energy transfer-rate studies 


92 



In Fig. 4.4 we plot the energy transfer rates between u-shells and fr-shells. It is 


evident from the figure that the transfer rates between shells in the inertial rangc‘ are virtually 
independent of the individual values of the indices m and n, and only dependent on their 
differences. This means that the transfer rates in the inertial range are self-mmilar. The 
differences in 2^“ for various n are smaller than the standard deviation of the saniph>d data, 
indicating that the perceived self-similarity is statistically significant. 

discuss our simulation results for the non-homologous transfer rate.s between 
^-shells and ^^shells. In Fig. 4.4 we have shown these transfer rates from the n** u-shell to 
the m ^shell by plotting versus m for various values of n. We find that for all m. 


except for m n 1 and n, is positive. This implies that a w-shell loses energy' to all 
the ^.shells but gains energy from the (n - and fr-shells. The quantity is found 
small for m < n m comparison with the transfer for m > n. Consequently, energy 
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from a u-shell is mainly transferred to the ^shells at higher wave numbers. We had found 
earlier that the flu.x is much smaller in magnitude to the flux 11“;^ (/\). The flu.xes 

and the transfer rates are consistent with each other. From Fig. 4.4 we also find that the 
energy transfer rate from the u-shell to the 6-shell is positive, implying that the 

6-shell gains energy from the u-shells — this is a nonlocal transfer of energy from the 
large u-modes to the small 6-modes. The inverse cascade observed in Section 4.5 (see 
Fig. 4.2) is due to this nonlocal transfer. 

The homologous transfers are found to transfer energy from 6-shells to u-shells (see 
Fig. 4.4); this is in contrast to the net energy transfer and a majority of the non-homologoiis 
transfers which are from kinetic to magnetic. We find that the energy gained by a u-shell 
through the homologous transfers is larger than the total loss of energy by the u-shell through 
non- homologous transfers. Consequently, there is a net gain of energy by the u-shells in the 
inertial range. In Section 4.5 it was shown that the magnetic energy outside a 6-sphere is lost 
to kinetic energy outside the u-sphere. It is now clear that this transfer arises primarily due 
to the homologous transfers from magnetic to kinetic. 

The Fig. 4.4 shows energy transfer rates to 6-shells from u-shells indexed as n = 
5. 6,7, 8.9. In Fig. 4.5 we show the transfer rates from the n = 1 u-shell to all the 
6-shells. Recall that the first u-shell is comprised of small wave number modes, k = I to 16. 
Comparing the magnitudes of to the 6-shells from the first u-shell (see Fig. 4.5) and from 
the u-shells at higher wavenumbers (see Fig. 4.4), we see that the energy transfer rate from 
the u-shell dominates the transfers from all other u-shells. Hence, there is a large amount 
of non-local transfer from the first u-shell to the 6-shells. In Section 4.5 we had claimed 
that there is no net kinetic energy transfer into any inertial range 6-shell. Now we show this 
explicitly in Fig. 4.6 by plotting the net kinetic energy transferred into a 6-shell {=J2n ^mn) 
versus the shell index m (o in the figure). From the figure, the net energy transferred into the 



4.6 5hell-to-5hell energy transfer-rate studies 


95 


0.04 
0.035 
0.03 
0.025 
0.02 

Vt’ bu 

0.015 
0.01 

0.005 
0 

-0.005 
- 0.01 

0 2 4 6 8 10 12 14 1< 

m 

Figure 4.6: The diamonds (o) represent the net energy transfer into a i>-sheli from all the 
u-shells. The pluses (-f) represent the net energy transfer into a 6-shell from all u-shells 
except the 1®* one. 

inertial range 6-shell can be seen to be nearly zero. We also plot in Fig. 4.6 the net kinetic 
energy transfer into a 6-shell from all w-shells, except the first (-f in the figure). This quantity 
is now no longer zero and has a significant magnitude. Therefore, the non-local transfer from 
the first shell plays an important part in balancing the other kinetic energy transfers into an 
inertial range 6-shell. 

We have also computed the transfer rates of magnetic energy from to 6- 
shell. In Fig. 4.7 we have plotted the quantity versus m for different values of n in the 
inertial range. We find that the differences in T^n+An)n fixed An(=n — m) are smaller 

than the standard deviation of their statistical fluctuations computed from the sampled data. 
Hence we can conclude that r,^*’„'is self-similar in the inertial range, dependent only on the 
difference An and independent of the location of the shell n. We find that the transfer rates 
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Figure 4.7: The energy transfer rate from the n'* 6-sheIi to the ^shell. The loss of 
energ}' from the fr-shell to the 6-sheII is defined to be positive. 

r,“ are negative for m < n and they are positive for m > n (see Fig. 4.7). Hence a 6-shell 
gains energy from the 6-shells of smaller wave numbers and loses energy' to the 6-sheiIs of 
larger wave numbers. Since is self-similar, the energy lost from a shell (n - to n is 
equal to the energy lost from n to the shell (n -1- An). Thus, the net magnetic energy transfer 
into any inertial range shell is zero, and the energy cascades from the smaller wave numbers 
to the higher wave numbers. 

We now discuss the simulation results for the kinetic energy transfer rates 7““ from the 
n** to the shell. The results are shown in Fig. 4.8 where we have plotted T"" versus m 
for various values of n. We find that the most dominant transfers are from the i/-shell to 
m — n ± 1. From the sign of these transfers w'e see that kinetic energy is gained from (n - 
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Figure 4.8: The energy transfer rate from the w-shell to the M-shell. The loss 
of energy from the 6-shell to the 6-shell is defined to be positive. The boxed points 
represent energy transfer from the u-shell to the u-shell. 

shell and lost to [n 4- 1)*^ shell. This means that the local transfers from the adjacent shells 
result in a forward ceiscade of energy towards the large wave numbers. The transfers from 
other shells are largely negligible except for the transfer to the shell m = 1 (shown boxed on 
the left of the figure), which represents a loss from high wave number modes to the m = 1 
shell. This non-local transfer (to the first shell) produces an inverse cascade to the small 
wave numbers, observed earlier in Section 4.5. Now it is clear that this inverse cascade is 
due to non-local transfers. The local and non-local transfers are seen to possess altogether 
different features: the former is a forward cascade which seems largely self-similar, while the 
latter is an inverse cascade only to the first shell. 

We schematically illustrate in Fig. 4.9 the energy transfer between shells. In this figure 
we show directions of the most significant transfers in the inertial rajige. The arrows indicate 
the directions of the transfers, and thickness of the arrows indicates the approximate relative 




small wave numbers inertial wave numbers 


Figure 4.9: A schematic representation of the direction and the ni.ii'niludM of energy lran.«fer 
between tt-shells and 6 -shells. The relative magnitudes of the differtmt t ran.sfers has been 
represented by the thickness of the arrows. The non-local with the T’ shell have 

been shown by dashed lines. 

magnitudes. Since the local transfer rates are self-similar, the transfers from any other shells 
m the inertial range will also show the same pattern. .An inertial range 6 -siu‘ll gains significant 
amount of energy from the smaller n-shells through both local and rron-local tran.sfers. and 
it also locally gains energy from the smaller 6 -shells {T^n)- "^^e energy gained by a 6 -shell 
from the smaller 6 -shells is exclusively lost to the larger 6 -shells and the energ\' gained from 
the u-shells is mainly lost through homologous transfer [T^) to the corresponding ii-shell k 
small fraction of the energy is also lost to the larger u-shells. In addition to the energy from 
the 6 -shells 5 a u-shell also gains energj' from smaller u-shells by local transfer. The energy of 
u-shells is mainly lost locally to higher 6 -shells and u-shells, but a significant amount is also 
transferred to smaller u-modes by a non-local inverse cascade. 

- -As illustrated in Fig. 4.9, there is a transfer of energy from the first u-shell to the first 
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6-shell. This is the most significant gain of magnetic energy from tjie kinetic energy. Under 
steady state the magnetic energy gained by the first shell gets transferred to the higher 6-shells. 
This transfer yields the forward cascade of magnetic energy discussed in Section 4.5 

To summarise, we find that there are various types of energy transfers in MHD turbu- 
lence: local, nonlocal, homologous, and non-homologous. The main local transfers are the 
forward magnetic energy transfer and the forward kinetic energy transfer, the homologous 
transfer from the magnetic to the kinetic energy (all in the inertial range). There is nonlocal 
transfer from inertial range u-shells to the first 6-shell and the first u-shell. and from the first 
u-shell to the 6-shells. There is a transfer from kinetic energy to magnetic energy in the first 
shell. 

In this section we had extensively described various cascade rates (fluxes) and shell-to- 
shell energy transfer rates. It is clear that the complete picture is quite complex. Infact, some 
of the features we have seen contradicts earlier conjectures and results. These differences are 
discussed in the next section. 

4.7 Discussion 

\\b have investigated the features of kinetic and magnetic energy transfer between scales 
at low values of cross-helicity in a quasi-steady state of forced 2-D MHD turbulence. .A 
summary of the results is given below. For the following discussion refer to Figs. 4.3 and 4.9. 

1) There is a net transfer of energy from the kinetic to the magnetic. 

2) There is an energy transfer to the large-scale magnetic field from the large-scale 
velocity field (w-sphere to 6-shere flux Hj^ in Fig.4.2), and also from the small-scale velocity 
field (flux from 6-sphere to modes outside the u-sphere !!„> in Fig.4.2). The forrher transfer is 
of a greater magnitude than the latter. The latter transfer has also been observed in numerical 
simulations and EDQNM calculations of Ishizawa and Hattori [19, 20]. The magnetic field 
enhancement is primarily caused by a transfer from the «-sphere to the 6-sphere. Indeed the 
first few 6-modes get most of this energy. 




3) Other significant transfers, not noted in earlier work, are from the .li,. vHocitj, 

field to the small-scale magnetic field (w-sphere to outside the /Asjdieici .uif. at, interesting 
reverse transfer from the small-scale magnetic field to the .wmah •>* ale 

modes outside the 6-sphere to modes outside u-sphere). 

4) There is an inverse cascade in the velocity field - this is etcssisfen! %-.ith tin. 
vation of Ishizawa and Hattori arising from numerical siiiiulat .‘dO hih! is also > -inviqent 
with the EDQNM closure calculations [18, 19]. This inver.sc* cascade of kinetic .-nergr is 
dri\en by the reverse transfer of energy from magnetic to the velcH-it}' field at iIk* small 
scales .A.lthough the flux study points to an inverse cascade of kiiH*tie energ)'. the slicll-to- 
shell energy transfer rates reveal the following feature. There exi.sts hrjth an ii!vers(’ and a 
forward transfer of kinetic energy. The inverse transfer is priina:il>' nmjorw/. coming from 
the larger wave numbers to the first few modes. The forward ca.scade i.s local, i.e.. between 
same sized eddies and is self-similar in the inertial range. 

5) There is a forward cascade of magnetic energy toward.s the small scales. This is 
consistent with other recent numerical simulations [20]. EDQNM closure- calculations also 
\ield a magnetic energy transfer to the small-scales [19]. The- magnetic eme-rgy transfer is 


primarily local and is found to be self-similar in the inertial range. 

6) In the inertial range we found a dichotomy between the hornologtms transfers (defined 
as the transfers between kinetic and magnetic energy shells of same wave numbers ) and non- 
homologous transfers (which are the transfers between shells of different wave numbers). 
These transfers are more complex than had been anticipated [18] and are summarised and 
discussed in detail in Section 4.6 (also see Fig. 4.9). From point (1) of this section we recall 
that the magnetic energy enhancement is mainly caused by the energy transfer to the large- 
scale magnetic field from the large-scale velocity field. This observation is important in lieu 
of the interest that exists regarding the physical mechanism responsible for the enhajicement 
of magnetic energy. Pouquet [18] found that the large-scale magnetic energy is destablised by 
the small-scale magnetic energy. In context of 3-D turbulence, Pouquet and Patterson [21] 
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conjectured that energy is redistributed between the small-scale magnetic and kinetic energies 
and that the the enhancement of large-scale magnetic energy is due to an inverse cascade of 
magnetic energy from small-scale to large-scale magnetic field — we, infact find a forward 
cascade of magnetic energy. In a more recent work, Ishizawa and Hattori [19] attributed 
the enhancement of large-scale magnetic energy to energy transfer from small-scale velocity 
field to the large-scale magnetic field. We have also found such a transfer in our simulation. 
However, the strength of this transfer is much smaller as compared to the strength of energy 
transfer from large scale velocity field to the large scale magnetic field (from u-sphere to 
6-sphere). 

The EDQNM closure calculations yield a net transfer of energy to the large-scale mag- 
netic field from the small-scales (magnetic-t- kinetic) [18, 19]. In our simulations, the energy- 
lost by the large-scale magnetic field to the small-scale magnetic field is much larger than the 
energy gained by the large-scale magnetic field from the small-scale velocity field. Hence, 
contrary to the predictions of EDQNM calculations, the large-scale magnetic field loses en- 
ergy to the small-scales. A similar observation was also made by Ishizawa and Hattori in 
their simulations [20]. Thus, it appears that EDQNM calculations do not yield the correct 
strengths of the various transfers. 

The picture obtained by us for the energy fluxes and the shell-to-shell transfers are 
consistent and complement each other. The fluxes and the transfer rates discussed here 
could find applications in the dynamo problem. In astrophysical objects, like galaxies and the 
sun, the magnetic field is thought to have arisen due to the amplification of a seed magnetic 
field. It must be borne in mind that our caclulations are two-dimsensional and devoid of 
magnetic helicity and kinetic helicity. A three-dimensional calculation (with the inclusion 
of magnetic and kinetic helicities) of various fluxes and shell-to-shell energy transfer rates 
will yield important insights on some unresolved issues concerning enhancement of magnetic 
energy. 

In all the numerical simulations of previous chapters, an initial state was time evolved 
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using MHD equations. In the next chapter we will discus.'* the irdiurnt <• uf !hr of jije 

initial state on the time evolution of some global quantiti**!-. 



Chapter 5 

Initial condition sensitivity of global 
quantities in magnetohydrodynamic 
turbulence 

5.1 Introduction 

In earlier chapters we had studied energy spectra and energy transfers using direct numerical 
simulations. Global quantities could affect the energy spectra and transfers. In Chapter 2. 
we had also considered the effect of global quantities like normalised cross-helicity on the 
energy spectra and certain energy fluxes. In this chapter we will study the influence of some 
fine details of the initial condition (to be specified below) on evolution of global quantities. 

In fluid turbulence, the evolution of the velocity u(x) at a given position or of a given 
Fourier component u(k) is known to be sensitive to the details of the initial conditions, e.g.. 
phases of ufk) (to be defined rigorously later). However, the evolution of global quantities 
like total energy are generally presumed to depend only on initial values of total energy and 
energy spectrum. The averaging over many modes appear to wash out the effects of the initial 
phases after a reasonably long time. The total energy in simulations with the same initial 
energy and spectrum, but with the modes chosen randomly, evolve along nearby trajectories: 
this is demonstrated in Fig. 5.1. It has been a common belief that in MHD turbulence also, 
the evolution of global quantities, e.g., total energy and cross helicity, depends only on the 



5.1 Introduction 



initial values of global quantities and their spectra. The ev(»Iufion uf 

earlier been studied by Ting et al. [31], Matthaeus et al. [HH;. H ’..v.p kiu! WVIter 

Pouquet et al. [9] in which they found dynamic alignment and ^anun<. Mt'h.-: 



Figure 5.1: Energy evolution for fluids. The two cases shown hen* differ ciily in the phases 
of the initial conditions. The modes in the runs were |u{k)| with A = fl.O and 0.4. 


In this chapter we show numerically that under certain conditions in 2D MHI) tur!>ulence. 
the evolution of global quantities may not depend solely only on their initial value.s. but may 
depend significantly on more subtle features like the phases of complex Fourier modes (to be 
defined below) of the initially prescribed fields of the dynamical variables. In other words, a 
knowledge of the gross initial features as specified by the global quantities is not sufficient 
under all conditions to determine the evolution of global quantities. We only choose phases as 
a con\enient way of demonstrating the inadequacy of specifying the initial global quantities 
and spectra alone for computing the evolution of certain global quantities. We find that the 
e\olution of cross helicity show's sensitivity to the initial phases in simulations with small 
\alues of initicil cross helicity. 
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The relevant global quadratic quantities for this chapter are 


(-5-1) 

(5.2) 


the magnetic energy, 

E, ■■ 


1 / b^dv 

2 J umt vol. 


= + h (k)f-2i?e(z+(k)-z (k)) 


(5.3 i 


and the fluid energy, 

= 


\L 


u^dv 


unit vol. 


8 


^ z'‘‘(k) ^ + z (k) ^ + 2i?e (z‘''(k) • z (k)) 


(5.4) 


where and z~ are complex conjugates of z"*" and z“ respectively. The total energy is 
E = {E'^ + E~)/2 and the cross helicity is He = {E'^ — E~) 12. 

The plan of the chapter is as follows. In Section 5.2 we define the problem which we 
tackle in this chapter. In Section 5.3 we describe the numerical procedure and the method 
of generating initial conditions for the different runs whose evolution is then compared in 
Section 5.4. The discussion follows in Section 5.5. The contents of this chapter are taken 
from the paper by Dar et al. [96]. 


5.2 Initial condition dependence. 

We denote the complex Fourier modes z'^(k) by |z=‘^(k)| exp(f^k ), where 9^ are the phases of 
the modes. All the three global inviscid invariants E, He, and A are independent of phases, 
while the Alfven ratio depends on the phase difference 9^—9^. Ting et al. [31] found that 
the Alfven ratio affects the evolution of global quantities; it follows from their observations 
that the initial phase difference 9^ — 9^ would affect the global evolution. In this paper 



5 3 Simulation approach and initial conditions 


106 


we demonstrate numerically that even keeping this initial pliii.se c.ffe;-!,. t;xrd, rhangeof 
absolute value of the initial phases 9^ aflfects the global evulution. 

In our simulations we investigate the effects of the initial ph.-.s*'.- <>j. :h>‘ ;b total 

energy E, normalised cross helicity cTc, and Alfven ratio r^. Itie it'inp.c.i; of 

has been the subject of investigation in a number of earlier .stmiies !>>. !i, li;i. .!! , hi several 
of these studies, a, has been observed to increase with time ;8. H. . a hrlwi iour termed a.H 
dynamic alignment. However, Biskamp and Welter [ 8 ] observed in Hieir >:nr.!la!Hins that the 
tendency towards dynamic alignment decreases with the increase in Ih". number, and 
cT(. could even decrease at high enough Reynolds number [ 8 ]. I nig ft nl. j-li j 

few cases of decreasing ctc for small values of initial Oc and EjA. In flie^.e eaiher studies the 
effects of absolute phases had not been studied. 


5.3 Simulation approach and initial conditions 

We solve the 2 -D incompressible MHD equations with hypervisrosit y in Kqs. (2.7) of {'Iiapter 
2 . The numerical algorithm applied to solve these equations ha,s iiotui discns-sed <‘arlier in 
Section 2.3 of Chapter 2 . For the simulations discussed in this ciiapter we chmise »/ = p = 
5 X 10“^ for runs on a grid of size 512 x 512, and we choose u = /u - lO"’’ fur runs on a grid 
of size 256 x 256. The hj'perviscosity related parameter is chosen to be 20 for runs on 
both the grids. The time step dt used for these runs is lO"^. The .simulations are carried up 
to the final time f f,nai = 50. 

The simulations are performed for various initial sets of and r .4 values, 'i’he initial 
conditions are generated by first fixing r^, E. The chosen value of ta cleterminu.s the phase 
difference 9-^ —9^. The initial E and a,, determine jz*(k)|. Note that the absolute phase 
is still a free parameter. Only modes within the annular region 1/2 < (k[ <3/2 are non-zero 
and each of the modes within this region receives equal energy (i.e., |z=*=(k)|^ = /M, where 

M IS the number of modes m the shell). The initial states are generated thus for only half the 
modes - the remaining half are conjugate to them. 
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Run 

N 

Seed 

rA 

A 

(Jcit = 0) 

ac{t = 50) 

(Jc increases/decreases 

mhdl 

512 

50 

1.5 

0. 

0.1 

0.06 

decreases 

rnhdl” 

512 

50 

1.5 

0.4 

0.1 

0.13 

increases 

mhdl"" 

512 

575 

1.5 

0. 

0.1 

0.20 

increases 

mhd2 

512 

50 

5.0 

0. 

0.1 

0.22 

increases 

mhd2’' 

512 

50 

5.0 

0.3 

0.1 

0.05 

decreases 

mhd2”' 

512 

575 

5.0 

0. 

0.1 

0.13 

increases 

mhd3 

256 

50 

5.0 

0. 

0.5 

0.87 

increases 

mhdS” 

256 

50 

5.0 

0.4 

0.5 

0.88 

increases 

mhd4 

64 

50 

1.5 

0. 

0.1 

0.02 

decreases 

mhd4“ 

64 

50 

1.5 

0.6 

0.1 

0.34 

increases 

mhd4” 

64 

50 

1.5 

1.0 

0.1 

0.17 

increases 

mhd5 

64 

50 

2.0 

0. 

0.1 

- 0.02 

decreases 

mhd5“ 

64 

50 

2.0 

0.2 

0.1 

- 0.01 

decreases 

mhd5” 

64 

50 

2.0 

0.4 

0.1 

0.22 

increases 

mhd6 

64 

50 

1.0 

0. 

0.5 

0.79 

increases 

mhd6“ 

64 

50 

1.0 

0.2 

0.5 

0.74 

increases 

mhd6“” 

64 

50 

1.0 

0.4 

0.5 

0.72 

increases 


Table 5.1: Initial values of the random number generator seed A, and for runs performed 
on grid of size N x iV. The initial and the final values (at t final = 50) of ctc are also shown. 
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The phase sensitivity of the evolution of ctc, E, and ta are studied by comparison of pa 
of simulations in which initial are different. We change the initial phases in two ways, 
one case we change uniformly for all the modes by an amount A, while m the other Cc 
the phases are changed by using a different random seed in the random number generate 
The initial global quantities E, He, r^, and their spectra remain unchanged under these phe 
changes. The evolution of cTc for a variety of initial ac, ta and A values are shown in Tal 
5.1 (pairs of simulations are shown together; for example, mhdl differs from mhdl* only 
that its initial z*(k) fields have to be shifted from the latter’s by A = 0.4). 

The N = 512 runs with t final — 50 are very time intensive. Hence only the small init 
(Jc runs, which we found to be sensitive to the phases, were carried out for N=512. A lar 
number of runs were performed on N=64 to explore a wider range of initial conditions. P 
these results showed behaviour consistent with the results discussed below which are bas< 
on the high resolution runs N=256 and 512. 

5.4 Results 

In this section we will discuss the results of the simulations for the initial conditions descril 
in Table 5.1. 

In our simulations, the small cr^ runs showed the most significant dependence on initi 
phase d^. For ae = 0.1 and = 1.5, we choose A to be 0.0 and 0.4 in mhdl and mhd 
respectively. It is seen in Fig. 5.2 that phase shifting has a marked effect on the evolution . 
(Te- For A = 0.0 (mhdl) increases from its initial value of 0.1 to its final value of 0.12! 
whereas for A 0.4 (mhdl ), decreases to a final value of 0.06. The total energy (Fig. 5.f 
and the Alfven ratio (Fig. 5.4) do not appear to be affected much by the phase shift. We als 
compare two simulations (mhdl and mhdl** in Table 5.1.) in which the initial phases ai 
generated using different random seeds. For these cases also the effect on the evolution c 
CTc (Fig. 5.2) is significant, but the corresponding effects on the evolution of the total energ 
(Fig. 5.3) and ta (Fig. 5.4) are not noticeable. 
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Time 

Figure 5.2; .N'ormaliseci rross Felicity (cTc) evolution for initial ac 
shown correspond to nnhdl. mhdl”, and mhdr* in Table 5.1. 
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Figure 5.3; Evolution of total energy (E) for same initial conditions as in Fig. 5.2. See Table 
5.1 for description of mhdl, mhdl*, and mhdl**. 
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Figure 5.4: Evolution of Alfven ratio (r^) for same initial conditions as in Fig. 5.2. Look at 
Table 5.1 for description of mhdl, mhdl*, and mbdl**. 



Figure 5.5: Normalised cross helicity (ctc) evolution for initial = 0.1, = 5.0. The curves 

shown correspond to mhd2, mhd2*, and mhd2** in Table 5.1 





Figtin' ■>.(>: Evolution of total energy [E) for same initial conditions as in Fig. 5.5. Look at 
Table 5.1 for dchcription of mhd2, mhd2*, and mhd2**. 
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Figure 5.7: Evolution of Alfven ratio (r^) for same initial conditions as in Fig. 5.5. Look at 
Table 5.1 for description of mhd2, mhd2*, and mhd2 *. 





re 5.8: Normalised cross helicity {ac) evolution for initial cTc = 0.5 and = 5.0. The 
ss shown correspond to mhd3, mhd3* in Table 5.1. 
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e 5.9: Evolution of total energy [E) for same initial conditions as in Fig. 5.8. Look at 
5.1 for description of mhd3. mhd3*. 
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Figure 5.10: Evolution of Alfven ratio for same initial conditions as in Fig. 5.8. Look at 
Table 5.1 for description of mhd3, mhd3*. 

V\e also .studied the effects of initial phases for the same initial (Tc, but with a large initial 
r..i(5.0). The runs mhd2 and mhd2* show the effects of changing A, while mhd2 and mhd2“" 
show the effect.s of different random number generator seeds. The results obtained for this 
ca.se are similar to the run for initial condition with r 4 = 1.5. In Fig. 5.5 it is seen that for 
initial value of A = 0.0 (mhd2), increases and for initial A = 0.3, decreases. The effect 
of changing A on the total energy (Fig. 5.6) and (Fig. 5.7) is seen to be small. Similar 
results are obtained if we change the seed of the random number generator (compare mhd2 
and mhd2'“). Hence, (Fig. 5.5) is sensitive to the change in the initial phases whereas 
the total energy (Fig. 5.6) and (Fig. 5.7) are not sensitive. Earlier, Ting ti al. [31] had 
observed a decrease in cTc for small initial (Tc. 

We also perform runs at higher initial values of Uc (mhd3 and mhd3 in Table 5.1). The 
effect of phase shifting for initial <t, = 0.5 and r 4 = 5.0 is shown in Figs. 5.8, 5.9, 5.10. It is 
seen in Fig. 5.8 that the changes in evolution caused by phase shifting are relatively smaller 
for high initial Cc values as compared to small initial Cc discussed above. From Fig. 5.9 it 
can be seen that the effect of A on total energy is also small and (Fig. 5.10) remains 
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:nsitive to the change in A. We have performed more runs than have been shown here 
in all cases ( 7 ^ was seen to increase for high Uc values. This increase in cFc is consistent 
ti simulations performed earlier for high values [8, 9, 29-31]. 

3 Discussion 

in the numerical results presented here we conclude that phases of the initial modes pla}' an 
'ortant role in the evolution of Uc, atleast for cases with small initial Gq values. For higher 
les of ( 7 c, phases do not appear to affect the evolution oi Gc by any significant amount. In 
he runs the total energv and were seen not to have any significant dependence on the 
ses. 

The origin of the phase effects discussed here is not clear at this moment. We need to 
mine the evolution more carefully before reaching any definite conclusion. These studies 
Id find applications in understanding the solar wind observations in which Gc has been 
Jrved to decrease [14, 15, 97, 98]. 

It has been demonstrated in the paper that the evolution of normalised cross helicity 
ignificantly affected by subtle features of the initial condition especially at low initial 
s helicities. This observation will require us to be more circumspect in drawing conclu- 
s based on arbitrary initial conditions and to exercise more care in choosing the initial 
litions in MHD turbulence. 



Chapter 6 
Conclusions 


This thesis dealt primarily with two aspects of turbulence : 1) The scaling behaviour of energ\' 
spect ra, and 2) the features of energy transfer between the velocity scales, the magnetic scales, 
and between tlie velocity and the magnetic scales. We briefly summarize our findings below 
and tiieii we discuss a few^ future research directions. 

6.1 Summary of results 

We compared the predictions of the Kraichnan-Iroshnikov-Dobrowolony’s phenomenology 
and the Ktdmogorov-like phenomenology in our simulations. It was difficult to distinguish 
between the scaling e.xponents —3/2 and —5/3 from the plots of the energy spectra obtained 
from our simulations. So, following the approach suggested earlier by Verma et al. [Ih 
we compared the fluxes n"*" and 11“ since the KID phenomenology and Koimogorov-like phe- 
nomenology make different predictions for the ratio of ll'^'/n . In our simulations we averaged 
our results over a steady state obtained by forcing small wavenumber modes, therefore giving 
better accuracy compared to Verma et aPs simulations. We also compared the results of 
Grappin’s extension [4] of KID phenomenology with our simulations. We found that the KID 
phenomenology aad its extension by Grappin are inconsistent with the simulation results, 
but the predictions of Kolmogorov-like phenomenology are in reasonably close agreement. 
Hence, our results favour the Kolmogorov-like phenomenology where the energy spectra is 
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proportional to k In addition to our own results, the Kolmogorov-like phenomenology is 
also supported by the Solar Wind data [14—16], higher Reynolds number simulations in 2-D 
and 3-D [12, 13], and a recent RG calculation [26]. We have studied the dependence of the 
Kolmogorov constants with (7^ (or f E~) and found that there is an approximate linear 
relationship between C~ /C'^ and [E'^ / E~). The constant are independent of <7^ over a 
wide range, and have a dependence on it only at low values (for C"^) and at high values for 
(for C-). 

An important component of turbulence phenomenologies and its description is the energy 
transfer between w’avenumbers. The MHD and Navier-Stokes equations show that energy is 
transferred, in spectral space, to a mode k from modes p and q such that the three wavenum- 
bers satisfy the condition k -t- p -|- q = 0. The method for computing the combined energy- 
transfer into mode k from the other two modes is well known. We devised an approach for 
describing the ‘mode-to-mode’ energy transfer between two modes of this triad, mediated by 
the third mode. We found that the mode-to-mode energy transfer is expressible as a sum of an 
‘effective transfer’ and a ‘circulating transfer’. These two have the property that the effective 
transfer contributes to the modal energy change, whereas the circulating transfer does not 
alter the modal energy because the energy transferred from k to p gets transferred from p to 
q and then back to k from q. By imposing the requirements of rotational invariance, galilean 
in\-ariance, finiteness, and retaining only certain linear terms in the series we showed that the 
circulating transfer is equal to zero. However, in absence of a complete proof, the value of the 
circulating transfer remains uncertain. Still, we have drawn important conclusions in terms 
of effective energy' transfer which is independent of the circulating transfer. We introduced 
the idea of effective mode-to-mode transfer for energy' transfer between a (1) pair of velocity 
modes, (2) a pair of magnetic modes, and (3) between a velocity and a magnetic mode in 
the triad. Using effective mode-to-mode transfer we defined energy fluxes and shell-to-shell 
transfers which cannot be defined in the older approach where only the combined transfer is 
known. The energy fluxes and shell-to-shell energy transfer rates were studied in a numerical . 
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sin.., la, Urn .,f 2-D MHD turbulence and important conclusions were drawn. 

Un s.,.-li,.,l the energr- duxes and shell-to-shell energy transfer rates in a numerical 
sniiulanuii of .’.D MHD turbulence where white-noise velocity forcing is included at low 


wa\'«’nt!!uiM*r.s. The results of 
ohiained hy over a 


energy cascade rates and the shell-to-shell transfer rates were 
quasi-steady state. Several interesting observations were made 


in tmr sinmiatiDiis. 

l-;«)in tile ^iinr.latiuiis we drew conclusions about the direction of energy transfer between 
dilferent scales, anti their relative magnitudes. The details of the energy transfers are ex- 
piaineit in the tontlusion of Chapter 4 and are pictorially represented in Figs. 4.3 and 4.9 
of tlie sanu’ (. hapter. 1 he most important outcome of our study can be stated in the follow- 
ing pt)iuts ; 1 ! i Enhancement of magnetic energy occurs due to the energy transfer to the 
large'seal(‘ magnetic field from the large-scale velocity field, (2) there is a forward cascade of 
magnt'tic eiierg.v from the large-scale magnetic field to the small-scale magnetic field. Our res- 
ult.s ai'ff important since they clarify some of the proposals made earlier regarding the physical 
mecluinism for tlie generation of large-scale magnetic field. Pouquet [18] had found that the 
larg<'--'Culf m.tunctic energy is destablised by the small-scale magnetic energy. Ishizawa and 
Hattori [Hf] [>njpusf>(l that the large-scale magnetic energy is enhanced due to energy transfer 
from the .small-.scaie velocity field to the large-scale magnetic field. Although we find such 
a transfer, the dominant transfer to the large-scale magnetic energy is from the large-scale 
kinetic energy. Pouepn't. and Patterson [21] proposed in context of 3-D turbulence that there 
an inverse cascade of energy from small-scale to large-scale magnetic field. However, our 
numerical calculations give a forward cascade in 2-D turbulence. A simulation should be 
done to verify whether this feature is also present in 3-D turbulence. 

Interestingly, we found an inverse cascade of kinetic energy to large-scale velocity field 
from small-scale velocity field, which was driven by the transfer of energy to the small- 
scale velocity field from the large-scale magnetic field. The inverse cascade was a result of 
non~local transfer to the large wavenumber velocity modes. Interestingly, the local transfer 
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had a forward cascade from the large-scale velocity field to the small-scale velocity field. By 
studying the energy transfer to a wavenumber-shell from all other wavenumber shells, we found 
that the ^homologous' (energy transfer between shells at same wavenumbers) transfers w^ere 
responsible for the reverse transfer from magnetic to velocity field. The non-homologous 
(energy transfer between shells at different wavenumbers) transfers were from the velocity to 
the magnetic field (see Chapter 4). 

Numerical simulations of MHD turbulence are performed by time evolving an initial 
state. The initial condition is generated by fixing the global quantities, and the phases are 
generated randomly. It has been assumed till now that the random phases will have no 
influence on the evolution of global quantities. We found that evolution of <Jc can depend 
significantally on the random set of initial Fourier phases, atleast when ac Is low. Hence w^e 
show' that one should be careful in drawing conclusions based on simulations with arbitrar\ 
initial conditions and care must be exercised in choosing the initial condition. 

6.2 Future research directions 

In this thesis we have not considered the effect of a mean-magnetic field on the energy spectra 
and the fluxes. The turbulence phenomenologies assume isotrop}'. However, Bq is expected to 
introduce anisotropy in the energy spectra. A mean magnetic field is known to affect energy 
cascades in turbulence [84. 85]. However, the effect of Bo on the various energy cascades 
and shell-to-shell transfer rates was not investigated in this thesis. Since, most astrophysical 
plasmas are immersed in an external magnetic field it is important to study the effect of such 
a field on the energy transfers. 

The dependence of the Kolomogorov constant only on ac was investigated. However, 

are also expected to depend on In future simulations, the effect of should be 
studied. 

The chief source of magnetic energy enhancement in 2-D is large-scale kinetic energy- 
transfer to large-scale magnetic energy (H^^). Our conjecture is that in 3-D too one probably 
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will of><.Tve snrh l)ehaviour. Considering the implications of the results to dynamo theories, 
we hdiev. that similar analysis for 3-D simulations will yield many important insights. We 
u)u!'< .lu. lilt th» abtne analysis because of lack of computational power . We hope that the 
probh in h« takfn up in the future. Our study has been performed over a quasi steady 
std.e |>« liod with A low .Aifxen ratio. In order to understand the build-up of magnetic energy 
starting truni a .setul value it is important to perform a similar study at high Alfven ratio. We 
ha\ ( *iu,a "on.e p.' !iiuiiiar\ decaying simulations for large Alfven ratio in which the magnetic 
fifh* jia .< .is« s. then saturates and finally decays. The growth of magnetic energy in the earlv 
phase i.s paiii* nlarh interesting. A careful study in this direction is being carried out and 
will he rrp«>rt»‘d at a later date. 

In some of iIk’ popular models, like the o-dynamo [81], the mean magnetic field gets 
lUnpIi.hei! in presence of helical fluctuations. The 2-D turbulence simulations discussed here 
doe.s not have any helicity. The numerical verification of a-dynamo and similar models 
recpiires 3-1) and ^uinewhat complex field and geometry configurations. Still, we clearly 
see an incr«Nase in magnetic energy in the early phases. This mimics the magnetic energy 
enhanceiinuit in galaxies and the universe. Hence, our study presented here could possibly 
find application in the galactic dynamo. 

Thv fluxes also find important applications in various phenomenological studies. For 
exainplt', Wrma (t al. [99] estimated the turbulent dissipation rates in the Solar Wind and 
obtained the it>iu!>cratur<' \ariation of the Solar Wind as a function of distance. The various 


cascade rat«*s disc 


nssed here could be useful for various astrophysical studies. For example. 


ni‘^ and Ill'y can be used for studying the variation of of the Solar Wind. 

Recent numerical and analytical results appear to show that Kolmogorov-like phenomen- 
ology is applicable to MHD turbulence. Our results support this belief. Still higher resolution 
■2-D and 3-D numerical simulations and more exhaustive analytic calculations are necessarv 
for drawing definite conclusions. 

We have shown that very rich physics is present in energy transfer. We have investigated 
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some of its aspects. Full 3-D simulations will yield many valuable insights which will help us 
understand many complex processes, e.g.. MHD dynamo, local vs. nonlocal transfers, etc. 

To conclude, we have discovered some important results in MHD turbulence. Yet. many 
more questions remain to be answered, and some new questions have been posed. 



Appendix A 

Derivation of the circulating transfer 
between the velocity modes in a triad 

la St'ction 3.4 of Chapter 3 we showed that the mode-to-mode transfer is given by = 
+ Xs- interpreted AA as a circulating transfer and [Eq. (3.12)] as an effective 
tiKxk’-to-itiode transfer. In this appendix we shall determine the circulating transfer Xj\ for 
fluirl and .MHD turbulence, using symmetry considerations. In the first section we shall 
<iisru.s.s .Y^ in fluid turbulence and in the next section we will discuss the MHD turbulence. 

A.l Xa in fluid turbulence 

Tlu‘ energy transfer between pair of modes should be invariant under a rotation (because it 
is a .scalar), and galilean transformation, and it should be finite for all values of k, p, q, 
u(k), u(p), u(q). We will determine the form of the circulating transfer by imposing these 
and a few other symmetry constraints. 

We construct a scalar expression for .Y^ dependent on k, p, q, u(k), u(p), u(q). and 
which is cubic in u and linear in k, p or q, to satisfy the for dimensional reasons. We write 
the general expression for such a scalar as 

Aa = Rei [k.u(q)] {[u(k).u(p)]a + [u(k).u(q)]/3 + [u(p).u(q)]7 + 

[u(k).u(k)]^ + [u(p).u(p)]i/ + [u(q).u(q)]/i} + 

[k.u(p)] {[u(k).u(p)]i7 + [u(k).u(q)]C + [u(p).u(q)]c + 
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[u(k).u(k)]« + [u(p).u(p)]A + [u(q).u(qM + 

[p.u(k)] {[u(k).u(p)]V’ + ••• + [u(p)-u(q)]tt^ + + [u(q)-u(q)]x} + 

[p.u(q)] {[u(k).u(p)]r + ... + [u(p).u(q)]^zi> + — + [u(q).u(q)]yj} + 
[q.u(k)] {[u(k).u(p)]<T + ... + [u(p).u(q)]e + - + [u(q).u(q)]s'} + 
[q.u(p)] [[u(k).u(p)];r + [u(k).u(q)]/5 + ... + [u(q).u(q)H) (.A.l) 

where all the coefficients a, (3^ etc. are atmost non-dimensional functions of the wavenumbers 
k, p, q and of u(k), u(p), u(q). in the expression in Eq. (A.l) indicates the real 
part of the expression. We should have infact written Xa ss a linear combination of the 
real and the imaginary parts of the entire expression in the above equation. However, we 
will only consider the real part explicitly. The imaginary part can be treated in the same 
manner and the following conclusions, although explicitly stated for the real part will also 
be valid for the imaginary part. The most general expression should also contain terms of 
the kind i?e[k.u(q)]i?e[u(k).u(p)], i?e[k.u(q)]/m[u(k).u(p)], and /m[k.u(q)]/m[u(k).u(p)]. 
However, we will explain below that such terms should not be present in Xa on account of 
galilean invariance. Hence, from the outset we do not include such terms. The ta^k now is 
to determine the coefficients in Eq. (A.l). 

The interactions involves all three modes of the triad, i.e., if the Fourier coefficients of 
an}- of the wavenumbers is zero, interaction between the modes is turned off. The mode-to- 
mode transfer should respect this fundamental feature of the triad interactions. Therefore, 
if u(k), u(p) or u(q) approach zero, Xa should approach zero as well. We have just 
argued above that q,/ 3, etc., in Eq. (A.l) are independent of the magnitude of u(k), u(p), 
u(q). By inspecting Eq. (A.l) for Xai we find that a few terms in the expression are 
independent of the magnitude of atleast one of the Fourier components. For example, the 
term [k.u(q)][u(k).u(q)]/? is independent of u(p). Such terms should be dropped from the 
expression of circulating transfer. After dropping such terms, we obtain 

Aa = Re{ [k.u(q)][u(k).u(p)]a + [k.u(p)][[u(k).u(q)]C + 
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[p.u(k)][[u(p).u(q)]u; + [p.u(q)][[u(k).u(p)]T + 
[q.u(k)][[u(p).u(q)]^ + [q.u(p)][[u(k).u(q)]p) (A. 2) 

I he recjui lenient of finiteness of the mode-to-mode transfer imposes restrictions on the 

form of A'^. hi Eq. (A.l) the dimensional terms are [k.u(q)][[u(k).u(p)], [k.u(p)][[u(k).u(q)]. 

etc-, and th(>y are finite for all values of k, p, q, u(k), u(p), u(q). The coefficients a, 3. 

etc. are ■diiiuuisionless scalars’. Therefore they can be written as 

^•u(p) u(k).u(p) u(k).u(p) 

kp kq' ’^'|u(p)|’ A:|u(q)|’‘"’ |u(k)||u(p)|' lu(k)||u(q)|””^ 

where tlie arguments of the function have been non-dimensionalised. Consider triads in 

wiiieh one of the wavenumbers, say q, tends to zero. Then the terms like k.p/Aig with q in 

the denominator diverge. Similarly if one of the Fourier coefficients, say u(q). is equal to zero. 

thi'ii tlie argnnient.s and etc. diverge. Hence terms depend on the magnitude 

♦ 

of vector.s must be dropped. Only the arguments which depend on the angles between vectors, 
e.g.. . i " . should be included. Therefore, all the dimensionless coefficients should 

be a funrtion of only these arguments, i.e., q = /[^, ..., K§nS|- 

We now write the coefficient a, C, etc. in the form 

a = (A.4) 


wh<*re is a constant, and 


kp kq pq 


k.u(p) (a, k.u(q) (2) p-u(k) (2) p.u(q) J2) q-'i(k) 

Of' ' = q) ", TT + Q 2 7Tr7rTT'h% , ;rTT + “4 m ^“5 1 /i_m “6 


(A.5) 

( 2 ) q-u(p) 


A:lu(p)l 

v(3) 


A:lu(q)| 


p|u(k)l Piu(q)| 


q|u(k)l 




, 3 ) u(k).u(p) , 131 u(k).u(q) ( 3 ) u(p).u(q) 

|u(k)||u(p)| “ lu(k)||u(q)| = Kp)|Mq)l 


<?iu(p)r 

(A.6) 

(A.7) 


The quantities ^(3) ^re ‘linear’ in the scalars of the type k.p/fcp. k.u(p)/A:iu(p)|. 

and u(k).u(p)/lu(k)|lu(p)| respectively. The quantity contains terms like (k.p) /k p , 
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[k.u(p)]^/Ar^([u(p).u(p)], etc. which are higher order in these scalars. All the other coeffi- 
cients etc. in Eq. (A.l) can be similarly divided into linear and non-linear terms. 

We now study the consequences of galilean invariance of Aa on Eq. (A.l). Consider 
two frames of reference S and S' with S moving with respect to S' with velocity Uq. We 
denote the Fourier coefficient of velocity in frame S and S' by u(k) and u'(k) respectively. 
The relationship between u(k) and u'(k) is given by 

u'(k) = u(k) exp‘^*‘ +Uo^(k). (A. 8) 

The second term is non-zero only for the k = 0 modes. In the subsequent discussion relating 
to galilean invariance we can drop this term for sake of simplicity without affecting the 
derivation. The transformation rule for scalars of the type k.p, k.u(p), u(k).u(p) between 
the two frames of references S and S' can be deduced from the transformation property of 
u(k). 

k'.u'(q') = k.u(q) exp’*'* (A. 9) 

u'(k').u'(p') = u(k).u(p) exp'^^'+P^’^®' (.A.. 10) 

It follows from the above two transformation rules that k.u(q) and u(k).u(p) are not invariant 
under a change of reference frame. However, a product of these quantities can be seen to be 
invariant 

[k'.u'(q')][u'(k').u'(p')] = [k.u(q)][u(k).u(p)]exp‘*‘'+P+‘*)-’^°' 

= [k.u(q)][u(k).u(p)] (A.ll) 

using k-f p-f-q=0. However, 

i?e[k'.u'(q')]i2e[u'(k0.u'(p')] / i?e[k.u(q)]/Ee[u(k).u(p)], (A.12) 


i?e[k'.u'(q')]/m[u'(k').u'(p')] ^ i?e[k.u(q)]/m[u(k).u(p)], 


(A.13) 
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aiui 

lm[k .u (q )|/m[u'(k').u'(p')] ^ /m[k.u(q)]/m[u(k).u(p)] (A. 14) 

I lius. wt> shown that the terms of the kind given in the above three equations are not 
galiiean invariant. 

( ItMil}. the iiuariance under Galilean transformation demands that all of the vectors 
u(k), u(p), u(q) appear either once, twice, or n times in the scalar. A scalar of the type 
ik.ul q )j ;iii k j.uiqtl having two or more repeated wavenumber indices is not invariant since 
the term appearing in the exponential is not zero for such a scalar. The transformation 
properties of any quantity formed from the scalars like k.p, k.u(p), u(k).u(p) can be 
ol>t;tiued in a similar manner. If the quantity is not galilean invariant then it should not be 
inciufled in the expression for mode-to-mode transfer. 

We ha\'e shown above that scalars of the type k.u(q) and u(k).u(p) are not Galilean 
invariant [Ecjs. (.A. 9) and (A. 10)]. It follows that etc., etc. are not 

Galilean invariant. To preserve the Galilean invariance of Aa, the coefficients and 
should be drt>p[>ed from Eq. (A.4). However, the nonlinear quantity can be constructed 
in a Galilean invariant form. Currently, we have dropped these higher order terms. Our 
derivation of A'a therefore has the limitation of being strictly valid only upto the linear order. 
.A complete proof containing all order of terms in Xa is beyond the scope of this thesis and 
i.s left for future work. After dropping the terms the Eq. (A.4) reduces to 

a = (A.15) 

Using equations (A.4) and (A.15) we now write X in Eq. (A.2) as, 

A'a = Re{ [k.u(q)][u(k).u(p)]{ao + + "3^} + 

[k-u(p)][u(k).u(q)]{Co + Ci^ + + Cs^} + 

[p.u(k)][u(p).u(q)]{a;o + 
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k.p , k.q 


p.q. 


[p.u(q)][u(k).u(p)]{ro + + 

k.p k.q p.q. 

[q.u(k)][u(p).u(q)]{6 + 

k.p k.q 

[q.u(p)][u(k).u(q)]{po + ^ 


(A.16) 


W’e have dropped the superscript '1' from ctj \ c ^2 incompressibility rela- 

tionship q.u(q) = 0 we obtain p.u(q) = — (k -f q)-u(q) = — k.u(q). Therefore, fourth term 
in the above equation can be absorbed into the first term. Similarly, using the relationships 
q.u(p) = -k.u(p) and q.u(k) = -p.u(k), the last term can be combined with the second 
term and the fifth term combined with the third term. Now, we can write 

A A = Re{ [k.u(q)][u(k).u(p)]{ao -h + 02 -^ + ^3“} + 

[k.u(p)][u(k).u(q)]{Co + Ci^ + ( 2 ^ + Cs^} + 

[p.u(k)][u(p).u(q)]{aJo + u;i^^ + a;2^ (A. 17) 


by making the replacements (a — r) — >■ a, (C “ /?) “t ( and (w — u. 

The quantity represents the energy circulating in the anti-clockwise direction, i.e.. 
A'a flows along p — > k -4 q -4 p [see Fig. 3.4 in Section 3.4.2]. We will now denote the 
circulating transfer along p -> k by A’(klplq), transfer along k -4 q by A''(q|kip) and so 
on for the circulating transfer between each pair of modes. The above A''a is the circulating 
transfer from p to k, i.e., A''(k|p|q|) = A’a- Then by symmetry the circulating transfer from 
q to k can be written as 

-'^'(k|q|p) = Re{ [k.u(p)][u(k).u(q)]{ao + -f- + Q 3 — } + 

kq kp pq 

[k.u(q)][u(k).u(p)]{Co + Ci^ + + Cs— } + 

fc? kp pq 

[q.u(k)][u(p).u(q)]{a;o -f 012— 4-013— }), (A. 18 ) 

kq kp pq 
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and the cirrulating transfer from p to k is 

A'lpikiq! =ft( [p.u{q)I[u(fc).u(p)l{a„ + a.ia + „,P:S 

kp pq kq^ 

[p.u(k)|[u(p).u(q)l{C„ + + 

kp pq kq 

[k.u(p)][u(k).u(q)]{t<;o + uji-~ + u; 2 ^-^ + (A. 19) 

kp pq kq 

I'luin f* ig. d. l in Chapter 3. A(kjp|q) flows along the anti-clockwise direction while A^klqlp) 
and A (p;kjqi flows along the clockwise direction. The former should therefore be equal and 
opp<Jsit(* To the latter two transfers, i.e., 

^(klpiq) + -^(kiq|p) = 0. (A.20) 

and 

^"(k|p|q) + A:(p|k|q) = 0. (A.21) 

Substitution of expre.ssions A'(k|p|q), .A(k|q|p), X(p|kjq) [Eqs. (A.17)-(A.19)] into the first 
of these two relationships yields 

[k.u(q)][u{kl.u(p)]{(ao + Co) + {<^i + + (03 + 

[k.u(p)](u(k).u(q)]{(Co + oo) + (Ct + ^2)^ + (C2 + + (Ca + 

[p.u(k)][u(p).u(q)]{(wi - ^ 2 )^ + - ‘*'’ 1 )^} = 0- 

(A.22) 

and 

[k.u(q)][u(k).u(p)]{(a 2 - + (as - “ 2 )^} 

-|-[k.u(p)][u(k).u(q)]{(Co d-t^o) + (Cl + (C2 + ^3)-^+ (C 3 + <^2)^} 

-h[p.u{k)][u(p).u{q)]{(u;o + Co) + (‘^i + Cx)-^ + (^2 + C3)-^+ (‘^3 + C2)^} =0. ; 

(A.23) ! 
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Each of the three terms in the two equations above has a different functional dependence on 
u(k), u(p )5 u(q). If we only rotate u(k) about k without changing its magnitude then each 
of the terms will change by a different factor depending on the direction of u(p), u(q) which 
can be independently chosen. Since the sum of three terms should be zero and each tfn*m 
can change a different factor, each of the terms should be individually zero. In addition, this 
should hold for all wavenumber triads. Therefore each term inside the curly bracket should 
also be individually zero. This condition gives the the following relationships between the 
constants: qq + Co = + C 2 = 0?<^2 + Ci = 0 , 03 + Cs = 0,a;i — 0^2 == 0 from Eq. (A. 22), 

and wo'o + Co = 0 , 02 ““ <^3 ~ Oi Ci "h ~ C 2 “h ^3 — Cs + ^^2 ~ 0 from Eq. (A. 23). 

We now consider the following geometry - the wavenumber triad form an isosceles triangle 
with jp| = |q|, and the Fourier modes u(p) and u(q) are perpendicular to the plane of k, 
p, q and are equal in magnitude. For this geometry the clockwise and the anti-clockwise 
direction are equivalent. There is no preferred direction along which the circulating transfer 
ma\’ flow. Therefore, for the geometry shown the circulating transfer should be zero. We 
may now write the circulating transfer [Eq. (A. 17)] for this geometry as 

-Y =Re( [p.u(k)][u(p).u(q)){(i^.+L;,)|B + ^3Ef! + a}) 

Kp pq 

= 0 (A.24) 

With u(p) and u(q) fixed perpendicular to the plane of the triad, X should be zero for 
all wa\'enumber triads which form an isosceles triangle. For .Y to be zero for ever\' such 
geometry, each term inside the curly bracket should be equated to zero giving the additional 
relationships: u;i = — u; 2 ,t <;3 = 0,a;o = 0. 

We have obtained 17 relationships between the constants a, etc. These equations may 
be sohed to get the values of all these constants. It can be verified that the only solution of 
these equations is 0 '’s= 0, C’s= 0, tu’s= 0, a = 6 = c = 0. 

Therefore we get JY = 0. It follows that the mode-to-mode transfer 

We have thus determined the mode-to-mode transfer. The proof used rotational in- 
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v^>ria.,r.', galilaan invariance, finiteness of mode-to-mode transfer, and symmetry of X with 
rt s|K < . k. p. q. uf k),u(p), u(q). However, in the series we have taken only linear terms 
in k.p. k.U(, p I. u( k).u(p]. A complete proof containing all orders of terms is beyond the 
scop.' of this thesis. It may be possible to obtain a completely general proof by finding 
.some 11 '.-nr tries which we have ignored in our proof. Otherwise a novel approach may be 
rr* L 

A. 2 A'a in MHD turbulence 

The circniating transfer between the velocity modes in MHD can also be shown to be equal 
to itero to linear order terms in k.p.k.u(p),u(k).u(p) . The arguments leading to this are 
nearly the same a.s that for the fluid case. However, we additionally need to show that the 
general rxprrs.-vion for the circulating transfer Xa for MHD should not contain b(k), b(p), 
b(q) aiul can also be given by Eq. (A.l) 

In the previous section we showed that each term in Xa should depend on u(k), u(p), 
u(q) ”■ thi.H followed from the fundamentally three-mode feature of the interactions. The 
expression for the combined energy transfer to u(k) from u(p) and u(q) is identical for the 
fluid ami for the xMHD case. Hence, in the case of MHD also, each term in the circulating 
transfer between the velocity modes should be dependent on u(k), u(p), and u(q). Thus, 
the general (xxpression for the circulating transfer between the velocity modes in MHD is the 
same as that for Fluid turbulence. The rest of the arguments leading to AIa = 0 in MHD are 
identical to those in Section A.l. 



Appendix B 

Derivation of circulating transfer 
between the magnetic modes in a triad 

In Section 3.7.2 we had shown that the mode-to-mode transfer from b(p) to b(k) can be 
written as where given by Eq. (3.49) is the effective transfer between 

the modes and Va is the circulating transfer. In this section we will derive the expression for 
the circulating transfer Va- The derivation for I a is similar to the derivation of A'a for fluid 
turbulence given in Appendix A. 

Since the mode-to-mode transfer is a scalar, the circulating transfer Fa will also be a 
scalar with a dependence on the wavenumbers k, p, q and the Fourier coefficients u(k), 
u(p), u(q), b(k), b(p), b(q). It should be linear in k, p or q and cubic in the Fourier 
coefficients u(k), u(p), u(q), b(k), b(p), b(q) for dimensional reaisons. 

The most general form of Fa that satisfies these properties can be written as 

Fa = Fi') + Ff ) + Fi^) + Fi") (B.l) 


where 

>A ’ = [k.u(q)] {[u(k).b(k)]... -t- [u(k).b(p)]... + [u(k).b(q)]... + 

[u(p).b(k)]... -b [u(p).b(p)]... + [u(p).b(q)]... -b 
{[u(q).b(k)]... -b [u(q).b(p)]... + [u(q).b(q)]...} -b 
[k.u(p)] {[u(k).b(k)]... + .... + [u(k).b(q)]... + .... 4- [u(q).b(k)]... -b .... + [u(q).b(( 
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Ip.u. k .] I iufk ,.b(k)]... + .... + [u(p).b(q)]... + .... + [u(q).b(p)]... + .... + [u(q).b(q)]...} 

(B 


= ik.utq): {;b(k).b{k)] + [b(k).b(p)]a + [b(k).b(q)]...+ 
ibip!.b(p)] + [b(p).b(q)]... + [b(q).b(q)].'..} + 

[k.lM P)1 {ib(k ).b(k)]... + .... + [b(k).b(q)]/3 + .... + [b(q).b(q)]...} + 

[p.ui k ): {?b(k).b(k)]... + .... + [b(p).b(q )]7 + .... + [b(q).b(q)]...}) 

{B.3) 


R(i [k.b(q(] {;U(k).b(k)]... + [u(k).b{p)](5+[u(k).b(q)]... + 

'u(p).bfk)]7? + [u{p).b(p)]... + [u(p).b(q)]... + 

{!u(q).b(k)]... + [u{q).b(p)]... + [u(q).b(q)]...} + 

[k.b(p)] {fu{k).b(k)]... + .... + [u(k).b(q)];/ + .... + [u(q).b(k)]i/ + .... + [u(q).b(q)]...} ^ 
|p.b(k)] {{u(k).b(k)]... + .... + [u(p).b(q)]^^ + .... + [u(q).b(p)]x + [u(q).b(q)]...}) 

fB.- 

= R€{ [k.b{q)] {(u(k).u(k)]... + [u(k).u(p)]^ + [u(k).u(q)]... + 

[u(p).u(p)]... + [u(p).u(q)]... + [u(q).u(q)]...} + 

{[u(q).b(k)]... + [u(q).b(p)]... + [u(q).b(q)]...} + 

[k.b(p)| {[u(k).u(k)]... + [u(k).u(q)]<?i + .... + [u(q).u(q)]...} + . 

[p.b{k)] {[u(k).u(k)]... + [u(p).u(q)]e + .... + [u(q).u(q)]...]) 

(B.5) 
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Following Appendix A we consider explicitly only the real part, which is denoted by Re in 
the above equations. If we take a linear combination of the real and the imaginary parts 
instead, the derivation and the conclusions regarding Va change. In addition, terms 

of the i?e[k.u(q)]/2e[u(k).u(p)], Re[k.n{q)]Re[u(k).h{p)l etc. are not included as they 
violate galilean invariance - we had explained this in Appendix A in context of A a- The 
three dots ‘...’ in the above equations in front of some of the terms indicate that a coefficient 
multiplies each of these terms but which we have not shown explicitly as they do not enter 
most of the future discussion. The coefficients a, /5, etc. in the above expressions for 

are non-dimensional functions of the wavenumbers k, p, q and of u(k), u(p), u(q), 
b(k, b(p), b(q). To determine Va we should now determine these coefficients. Note that 
besides the terms of the kind [k.u(q)]{...}, [k.u(p)]{...}, [p.u(k)]{...} in Eqs. (B.2)-{ B.o), 
there are also terms of the kind [p-u(q)]{...}, [q.u(p)]{...}, [q.u(k)]{...}. However, using the 
incompressibility constraints k.u(k) = 0 and the constraint k -f p -h q = 0, we have absorbed 
the latter three terms into the former ones. Now we will impose a few restrictions on the 
form of 1 a. 

The interactions between the magnetic modes fundamentally involve three-modes — the 
combined energy transfer to a magnetic mode, say b(k), depends on the Fourier coefficients 
at all three wavenumbers (k, p, q). For example, the combined energy transfer to b(k) from 
b(p) and b(q), given by 5'’^(k|p,q) = -ile(f[k.u(q)][b(k).b(p)] -|- f[k.u(p)][b(k).b(q)]), 
depends on u(p), u(q), b(p), b(q), and b(k). The mode-to-mode transfer should not 
break the fundamentally three-mode feature of these interactions, i.e., the mode-to-mode 
transfer should have Fourier coefficients from all wavenumbers k, p. q. Hence we will drop 
all those terms from ^ to ^ which are independent of atleast one of the wavenumbers 
k, p,q. Then, we get 

^ [k-u(q)] {[u{k).b(p)]... + [u(p).b(k)]...} -t-[k.u(p)] {[u(k).b(q)]... -|- [u(q).b(k)]...} 

[p.u(k)] {[u(p).b(q)]... -f- [u(q).b(p)]...}) (B. 
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= Rn [k.u(q))(b(k).b(p)]a + [k.u(p)][b(k).b(q)]/3 + [p.u(k)][b(p).b(q)h) 

(B.7) 

= i?t( [k.b(q)] {[u(k).b(p)]5 + [u(p).b(k)], 7 } + 

[k.b(p)] {[u(k).b(q)]/i + [u(q).b(k)]i/} + 

[p.b{k)] {[u(p).b(q)]^’ + [u{q).b(p)]x}) (B-S) 

= Rti [k.b(q)j[uik).u{p)]^ + [k.b(p)][u(k).u(q)](^ + [p.b(k)][u(k).u(k.;... ) 

(B.9) 

j)iif the' following additional restriction on the form of I'a. If none of the magnetic 
nK)des in the triad are gaining/losing energy to the other two modes, i.e., if the combined 
energy transfer to each one of them is equal to zero, then even the circulating transfer should 
be etjnal to We observe the following additional feature of the interactions. The com- 

bined energ\' transfer to every magnetic mode in the triad from the other two magnetic modes 
is equal to zero, if any two of b(k),b(p),b(q) are zero (see Eq. (3.36 in Chapter 2). The 
morle-t 4 ,-im>ile transfer should also respect this feature of the interactions, i.e.. the circulating 
transfer should be zero if (bCk). b(p)=0}. or {b(k), b(q)=0}, or {b(p), b(q)=0}. Now. 
let us put b(k)=0, b(p)=0 in fI” to Ky>. It can easily shown that Fa will be zero, for 
arbit rary values u's and b's. only if the all coefScients in > and Ff are themselves equal 
to zero. Hence we drop F]” and F^’ from Fa [Eq. (B.l)]. We get. 

Fa = Ff’ + Ff 

and 

Yl^^ = Re{ [k.u(q)][b(k).b(p)]a+ [k.u(p)][b(k).b(q)]^ -f- 

[p.u(k)][b(p).b(q)]7 ) 

rf ) = Rei [k.b(q)] {[u(k).b(p)]<5 + [u(p).b(k)]7?} + 


(B.ll) 
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[k.b(p)] {[u(k).b(q)]/z + [u(q).b(k)H + 
[p.b(k)] {[u(p).b(q )]0 + [u(q).b(p)].\}) 


(B.12) 


In Appendix A we had considered the derivation of circulating transfei between 
the velocity modes. By imposing the condition of finiteness of the circulating transfer we 
had concluded that the coefficients should be independent of the magnitude of the wavenum- 
bers and the Fourier coefficients u(k), etc. On imposing the constraint that should 
be finite and following the reasoning in .Appendix A, we can conclude that the coefficients 
a, 3. etc. should be independent of the magnitude of the wavenumbers and of the Fourier 
coefficients of the velocity and the magnetic fields. Therefore, we may write the coefficient 

— /TLE J^-b(p) u(lt)-u(p) b(k).b(p) u(k)-b(p) -i rp, othf>r rnpffiripnfc; ran 

" - Jikp^ •••’ i-|u(p)| • A:|6(p)| ’ •••’ |u(c-)llu(p)p |6(A-)l|6(p)|’ |«(i:)||6(p)p ^ coemcicnts Can 

also be written in the same manner. 

We will now write the coefficients a, 3, etc. as 

a = + 0^“*^ (B.13) 


where is a constant, and 


o(i)_ , ai)P:3 


(B.14) 


a(2) = 


^,( 2 ) Wp) , ^.(2) k.u(q) ^(2) p.u(k) (2)P-u(q) (2) q-u(k) (2)q.u(p) 

■ A:Kp)r“^ *^|u(q)r“^ pKk)|+“* p|u(q)|+“' 5 Wk)i 


(2) k.b(p) , _(2) q-b(p) 

" kib(p) A ib(p)i’ 


(B.15) 


a 


(3) 


= „(3) , ^(3) u(k).u(q) , ( 3 ) u(p).u(q) 


u(k)||u(p)| |u(k)||u(q)| """ |u(p)|iu(q)| 

b(k).b(p) (3) b(k).b(q) . (3) b(p).b(q) 


Q 


(3) 


|b(k)|lb(p)r"® |b(k)||b(q)| 


+ Oi 


+ 


+ 


(3) u(k).b(k) ^ (3) u(k).b(p) 

' |u(k)||b(k)r“» Kk)||b(p)I + - + “'=Kq)||b(p)|- 


“ fb(p)l|b(q)| 

(3) u(q)-b(p) 


(B.16) 


where etc. are all unknown constants. The quantities qO) are ‘linear’ 

in the scalars of the type k.p/fcp, k.u(p)//:lu(p)[, k.b(p)/Jb|b(p)|, u(k).u(p)/|u(k)l|u(p)l. 
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b(k).b(p1/ib,k)iib(p)!,andu(k).b(p)/|u(k)l|b(p)|. The quantity a(") contains higher order 
terms in these scalars like {k.p)^ ky, [k.u(p)]VP([u(p).u(p)], etc. All the other coefficients 
et( . in F.qs. (B.11)-(B.12) can be similarly divided into linear and non-linear terms. 

I ht' m()fic-tt>-nu)de transfer should be invariant under a galilean transformation. In 
Apj>en<iix A, we had given the transformation rules for the scalars of the type k.p. k.u(p). 
u(k).u(p) and for a combination of these like [k.u(p)][u(k).u(q)]. For the present case we 
also nee*(i to con.sicier the transformation properties of the magnetic field. MHD equations 
hold for low velocities and only quantities upto the order of vjc are included in the equations. 
Tiuler this approximation the magnetic fields are equal in the two frames S and S' with a 
relative velocif}' Uo- Then, the Fourier coefficient b(k) in the frame S and b'(k) in the frame 
S' are relattni as 

b'(k) = b(k)exp’(‘'-^°)‘ (B.17) 

This differs from the galilean transformation of u(k) by the term Uo5(k) only. Using the 
above relationship between between b(k) and b'(k) the transformation rules for k.b(q). 
b(k).b(p). u(k).b(p) are 


k'.b'(q') = k.b(q)exp‘(‘i-^°)‘, 

(B.IS) 

b'{k').b'(p') = b(k).b(p)exp‘(>'+P^-U‘>‘, 

(B.19) 

u'(k').b'{p') = u(k).b(p)exp’(‘'+P)-^°‘ . 

(B.20) 


These scalars are not galilean invariant (in Appendix A, we had likewise shown that the 
scalars of the type k.u(p), u(k).u(p) are not galilean invariant). It follows that 
etc., etc. are not invariant under a galilean transformation. To preserve the galilean 

invariance of Va, the coefficients qP>. o^, and other such coefficients should be 

dropped from Eq. (B.13). The nonlinear quantity a*''*, /S*"*', etc. can be constructed in a 
galilean invariant form. Even then we will drop this term. We had also dropped this term 
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from the expression for Xa in Appendix A. Our derivation for Ya. like the derivation of 
in Appendix A, is is thus valid only upto the linear orders. After dropping the coefficients 
q(2)_ 0(3)^ 0(4) other coefficients of this type, Eq. (B.13) reduces to 

0 = Q(°) + a('l (B.21) 

Using Eqs. (B.21) and (B.14) we can now write and Y^^ in Eqs. (B.ll )-(B.r2) as 
Yi^'> = Re{ [k.u{q)][b(k).b(p)]{Qfl°l ^ + 02^ + 03^} + 

[k.u(p)][b(k).b(q)]{/ 3 (°l + ^ + .?3^} + 

[p.u(k)][b(p).b(q)]{7l°l + qi^ + 72 ^ + 73^} ) (B.22) 

= Re{ [k.b(q)][u(k).b(p)]{( 5 <°l q. + S2^ + <^3^}+ 
[k.b(q)][u(p).b(k)]{r;(°) + 71^ + ^2%? + ^ 3 ^}+ 
[k.b(p)][u(k).b(q)]{/il°l + + P2^ + /^3^}+ 

[k.b(p)][u(q).b(k)]{i/l°l + + ^^3^}+ 

[p.b(k)][u(p).b(q)]{^l°l + 4>i^ + <f>2^ + <?i>3^}+ 
[p.b(k)][u{q).b(p)]{x<°) + Xi^ + X2^ + X3^} ) (B. 23 ) 

where the superscript ‘T has been dropped from etc. 

The quantity Ya = Y^^ + represents the energy circulating in the anti-clockwise 
direction, i.e. Ya flows along p -4 k -)• q -4 p [see Fig. 3.10 in Section 3.7.2]. We will now 
denote the circulating transfer along p -4 k by yA(k|p|q), transfer along k -4 q by >'A(q|k|p) 
and so on for the circulating transfer between each pair of modes. We take V'a = Ya^ + >1^’ 
from Eqs. (B.22)-(B.23) to be the circulating transfer from p to k, i.e., >A(k|plq|) = I'a In 
the same manner yA(k|qlp)) = ll^^(k|qlp) -|- yi^^(k|qlp) can be written as 

>'?^(klq|p) = Re( [k.u(p)][b(k).b(q)]{a(°) + + «3^}+ 

[k.u(q)][b(k).b(p)]{^(o) f 3 i^ + ^2^ + /?3^}+ 
[q.u(k)][b(p).b(q)]{7(°)-h7i^-h72^ + ^3^} ) (B.24) 
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>‘r’(}<:qip) = Re{ [k.b(p)][u(k).b(q)]{^(0) + <Ji^+<52^ + <j 3 M}+ 
[k.b(p)][u(q).b(k)]{77^°) + + ‘f] 2 ^ + 73 ^}+ 

[k.b(q)][u(k);b(p)]{^(°) + /fi^ + ^ 2 ^ + /^ 3 ^}+ 
[k.b(q)][u(p).b(k)]{z/(°) + 4- 1 / 2 ^ + 

[q.b(k)][u(q).b(p)]{^(°) + 4 >i^ + <i> 2 ^ + 03 ^}+ 
[q.b(k)][u(p).b(q)]{x(°) + XiT? + X2^ + \3^} ) (B.25) 

and the circulating transfer from p to k, yA(p|k|q)) = y^^^(p|klq) + >'r*(p|k|q) can be 
vvritt(Ui as 

>f '(Pikiq) = i?e( [p.u(q)][b(k).b(p)]{a(o)+ai^ + a 2 ^ + a35^}+ 
[p.u(k)][b(p).b(q)]{/?(°) + +/^ 2 ^ +/^ 3 ^}+ 

[k.u(p)][b(k).b(q)]{ 7 (“) + 71 ^ + 72^ + 73%?} ) (B.26) 

>f (p|k|q) = Re{ [p.b(q)][u(p).b(k)]{5(°) + 5,^ + 82 ^ + <^ 3 %?}+ ^ 

[p.b(q)][u(k).b(p)]{r/(°) + + 72 %? + 73%?}+ 

[p.b(k)][u(p).b(q)]{p^°’ + + P2%? + P3%?}+ 

[p.b(k)][u(q).b(p)]{r/(°) + + ^^ 2 %? + ^^3%?}+ 

[k.b(p)][u(k).b(q)]{(?i><°> + 4>i^+ <?^2%? + <?!>3%?}+ 
[k.b(p)][u(q).b(k)]{x*°) + Xi^ + X 2 %? + \3%?} ) (B.27) 

yA(k|piq) flows along the anti-clockwise direction while VA(klqip) and yA(p|k|q) flow along > 
the clockwise direction. The former should therefore be equal and opposite to the latter two : 

transfers, i.e., i 

TA(k|p|q) + yA(k|q|p) = 0, (B.28) - 


I 

i 
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and 


l"i(k|p|q) + Ki(p|k|q) = 0- 


(B.29) 


Replacing the expression for hA(k|p|q) from Eqs. (B.22)-(B.23) and for l^lklqlp) from 
Eqs. (B.24)-(B.25) into Eq. (B.28) we get 
[k.u(q)][b(k).b(p)]{(a(°) + + (ai + ^ 2 )^ + (a^ + + (^3 + + 

[k.u(p)][b(k).b(q)]{(/ 3 '°^ + 0 ^( 0 )) + (A + ^^ 2 )^ + (A + + A + °‘3^}+ 

[p.u(k)][b(p).b(q)]{( 7 i - 72 )^ + (72 - 7 i)^+ 

[k.b(q)][u(k).b(p)]{((J*°^ + + (^1 + 1^2)-^ + (<^2 + ^i)'^ + (A + '“ 3 )^)+ 

[k.b(q)][u(p).b(k)]{( 77 ^°^ + + ^ 2 )^ + (72 + + (73 + 

[k.b(p)][u(k).b(q)|{(/u(°^ + <5^°^) + + (^2 + ^ 1 )^ + (/“3 + A)— }+ 

Kp fcq pq 

[k.b(p)][u(q).b(k)]{(i/^°^ + 7]^°^) + {ui + 772 )^ + i^2 + 7 i)^ + (^^3 + 73 ) — }+ 

kp kq pq 

Ip.b(k)][u(p).b(q)]{(<^(°^ - x'°^) + (A - Xa)^ + (A - Xi)^ + (A - Xs)— }+ 

kp kq pq 

lp.b(k)][u(q).b(p)){(x"’' - + (xi - + (X 2 - A)|3 + (X3 - *)^} = 0 

(B.30) 


and replacing the expression for yA(k|p|q) and yA(p|k|q) from Eqs. (B.26)-(B.27) into 
Eq.(B.29) we get 

[k-u(q)][b(k).b(p)]{(o 2 - 03 )^ + (q3 - ^ 2 )— }+ 

kq pq 

[k.u(p)][b(k).b(q)]{(/ 3 (°) + 7 ^°^) + (A + 72 )— + (A + 73 ^ + A + 72 —}+ 

Kp kq pq 

[p.u(k)][b(p).b(q)]{( 7 (°^ + + (71 + A)^ + (72 + A)^ + (73 + A)— + 

kp kq pq 

[k.b(q)][u(k).b(p)]{(^(°) - 7 ^°^) + (A - 7 i)^ + (A - 73 )^ + (A - 72 )—}+ 

kp kq pq 

lk.b(q)][u(p).b(k)l{(,<") - i(”)) + -h)^ + (m - ■5^)^}+ 
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ik.bfpIHul k).b(q)]{(p«l + ^O) + + 4 ,^)^ + + (f,3 + *) 

kq pq 

[k.b(p];!u(q).b(k)j{(^m + + (^, + ^,)|p ^ (3,^ ^ + (3,3 ;^3)M}+ 

kq pq 

[p.blkfjmi plbiq)]{{d(0) + ^ ^ ^ ^ 2 )—}+ 

^P kq pq 

:p.b.:kj][uiq).btp)l{(x'°* + 2 /'°^) + {xi + 2 / 0 ^^ + (X 2 + 2 / 3 )^ + (X 3 + 2 ^ 2 )—} = 0 

fep kq pq 

{B.31) 

Siiu'<‘ each term in the above equation has a different functional dependence on u(k), u(p), 
u(q), b(k), b(p), b(q). each term should be individually zero. It follows that the constants 
in the above equations should satisfy the following relationships : from Eq. (B.30) we get 
= 0.01 + J2 = 0,02 + A = 0)0:3 + A = 0,7i — 72 = 05^°^ + = 0,^i +P2 = 

0. <>2 -t Ml = 0, Sz + M 3 = 0. 7^°* + 2/^°^ = 0, 7i + 2/2 = 0, 72 + 2^1 = 0, 73 + 1/3 = 0, = 

0. Oi — \2 = O.O 2 - \ i = 0, 03 — X 3 = 0, and from Eq. (B.31) we get 02 — <^3 = 0,,5*°^ + 7 f°^ = 

0. -^1 + 7i = 0. .i 2 + 73 = 0, A + 72 = 0,5^°^ — 7 ^°* = 0,(5i — 7 i = OjA — 73 = 0, a — 72 = 

+ 0 ^*^** = O.Ml +01 = O.M2+^3 = O,M3 + 02 = 0,2/^°^+X^°^ = 0, 2/2 + Y 3 = O,M3 + 02 = 0. 
Tlu* above relationships can be arranged into the following groups: = 7 *°*}. 

= 0 ^°^ = = X^°^}) {0:1 = -A = 73}, {0:2 = -A = 7 i = 72 = 

03 = -A}. {A = -M2 = 03 = X3 = -2^2 = 7i}) {A = -Ml = 4 >i = X2 = -i"3 = 73}- 

{A = — M:i = 02 = \i = —2^1 = 72}- 

Let us now consider a triad with |p| = |q|, k.p = k.q and with the Fourier modes 

u(p).u(q),b(p),b(q) perpendicular to the plane formed by this triad. We argued in .V 

pendi.x .A. that the circulating transfer should vanish for such a geometry. For this geometry 
we can write 1'^ from Eq. (B.22) and ( B.23) as 

Vk = R 4 [P-u(k)][b(p).b(q)]{7(“) + (71 + 72)^ + 73^}+ 
[p.b(k)][u(p).b(q)]{0(°^ + (01 + 02)^ + 03^}+ 
[p.b(k)j[u(q).b(p)]{x^°^ + (xi + + Xs^} ) 
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= 0 (B.32) 

The above equality should hold for all triads with |p| = |q| 3-nd k.p=k.q. Therefore each 
term inside the curly brackets should be individually zero. This leads us to the following 
additional relationships between the various constants in the above equation: = O.qi + 

q2 = 0,73 = = 0,<^2 + d>3 = 0,d>i = = 0,xi + X3 = 0,\2 = 0. Similarly. b\ 

choosing |k| = |p|, and u(k), u(p), b(k), b{p) perpendicular to the triad plane we get the 
relationships = 0, ai = 0, cx2 + cxz = 0i = 0, <^1 = 0, <^2 + ^3 = 0^ = O- = 0. 1)2 + 

rjz = 0. and by choosing |k| = |q|, and u(k).u(q),b(k),b(q) perpendicular to the triad plane 
we get = 0, /32 = 0, 13 i + I 3 z = 0, = 0. p2 = 0, pi + pa = 0, 1/^°^ = 0, 1^2 = 0, z/i + 1/3 = 0 

Using these relationships the above groups take the following values: { 0^°^ = — = 

= 0 }, = -pW = <^(°) = 77 (°) = = 0 }, {ai = -/?2 = 73 = 0 }. 

{02 = — 3 i = 7i = 72 = 03 = — /33 = 0}, {<^1 = — /i2 = <f >3 = X 3 — ~^2 = Vl = 0}< 
{^2 = -fil = (f>i = X2 = —^3 = V3 = —^3 = IJ'3 = —<i>2 = “Xi = 1^1 = —772}- Hcnce. some of 
the constants are equal to zero, and all the non-zero constants are related to each other b\- 
just one undetermined constant which we wdll denote by 'C', i.e., {62 = —fJ-i = ... = C) 

Putting the values of these constants in the expressions for and in Eqs. (B.221- 
(B.23) we get 

= 0 (B.33) 

’ = Re{ [k.b(q)][u(k).b(p)]{C^ - C^} + [k.b(q)][u(p).b(k)]{-C^ + C^}+ 
[k.b(p)][u(k).b(q)]{-C|f + C^} + [k.b(p)][u(q).b(k)]{C^ - C^}+ 
[p.b(k)][u(p).b(q)]{C^ - + [P-b(k)][u(q).b(p)]{-C^ + C^} ) 

(B.34) 

with a single undetermined constant. 

For the determination of in Appendix A, the above considerations were sufficient 
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to (It’tormine Aj^. However, in the present case we are still left with one undetermined 
constant. \\e will now impose the constraint of locality which we state as follows: In the 
limit k. q :x: >> p, Ia should tend to zero. If Za remains finite then the noise, which is 
alway s iiH'vitahly present in any flow at the smallest scales, would interact directly with the 
lame s<'ales. Moreover, the locality constraint is also necessary in order that the circulating 
transfer is consistent with the effective transfer between b(k) and b(q) which goes to zero 
for in the above limit. Let us take u(k), u(q) to be perpendicular to the plane of the triad. 
Then we can write >a as 

Vy = ft, lk.b(p)lMk).b(q)l{-C^+C^) + lk.b(p)][u(q).b(k)l{C!^-C“} I 

(B.3-31 

which may be written in terms of the angles of the triangle formed by k, p, q [Fig. B] as 

k.b(p)|[u(k).b(q)]{-C'A:cos(t/)) - Cfccos(^)}+ 
[^k.b{p)][u(q).b(k)]{C^:cos(t(!') + Ckcos{&)} ) 

(B.36) 

We now take the limit k. q oo. keeping i' and p constant. In this limit 0. Making 
u.se of the relationsliips 0 = tt - 1 /. - p, fc/ sin 0 = p/ sin and taking the limit -> 0. we get 
for Va [k.b(p)](u(k).b(q)]{-Cpsin^(t/.)} + [k.b(p)][u(q).b(k)]{Cpsin^('0)}. However, 
in thi.s limit V'a should tend to zero in order to satisfy locality. From this constraint we get 

C = 0. Putting the value of C in Eq. (B.34) gives vl ^ = 0. 

Therefore we get Vk = 0. It follows that the mode-to-mode transfer Iji^ ■ 

In the above proof we used rotational invariance, galilean invariance, finiteness of mode- 

tenmode transfer, the symmetry of Fa with respect to k, p. q, u(k), u(p). u(q), b(k). 
b(p), b(q), and the assumption that the circulating transfer should vanish if two of the 
wavenumbers approach infinity. This proof is valid only upto linear orders, as explained 
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Figure B.l: In the wavenumber triad k, q — > oo >> p the circulating transfer, h ^ = 0 to avoid 
non-local transfer from very large wavenumbers to very small wavenumbers in comparison. 


earlier in this appendix. A completely general proof without taking recourse to the linearity 
assumption is lacking at present. 



Appendix C 


Derivation of the circulating transfer 
between velocity modes and magnetic 
modes in a triad 


III Section 3.7.3 we had shown that the energy transfer from u(p) to b(k) in the triad (k, 
p, qj ran be e.xpressed as '5l*’“(k|p(q) = ^^“(k|p|q) + Za. 5 ?^“(k|p|q) is given by Eq. (3.76) 
of Section 3.7.3 wa,s interpreted as the effective transfer from u(p) to b(k) and Za was 
iut<'ri)r('fed as a circulating transfer between the velocity and magnetic modes belonging to 
a triad. U'e will follow the approach of Appendix A and Appendix B to determine Za in 
this .Appendix. We will impose the following conditions on Za- It should not violate the 
fumla mentally triad nature of the interactions, it should be finite, rotationally invariant, 
galih'an invariant, and should vanish if fc/p -)• oo. 

The circulating transfer Za is a scalar with a dependence on atmost k, p, q, u(k), u(p), 
u(q), b(k), b(p), b(q). We know from Appendix B that such a function can be expressed 
in the form given bv Eqs. (B.1)“(B.5) of Appendix B. The interactions between the velocity 
and the magnetic modes are fundamentally three-mode, i.e., each term in the expiession 
for combined energy transfer to a velocity or a magnetic mode depends on all the three 
wavenumbers (k, p, q). the combined energy transfer to any of the velocity or the magnetic 
modes in the triad. Za should not violate this fundamental nature of the interactions. Hence 
we should drop all those terms from Eqs. (B.1)-(B.5) which are independent of atleast one of 
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the modes and then we get Eqs. (B.6)-(B.9) for Za- The combined energy transfer to velocity 
and magnetic modes in the triad are zero if {b(k), b(p) = 0}, {b(k)^ b(q) — 0}, or {b(p)j 
b(q) = 0}. The terms which are non- vanishing if one of the above set of Fourit^r coefficients 
are zero will be dropped from Eqs. (B.6)-(B.9), giving Eqs. (B.10)-(B.12) In Appendix B 
also we had imposed identical constraints on Eqs. (B.1)-(B.5) to get Eqs. (B.10)-(B.12), The 
dimensionless coefficients in these equations are k, p, q, n(k), ii(p)? b(k), b(p), 

b(q) dependent and we will now attempt to determine these. 

In -Appendix B we had considered the consequences of finiteness and galilcan invari- 
ance on Eqs. (B.10)-(B.12) of Appendix B. We will briefly summarise the conclusions now. 
To ensure the finiteness and galilean invariance of Eqs. (B.10)-(B.T2) the dimensionless 
coefficients in these equations should be independent of the magnitudes of the wax'enum- 
bers and the Fourier coefficients. Then each coefficient in the equations can be written as 
a — where is a constant, defined in Eqs. (B.14)- 

(B.16) of Appendix B. are linear in terms of the kind k.p, k.u(p), u(k).u(p) respectiveh. 
and is nonlinear in these terms, and 0*^1 are not galilean invariant and would 

therefore lead to the non-invariance of Eqs. (B.14)-(B.16) of Appendix B. Hence. and 
should not be included in a. In Appendix B we had also dropped Therefore, we get 
Q = q(°) -f- Following Appendix B we will now write [given by Eqs. (B.10)-(B.r2) in 
.Appendix B] 

= (C.l) 

where 

^ 2 A(k|p!q) = Re{ [k.u(q)][b(k).b(p)]{a(°) + 4^^ -b -b 4^^^}+ 

[k.u(p)][b(k).b(q)]{/?W + 4'^^}+ 

[p.u(k)][b(p).b(q)]{7(0) -b + 71^)^ -b ) (C.2) 

^3A(k|p|q) = Re{ [k.b(q)][u(k).b(p)]{(J(°) + -b 4^^}+ 
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[k.b(q)][u(p).b(k)]{^(o) + + 41)^}+ 

[k.b(p)][u(k).b(q)]{^(0) + +^a)^}+ 

[k.b(p)][u(q).b(k)]{i/(°) + + 4*’^}+ 

[p.b(k)][u(p).b(q)]{^(0) + + 41)^ + 4')^}+ 

[p.b(k)][u(q).b(p)]{x(°) + + X?^} ) (C.3) 

1!h* fjtiantitics Z^‘ = Z|^ + Z|^ represent the energy circulating from u(p) -> b(k) -> u(q! 
— * btp!-.. — > u(p) [see Fig. 3.12 in Section 3.7.3]. We will denote the circulating transfer from 
u^p l -4 bi kj by Z^^^lklplq), from b(p) — >■ u(k) by Z^*’(k|p|q) and so on for the circulating 
transfer between other modes. We take Z^“ = Z2A +*“ Zg^ in Eqs. (C.2)-( C.3) to be 
the circulating transfer from u(p) -4 b(k), i.e., Z^“(k|p|q) = Za- Then by symmetry the 
circulating transfer from u(q) to b(k) can be written as 

Zlr(k|q|p) = Z^Kklqjp) + Z3^1(klq|p) (C.4) 

ZlKkjqlp) = Rei {k.u(p)][b(k).b(q)]{a(°) + + 02?^ + ^3^}+ 

[k.u(q)][b(k).b(p)]{/?W + + ^3^}+ 

[q-u(k)][b(p).b(q)]{7^°) + 7i^ + 72^ + 73^} ) (C-5) 

Zjl{k|q|p) = Re{ [k.b(p)][u(k).b(q)]{5(°) + + <^2^ + <^3^}+ 

[k.b(p)][u(q).b(k)]{77(°^ +771^ + r/2^ + %^}+ 
[k.b(q)][u(q).b(k)]{^(°^ + ^ + ^2^ + /23^}+ 

[k.b(q)][u(k).b(q)]{i/(°) + 1 ^ 1 ^ + ^"2^ + J^3^}+ 
[q.b(k)][u(q).b(p)]{4°^ + + 4>2^ + <^^}+ 

[q.b(k)][u(p).b(q)]{x^°^ + Xi^ + X2|^ + X3^} ) (C-6) 

We have dropped the superscript 1’ in the above equations. Z^“(k|p|q) and ZX^p|k|q) are 
physically the same transfers, but the former represents the transfer from u(p) to b(k) and 
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the latter represents the same transfer from b(k) to u(p) The former should therefore be 
equzd and opposite in sign to the latter, i.e., 

4“ (k|p|q) + Z2‘(p|k|q) = 0 (C.7) 

Z^^(p!k|q) is thus given by the negative of Eqs. (C.1)-(C.3). Using the expression for 
‘^^l^lPi'klq) we can write -^A*’(k|p|q) to be 

Z^''(k|p|q) = Z^Kklplq) + ^3t(k|p|q) (C.S) 

(k|p|q) = Re{ -[k.u(p)][b(k).b(q)]{7(°) + 71^ + 72^ + 73^} 
+[k.u(q)][b(k).b(p)]{o(°) + 

-[p.u(k)][b(p).b(q)]{^(°) + /?i^ + /?2^ + ) 

(C.9) 

■^ 3 i\(k|p|q) = Re{ — [k.b(p)][u(k).b(q)]{(;^(°) + ^^i^ + 4>2^ + 4>3^} 
-[k.b(p)][u(q).b{k)]{x(°^ + Xi^ + X 2 ^ + X3^} 
+[k.b(q)][u(p).b(k)]{^(°) + <^i^ + <^ 2 ^ + <^ 3 ^} 
+[k.b(q)][u(k).b(p)]{77(°l + 7i^ + + ^ 3 ^} 

-[p.b(k)][u(p).b(q)]{;^(0) 

-[p.b(k)][u(q).b(p)]{l/(0) + + ^2^ + 1 ^ 3 ^} ) 

(C.IO) 

and Z^-Nk|q|p) is given by 

(k|q|p) = ^2i(k|q|p) + ^3t(k|q|p) (C.ll) 

^2A(k|q|p) = Re{ -[k.u(q)][b(k).b(p)]{7(0) + 71^ + 72^ + 73^^} 
+[k.u(p)][b(k).b(q)]{a(°) + + 02^ + ^3^} 

-[q.u(k)][b(p).b(q)]{/ 3 ( 0 ) + /?i^ + + /?3^} ) 
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(C.12) 

^'(Itklqlp) = Re{ -[k.b{q)][u(k).b(p)]{,^(o) + + 02^ + <^3^} 

-[k.b(q)][u(p).b(k)]{x(°) + Xi%? + X2^ + X3^} 

+[k.b(p)][u(q).b(k)]{<J(°) + + <^2^ + <^3^} 

+ [k.b{p)][u(k).b(q)]{77W + + + ^3^} 

-[q.b(k)][u(q).b(p)]{^(0) + 

-[q.b(k)][u(p).b(q)]{:/(o) + 1 ^ 1 ^ + 1 ^ 2 ^ + 1 ^ 3 ^} ) 

(C.131 

The sum of the circulating transfer to b(k) from u(p) and from u(q) should be zero. 

i.e.. 

Zi“(k|p|q) + Zi“(k|q|p) = 0 (C.Ut 

and similarl}' the sum of the circulating transfer to u(k) from b(p) and from b(q) should 
aisi> b<‘ e({ual to zero, i.e., 

ZX^kjpjq) + ZX^(k|q|p) = 0 (C.lo) 

Replacing -?i“(k|p!q) and (kjqjp) from Eqs. (C.1)-(C.3) and Eqs. (C.4)-(C.6) into 
Eq. ((M4) we get the following relationships between the constants: a,/3, etc. 0^°) + = 

0, cq + J-i = 0, a2 + /?i = 0, 03 + /?3 = 0, 7i - 72 = 0, = 0, + ^2 = 0, ^2 + Pi = 

0 , 63 + p3 = 0, ?7<°) + = 0, 771 + z/2 = 0, 772 + i/i = 0, 7/3 + 1/3 = 0, </)(°) - = 0, d>i - X2 = 

0, 02 - \ i = 0, d>3 - \3 = 0. After replacing ZX*(k|p|q) and Z^^(klq|p) from Eqs. (C.1)-(C.3) 
and Eqs. (C.11)-{C.13) into Eq. (C.15) we get the following relationships: 02 - 03 = + 

7'°’ = 0, ii + 7i = 0, /32 + 73 = 0, /33 + 72 = 0, - 7^°’ = 0, - 7i = 0, 62-^3 = 0, <^3 - »72 = 

0, ^(0) + p(0) = 0, Pi + (Pi = 0, P2 + <^ = 0, P3 + <^2 = 0, + X^°^ = 0, 1/2 + X3 = 0, P3 + <?^ = 0. 

Let us now consider the following geometry: the three wavenumbers in the triad form an 
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isosceles triangle, and we choose Founer coefficients at the two wavenumbers of equal lengths 
to be perpendicular to k, p, q- We had considered this kind of geometry in Appendix A and 
Appendix B. We had argued that because of the symmetry between the two wavenumbers 
of equal length, the circulating transfer should vanish for this geometry. I his constraint 
had given us an additional set of relationships between the constants in Y^. should 

also be equal to zero for this geometry. Since and have the same functional form, 
given by Eq. (B.22)-(B.23) of Appendix B and Eq. (C.1)-(C.3) of this appendix respectively, 
hence imposition of this constraint on would give the same set of relationships, i.e., 
.,,(0) = 0.71 +72 = 0>73 = = 0 , 4>2 + <h = = 0,\1 + \3 = 0, \2 = 0.a'°' = 

0. oi = 0. a2 + 03 = 0 , = 0, (Ji = 0 , 52 + 53 = 0, = 0 , 771 = 0, m + m = 0 , = 0, = 

0..5i + 33 = = 0,/r2 = 0 ,fii +;23 = = 0,r/2 = 0,i/i + 1^3 = 0. 

The above relationships between the constants in Eqs. (C.2)-(C.3) are identical to the 
relationships between the constants in Eqs. (B.22)-(B.23) of Appendix B. We know from 
.Appendix B that these relationships can be expressed more clearly as follows: { = 

^ ^(0) ^ 0}, = 0), {qi = -32 = b 3 = 0}, 

{q2 = = 72 = as = -03 = 0}, {5i = -fl2 = (1)3 = u = -7^2 = m = 0}- 

{<^2 = -/Wi =4>i = X 2 = = V 3 = -S 3 = /X 3 = - 4^2 = -X 1 = ^^1 = -%}• There is just one 

undetermined constant which we will denote by ‘C’, i.e., {^2 = = C’}- Expressing 

[Eqs. (C.2)-(C.3)] in terms of this single constant, we get 

Z2a=0 (C.16) 

= Re( [k.b(q)][u(k).b(p)j{C^ - + [k.b(q)][u(p).b(k)]{-C^ + C^}+ 

[k.b(p)][u(k).b(q)]{-C^ + C^} + [k.b(p)][u(q).b(k)]{C^ + -CM}+ 
[p.b(k)][u{p).b(q)J{C^ - C^} + [p.b(k)][u(q).b(p)H-C^ + C^} ) 

(C.1-) 


which is the same as the expression for for 1 a in Eqs. (B.33)-(B.34) of Appendix B. 
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If A'/p ‘T-nd qjp -4 oo, while p is kept constant the circulating transfer Za should tend to 
zero. This constraint was also applied on the circulating transfer between magnetic modes. 
>A in Appendix B. We had shown that the constant C = 0 for this constraint to be satisfied. 
Putting C — 0 in Eqs. (C.16)-(C.17) we get Z^" = equal to zero. 

Tht‘ circulating transfer between the velocity and the magnetic modes in the triad is 
ecpinl to zero and the mode-to-mode transfer "^^“(klplq) = ^^“(k|p|q). 

The derivation of Z^“. however took the coefficients to be linear. A completely general 
deri\'ation for Z^” will require imposition of additional symmetry and invariance or some 
other const raint.s. which we have ignored here. 


Appendix D 


Equivalence of Kraichnan’s formula and 
our formula for the Kinetic energy flux 


In this appendix we will show that the flux formula (3.30) based on the effective mode-to- 
mode Transfer is equivalent to Eq. (3.28) derived by Kraichnan. 

Kraichnan[9o] showed that the kinetic energy flux in terms of the combined energy trans- 
fer 5“'* given by 

n;5(A') = - E El'5“(k|p,q) (D.i) 

|k|<K P ^ 

Our formula for flux in terms of the effective mode-to-mode transfer is given by 

ng(A-) = - £ E r"(k|p|q). (D.2) 

|k|<K |p|>K 

The equation ( D.I) can be written in terms of ^’s as follows: 

n“>(A")= E 5I^(^““(k|p|q) + 5^““(k|q|p)) (D.3) 

|k|<K p ^ 

which can be further split up as 

n: 5 (A') = 

E E 5(r”(k|p|q)+^“(k|q|p)) 

|k|<K |p|,|q|<A' ^ 

+ E E 5(r*(k|p|q) + r”(k|q|p)) 

|k|<K |p|,|q|>A- ^ 
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+ E E ;(r“(k|p|q)+r"(k|q|p)) 

|k|<K |p!<K-,|q|>A' ^ 
lk|<K |pl>A:,iq|<A- ^ 

Tlu'.sr transfers are shown in Fig. D.l. 



Figure D.i: The various triads involved in the terms in Eq. (D.4). Triad of type A 'does not 
contribute to energy flux. Triad of type C and D are equivalent and contribute the same 
amount. 


Making use of symmetries between k,p,q in the summations we get the following rela- 
tionships. Since p, q are interchangeable we get 


y: s r“(k|p|q)= E E ir“(k|q|p), 

|k|<K ipI,lq|>A.' |k|<K |p|,|q|>A' 


since k, p are interchangeable, 


X: E r“(k|piq) = E E r“(p|k|q), 

iki<K |pl<AMq|>A' |k|<K |p|<AMql>A' 


since k, q are interchangeable, 


^ Y, r‘(k|q|p) = E E r“(q|k|p). 

|k|<K lqi<A:,|p|>A: |k|<K |q|<AMp|>Ar 


(D.5) 


(D.6) 


(D.T) 


The sum E|k|<K EH<A',|ql>A' 5r“(k!p|q) involves transfers from k p and also from p k. 
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which are equal and opposite. Therefore 

E t lr‘(Mp|q) = 0. 

|k|<K lp|<K,|q|>A' ^ 


Similarly, 


and, 


E E = 0. 

|k|<K |p|>K,|q|<A' ^ 


E E 5r“(k|q|p) = 0. 

|k|<K |pl<K,|q|<A- “ 


(D.S) 


(D.9) 


(D.IO) 


In other words all terms except both legs of Fig. D.l, p — )• k legs of Fig. D.l, and q — )■ k 
legs of Fig. D.l survive. We use these results to rewrite n“>(A') as 

nS(A-) = E E r*(k|p|q) + E E , 5>?“(k|q|p) 


|k|<K |p|,|q|>A' 

^ 1 


|k|<K |q|<A-,|p|>A' 2 


|k|<K |p|<AMq|>A' 

r“(k|p|q) 


(D.ll) 


Using symmetry between p and q we get, 

E E lr’‘(k|q|p)= E E 5 r“(k|p|q) 


|k|<K |p|<K,|q|>A- 
Hence. we can write, 


(D.12) 


|k|<K lq|<K,|p|>A 


ns = E E i?““(k|p|q) + E E r"(k|p|q) 

|k|<K |p|,|q|>A' |kl<K |q|<AMp|<A 

= E E r“(k|p|q) 

|k|<K |p|>A' 


(D.13) 


This is the same as n“< of Eq. (3.28). Thus, we have shown that Kraichnan’s formula for 
flux is the same as our formula based on the mode-to-mode transfer. 
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